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[Text] Let us study the dynamics of a new type of correctable 
gyrocompasses--with two-race tuned rotor gyroscope (DNG) as the main 
sensitive element. The layout of the instrument is shown in the figure. 
The platform is stabilized with respect to axes Ox and 0z by 
electromagnetic torque converters (PM) from signals of the DNG angle 
sensors (IU). The “pendulum” and damping moments are generated by the 
electromagnetic torque converters of the gyroscope from signals of the 
accelerometer, whose axis of sensitivity is parallel to axis Oy. 
Correcting signals, which are generated by the computer of the control 
circuit from available information on the observer’s latitude, the speed 
of the ship and on the angle of deflection of the platform with respect 
to axis Oy arrive av the torque converters of the DNG to compensate for 
the velocity error and the error due to inclination of the platform 
together with the object with respect to the axis of rotation of the 
gyroscope. Axis Oy is measured by the accelerometer, whose axis of 
sensitivity is parallel to the horizontal axis of rotation of the 
platform. 


1. Mathematical model of instrument. The equations of motion of a 
one-rotor corrected gyrocompass are found by having joined the method of 
derivation of the equations of the DNG [1], follow-up gimbels [2] and of 
the control circuit [3]. 


Let us write the linearized equations of the DNG with regard to the 
symmetry of the flexible suspension of the DNG in the form 


(B +4- b)a, + (A + 20), -+ Acer, -} he, 4+ hOB, + 
+- [A + 2(a — @)]Qo,a, ++ (A — B)w,B, +(A— 2B)o,B, = (1) 
= M, —- (B - b)o, — (A 4-4 4-b— 0) 20, — (A + a — B— dw, 0: 











(B + b)B, —(A + 2b) Qa, + Ach, + hp, — hQa, + 
+ [4 -+ 2 (a — ¢)] 20,8, — (A — B) oye, — (A — 2B) 0,2, = 
= M, — (B-!- b) , L(A -f- a -} Ob — 0) Qo, -4- (A fa — Re C) OY, . 


where A and a are axial moments of inertia, B, b, and c are equatorial 
moments of inertia of the rotor and races of the DNG, ? is the angular 
rotational velocity of the drive shaft of the DNG, Ac is the residnal 
stiffness of the gimbal suspension of the DNG, h is the viscous damping 
factor of rotation of the rotor, ar and #, are the angles of rotation of 
the axis of the gyroscope rotor with respect. to the body, M, and My are 
projections of the total control], correction and noise torque, acting 
from the direction of the body to the gyroscope rotor, onto axes 0z and 
Ox, and wy, wy, and w, are projections of the angular velocity of the 
platform onto axes Oxyz. 


















































KEY: 
1. Heading 3. Torque converter 
2. Angle meter 4. Accelerometer module 


Assuming that the kinematics of the platform is similar to that of the 
follow-up sphere of the corrected gyrocompass with liquid-torsion bar 
suspension of the sensitive element [3] and taking into account that 
angle 7 is deflection of the platform together with the object with 











respect to axis Oy, we write the projections of its angular velocity in 
the following form: 


, == o, +- p. -{- O),%. — (a. -|- ;) y ’ (2) 
0), = a O,% - (wo, -b a. )B. -|- y ’ 


, ~~ , -| 1, _ o,f. ++ (B. + O,) y, 


where we ay? and “¢ are projections of the angular velocity of the 
object onto the axes of geographic coordinate system 0f7( [3] and a, and 
be are the angles of rotation of the azimuth and horizontal follow-up 
gimbals. 


The equations of motion of the follow-up gimbals are [2]: 


Ia. -+ ha, = M, + My. (3) 


Ip. t- h,B, = M,. -{- M, ° 


where /, = /, -|-B 46; 1, ,-| B46, 1,1, are the moments of inertia of 


the azimuth and horizontal follow-up gimbals, h, and hy, are the viscous 


friction coefficients in the axes of the follow-up gimbals, M,. (i = 
1 


= 1, 2) are the electromagnetic moments of the torque converter of the 
stabilization system, and Mg (q = (, €) are the torques of th> 
perturbing actions. 


The electromagnetic moments are determined by the expressions 
M,. = WV, (p) &,., M,. er Vv, (p) Bs (4) 


Where W,,(p) (i=1, 2) are the transfer functions of the relief circuit 


regulators. 


The laws of position control of the sensitive element of the gyrocompass 
are taken according to [3]: 


M, —_— nd _ M; + M7; M, =e n§ -|- M: -/ M ’ (5) 








where é is the output signal of the horizon indicator, ny and nz are the 


scaling coefficients, Mk, are correction moments, and M", (q= Zz, x) 
are the perturbing moments that cause drift of the gyroscope. 


Let the output signals of the accelerometer module (BA) be described by 
the following relations: 


og v 4 er 
ager ht). aaah 


where 6 and , are the signals of the north and east horizon indicators, 
respectively, T, and T, are the time constants, and We and ", are the 


east and north components of acceleration of the object. 


Equations (1)-(6) are the input mathematical model of a one-rotor 
corrected gyrocompass with DNG. 


It has become traditional to use the precession model of the device when 
studying the dynamics and accuracy of corrected gyrocompasses. However, 
this approach, when the DNG is the sensitive element of the gyrocompass, 
limits the study of the dynamic features of the device. Let us 
construct an "integrate“' model of a gyrocompass, which takes into 
account the dynamics of the sensitive element and of the fo] low-up 
gimbals, from which the compass and vibration models of the device, 
respectively, can be found as special cases of analysis of low- and 
high-frequency vibrations of the object. 


The absolute angles of rotation of the sensitive element with respect to 
the geographic system of coordinates are expressed in the first 
approximation by the following relations: 


A=Ge -}+a,, P=f -}- fi. | (7) 


ti 


Separating the relative angles of rotation of the gyroscope 4, and §, 
from the last expressions, let us substitute them into the equations of 
motion of the follow-up gimbals (4) and (3). Then 


%_(P) = Da, (P)% + De (P) My. B(P) = Mp, (P) B+, (p) My (8) 


a, (p) = a, (pr) — OD, (p) AL. BCP) = Mp, (P) B— 1, (P) ML, (9) 








where 


Dy. (p) = We, (p) Up? hp A Ws, (PI, 


(10) 
Wp. (p) = We, (p) [1,p? + hp -t- Ws, (p)\': 
Pe, (p) = | — My. (/?), (Da (p) = |] — We. (p); (11) 
Ve (p) = [,p? + hp We, (PI, Ye (p) = Lp? + hp + We, (p) 7! (12) 


Let us write the initial equations of the sensitive element (1) in 
Laplace transforms at zero initial conditions with regard to (9): 


(/s? -f hs -- Ac) (Pq, — Dee) 4- (11s +- AQ) (Ye, f — DM) = 
= M, —/sw, — II"o,, 

(/s* -4- hs -- Ac) (Mp. — DAL) — (1s -+ AQ) (Dg, — DM) = 
= M,--/so, + /1*o,, 


(13) 


where J =-B+-b; H=(A+2b)Q; Ii*=(A+a+b—c)Q; s=jo. 


The equations of motion (13) are the “integrated” mathematical model of 
a gyrocompass for a wide frequency spectrum. 


It is known that the error of reproducing the amp]itud2 in the bandwidth 
of a closed system does not exceed the permissible value, while the 

amplitude-frequency characteristic decreases with a further increase of 
frequency. The following relations are valid for these frequency bands: 


1)0<¢we< g 7 ~bandwidth : 


Do (w) © 1 => Dy, (ww) > 0, Dy (0) & Kr", 


- (14) 
Dp, (0) 1 => Vp, (w) > 0, Dy (0) & K's 
2) «> v,~outside bandwidth : 
Dy (0) > O=> Da, (w) & 1, Dy (o) +0, 
(15) 


Wp. (@) > 0 => Dp, (wo) ~ 1, OD (@)> 0, 
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where Kj (i = 1, 2) are the total amplification factors of the 
static-type relie*® circuit regulators. 


Using relations (14) and (15) and omitting the intermediate 
transformations, let us write equations (13) for two frequency bands. 
The equations of a gyrocompass in the bandwidth of a closed system 
assume the following form in temporary form 


H* (6 + @, +0,% —ay — oy) = M, + Kz! (AcMy + hQM,), ae) 


— H*(a +0, —0,8 + By + 0,1) = My + Kz" (Aci, — hOM)). 


For vibrations at frequencies that do not exceed the bandwidth of a 
closed stabilization system, equations (13) are transformed to the 
following system of differential equations: 


loo + hee + Aca. 4- HB + hOB = ME — Ju, —H*u,, 


116 4- hf + Ach — Ha —hQa = MP — Iu, + Hu. 


Here ux and uz are the vibration components of the angular velocities of 


the platform #, and w#, and Me and Me are the high-frequency components 
of the perturbation moments. 


Let us call the motion of the sensitive element with frequencies 
belonging to the bandwidth of the closed system compass motion, and let 
us call those outside the bandwidth vibration motion. The follow-up 
gimbals track the position of the sensitive element with accuracy up to 
static errors in the case of compass motion. The follow-up gimbals, due 
to inertia, are unable to counter rapidly variable vibrations of the 
sensitive element in the case of vibration motion. 


Separation of the "integrated" model of motion of the device into 
compass (16) and vibration (17) permits one to determine the dynamic 
characteristics of a corrected gyrocompass with DNG. 


2. Analysis of compass motion. The compass model of the device (16) 
assumes the following form with regard to the equations of moments (5) 


H*(B +, + 0, — ay —o,y) = — 1,51. ME 4 Me 


4- KT" (AcM, +- hQM,), 


— H*(% + 0 — 0,8 + By + oy) = — 1,6 + Mi + Mr + 
+ Kz! (AcM, — hQM,). 
T 








Based on conditions of providing the position of dynamic equilibrium in 
the direction of the true meridian and in the plane of the horizon, we 
find ideal values of the correction moments 


Ate" -- ]{* (, — ay —_— ,Y) — M; — Kk; (AcM, -+- hQM,), 18) 
Wy" -= --- HY (, -|- by -+- oy) — M— Ke" (AcM, _- hQM,). 
The following correction moments are realized in practice ina 
gyrocompass : 
ME = Ho, 1 2d) — M; -|- AM,, 
(19) 


My == — Hw, — My -+ AMx. 
where AM, and AM, are the errors of the correcting moments due to errors 
of measuring the latitude and velocity of the object. 


The imperfection of the correction moments resu]ts in a methodical error 
of the gyrocompass: 


a, -- U-'sec p| — @ — Ma -+- WS (6a, — 5p.) + T=! 6p, — mba,)| , (20) 


where U is the Earth’s angular rotationaql] velocity, » is the observer’s 
latitude, m, = M)/H*,«,=—A}/H* are the azimuth and horizontal drift 
of the gyroscope, 6,. = Al;4;', 63. = M,K;! are the static errors of the 
follow-up gimbals, 7. -- //*(j8)~! is the time constant of the gyroscope, 
and m = nz/n,. 


Let us calculate as an example the error component of the gyrocompass 
(20), caused by the time constant of the DNG and by the static errors of 


the stabilization system. Let us assume that T, = 50 s and by = bq, =0,1'; 


m=(,05, Vcosp=7-10-§ s-!, We find in this case that ao = 28.6’. 


Accordingly, disregarding such features as the time constant of the DNG 
and residual stiffness results in significant errors of the device. 
Therefore, it is recommended that the last terms from (18) be taken into 
account in the expressions of the correcting moments (19). 








3. Oscillating deviation of gyrocompass. Let us use the compass model 
of the device (16), having first discarded terms that contain static 
errors of the follow-up systems for simplification of computations. Let 


there be rolling motion ()=0,,sin (wf+e). The pendulosity of attachment 
of the center of gravity of the device 1 causes rotational acceleration, 
projections of which onto axes 0f and Oy are: W;=0/ cos K; W,= —Olsin K, 
where K is the heading of the object. 


Omitting the intermediate transformations, we find the formula of the 
intercardinal deviation of the gyrocompass 


‘ . » ra) 
(2) = n,l02 o* sin 2K IT. COS (f2 — Fy) —- 7 ecw (Fy — &) 








a ” . 9\9 ol | 
x {4n,gU cosp | 1 -{-m?7) | (1 — 7,7 0°)? +- wT; | (21) 
| -—-T, To? | 
re e, = arclg-—,7 *— 5 ef, == arclg or 


z 


e, = arctg(— o7 ,.); 1, = , 
2 


The differences of deviation (21) from the known deviation [4] are 
manifested due to the structural features of the platform suspension: a 
two-axis suspension required introduction of algorithmic compensation of 
the gyrocompass errors, caused by heeling inclinations of the object. 
However, the constant value of perturbation <é(y + 7)>, as a result of 
calculation of which formula (21) is found, is manifested due to the 
inertia of the accelerometer--of the indicator of given inclinations. 


Thus, the frequency method of separation of motion permitted us to find 
the compass and vibration mathematical models of a corrected 
gyrocompass, convenient for analysis and which reflect the specifics of 
application of a DNG as the sensitive element. Disregard of the time 
constant and of the residual stiffness of the DNG in the case of static 
regulators of follow-up systems results in considerable errors in 
determination of meridian, which should be taken into account when 
working out algorithms for the control and correction circuit. 


Study of rolling deviations of the device showed that algorithmic 
correction of the error of the gyrocompass due to rolling inclinations 
of the object and due to a two-axis layout of platform suspension does 
not correct the intercardinal deviation due to the inertia of the 
accelerometer-indicator of the given inclinations. 
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[Text] The purpose of this article is to develop a method of 
determining the geometric dimensions of flexible elements (Figure 1) of 
tuned rotor gyroscopes (DNG) upon fulfillment of the following 
relations: 


1) the condition of strength upon exposure to impact loads 


~ A,B,’ (1) 


where [7] is the permissible tangential stress, p is the permissible 
impact load per flexible element, and Aj and Bj are the geometric 
dimensions of the flexible element (Figure 1) of variable (i = 1) and 
constant cross-section (i = 2); 


2) the condition of strength upon exposure to cyclic loading 


fy 
lo] 25,=+— ? (2) 


where [7] is the fatigue strength, 7 is the permissible angle of 
rotation of the flexible element, c, is the given angular stiffness 


coefficient of the flexible element with respect to the working axis, W 
is the moment of cross-sectional resistance of the flexible element, and 


11 








Vay | (3) 


3) the condition of conservation of the angular stiffness factor in 
a specific range is 


Cymin < Cy S Cymax: 


Let us consider a flexible element of variable cross-section. For this 
element, 





— 26 A5/2p p-i/2 
Cy = oq Ar BR (4) 
or 
9 
By =— set Ar R", (5) 


where E is Young’s modulus. 























Figure 1 


Substituting (3) and (4) into (2), we find condition 2 in the form 


12 








4E - 
lo] > 3, AU R'"G. (6) 


If B; = Bo,;, let us derive from (1), (5), and (6) the expressions that 
give the range of selection of dimensions A; and R with regard to 
constraints on the geometric dimensions of the flexible element: 


R>K,,A, (condition 2), 


» 1GE%p? . 
Ku = Qniqoi 


(7) 


R> K Ai (condition /1), 








4E%p? . (8) 
Kis = Bini? — [e? 
ymin 
KyAR<R<K Ai (condition 3), 
4E*Bo, 
MO Bla max (9) 
Kis = ie ’ 
81M ey min 
A; © Aoy (constraint on dimension A; ); (10) 
R ¢ Ro (constraint on dimension R). (11) 


Let us consider a flexible element of constant cross-section. For this 


element, 
cy = gr ABs (12) 
or 
» - bin Bi (13) ‘ 


13 








Substituting (3) and (12) into (2), we find conditon 2 in the form 


E —_ 
[0] >= ALL '9. (14) 


If Bo = Boo, let us derive from (1), (13), amd (14) the expressions that 
give the range of selection of dimensions A» and L with regard to 
constraints on the geometric dimensions of the flexible element: 


L>KyAq (condition 2), 








‘ (15) 
( 
L>KyA; (condition 1), 
EP (16) 
Ky = 12¢ min [t] , 

KAr<L< Ky A3 (condition 3), 17) 

Ko, 2 oot, Ku = EB ) 

at 12c nan. 25 12C renin” 

Ao £ A o9 (constraint on dimension Ao); (18) 
L ¢ Lo (constraint on dimension L). (19) 


To find the specific values of the dimensions from the above regions, 
0, 
let us introduce coefficients Mo; = nn (i= 1,2), that take into account 


imino 


the different structural and technological constraints on the values of 
the desired dimensions and that correspond to flexible elements with 
variable (i = 1) and constant (i = 2) cross-section. 


Let us write with regard to B; = Bo; and Bo = Boo 


Ku) _ AL Buy AY’? Ro? A 


Bap 


‘ EV Boo Arp! 
Koo - =e Ad > 
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hence, 


no 9 Shi ’ » (20) 
yk ae Fai or = lon 
9, - 2~P . 
ASL ' sap “EqBor Ko» — los. (21) 


It is easy to show that 


? 2 1G E2Py? ,-4,— 
(O1T))min = OpminTimin = Qn? 4 Bo'R '. (22) 








SW als 
(92T2)imin = OominTomin = + Boo'L 7 (23) 
It follows from (22) and (23) that 
L = | (25) 


Substituting (24) and (25) into (20) and (21) and making the necessary 


transformations, we find \ 
\ 

A, _ Is,’R3/3: (26) 

A, = loy°Ly”. (27) 


When selecting the values of coefficients Ko; and Ko2, one must take 
into account that 


Pa ) 93/2 . - 4E@B —1/2 43/2, 
wo Ro? Ain < Kor Fun Ro “Aimax: (28) 
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EqB —i 4? n Ao 
Fi bet 15 Aimin < Kon < po La Admax (29) 


The method of selecting the geometric dimensions of the flexible 
elements of variable cross-section is based on relations (7)-(11), (24), 
(26), and (28), while the method of selecting flexible elements of 
constant cross-section is based on relations (15)-(19), (25), (27), and 
(29). 
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Figure 3 


To illustrate the foregoing method, let us consider a numerical example. 
Let E = 2:10!! N/m’, p = 9.8N, y = 1.74-10°2 rad, if] 2 46.3-10? N/m?, 

= 7 7 2 = .o° -2 ‘mM, = "« . 4 = 
[oe] 72.5°10/ N/m-, Cymin 0.8-10°4 N-m — 10°-- Nem, Boy Bo? 
= 2-10°3 m, Ao; = Aoo = 1074 m, Lo = 2-10°3 m, Ro = 1073 m, Koy = 10, 
and Koo = 0.5. 
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As a result of the calculation, we find kK;; = 4.2, Kj». = 1 
Kig = 0.94-103 m-?, Kiq = 0.08-102° m-4, Kis = 0.13-102° m 
Koa = 45-103 m'!, Ko3 = 36-103 m-!, Ko4 = 3.34-109 m-°, and Ko5 = 4.17 x 
x 109 m-2, 


Let us use the derived values to construct regions for selection of 
dimensions A;, R (Figure 2) and Ao, L (Figure 3). These regions are 
shaded. Curve 1 in Figure 2 corresponds to the expression R = K;;Ai;, 
curve 2 corresponds to R = K;4A5;, curve 3 corresponds to R = K;5A5,;, 
while curve 1 in Figure 3 corresponds to L = K2;A9, curve 2 corresponds 
to L = Ke4A%o, and curve 3 corresponds to L = Ko5A%9. We find in the 
regions for selection of the dimensions points M, (0.023; 1) (Figure 2) 
and Mo (0.034; 2) (Figure 3), whose coordinates are dimensions A;, R and 
Ao, L, respectively. 


The studies permit the following conclusions: if the overall dimensions 
of tuned rotor gyroscopes are identical (A; = Ao, By; = Bo, and 2Ro = L) 
and if the strength and stiffness conditions are fulfilled, the angular 
stiffness factor of the flexible element is less ai constant 
cross-section than it is at variable cross-section. It follows from 
this that the accuracy of a tuned rotor gyroscope with flexible elements 
of constant cross-section is higher. 
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[Text] The method of induced rotation of the gimbal suspension of a 
gyroscope about the vector of the moment of momentum [1] is used for 
autocompensation of the instrument errors of free gyroscopes. The 
purpose of this article is study of the possibility of applying this 
method to two-degree gyrotachometers (GT). 


Let us link the orthogonal coordinate system £7(, axes 7 and ( of which 
are directed in the initial position of the gyrotachometer along the 
precession axis and the principal axis of the gyroscope, respectively 
(Figure 1). Let us link system xyz to the movable part of the 
gyrotachometer such that axes y and z will coincide with the suspension 
axis of the gyroassembly and the principal axis of the gyroscope, 
respectively. The body of the gyrotachometer rotates uniformly at 
angular velocity # about axis ( with respect to coordinate system (7( 
(Figure 2) and a = {It. Let us assume for analysis of the dynamics of 
this gyrotachometer that the base rotates about axis 7 at constant 
angular velocity @. 


The linearized differential equation of motion of the gyrotachometer has 
the following form with the respect to the precession axis 


IB + fh + 1, —1,) (Q? — o sin? Qh) + HQ + c]p = 
= Hw sin Qt + (1, -+ ly —/,) Qo sin QF. (1) 


Here Ix, Iy, and Iz are the moments of inertia of the movable part of 


the gyrotachometer with respect to axes x, y, and z, H = IP, is the 
moment of momentum of the gyroscope, IP, is the axial moment of inertia 
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of the gyroscope rotor, y is the angular rotational velocity of the 


gyroscope rotor (7 >), f is the damping factor, c is the angular 
stiffness of the flexible element of the gyrotachometer, and j is the 
angle of rotation of the gyroassembly about the precession axis y. 











Figure 2 
Let us study the dynamics of the gyrotachometer in three cases: 


Q2>>0?; Q°’<w’; and} = vw. In the first case, let us write differential 


equation (1) in the form 


B + 20,8 + wf = Aw sin Qf, (2) 


where tu, =f/ly; ag = [c -+ HQ + (1, —1,) 22) 17'; A = (H+ ty Ag) UZ". 


The partial solution of equation (2) is as follows: 
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p = = === sin (Q/ — 9), 
V/ (a) — 99)? +4. afteo9Q# (3) 








where 
p = arctg 26w,82/(«p — Q2). 


It follows from (3) that the considered device behaves like a vibratory 
gyroscope, since it responds to constant angular velocity of the base # 
by vibration with frequency {), equal to the angular rotational velocity 
of the gimbal suspension. One can use the resonant operating mode to 
improve the transfer factor of the device. Resonance begins at wy = Hf] 
Expression (3) then assumes the form 


‘ \ 


_ A@ . m\ Ho \ @ 
p= 2h woh _ (t— +} a —(> Fl,+ ty — Ig) cos M1. (4) 


\ 


Since 7 } fl, regardless of the damping factor, which is selected on the 
basis of the required nature of the transition process, one can provide 
the necessary value of the transfer factor of the gyrotachometer by 
selecting the ratio between H and {? upon measurement of small angular 
velocities w. The considered device differs advantageously by this from 
rotary vibratory gyroscopes, since its transfer factor and phase shift 7 
become less critical to small variations of #o or Il. 


The relation Noes = wo is fulfilled for the resonant mode of the 


gyrotachometer, i.e., 
Qhes = [e+ Mpes 1 (Ly — 1g) QR TZ" (5) 
Hence, we find at I, = Iy the condition of dynamic tuning 


H + [H? + 1¢(21.—] yi"? 
pineal. = 1 i (6) 
= 7’ 





If there is no flexible element in the gyrotachometer (c = 0), formula 
(6) is simplified considerably: Qh. =U; Qh». = 1 (2/,—/,)-'.. Formulas 


derived for a rotary vibratory gyroscope follow from expressions (5) and 
(6) as a special case at H = 0 [2, 3]. 
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When the device operates in the second frequency band ({)? ¢ #?), 
equation (1) assumes the form 


If +}. /f}-4- le -}- WQ--(, --- 1,) 0? sin? Q2] B = Aw sin Q/, (7) 


and ir the third case ({} = «#) 


Hf 1 1B | Ie} H&-+ Ul, —1,)0* cos? Qt B = Ao sin Q. 8) 


Using the trigonometric identities, equations (7) and (8) can be written 
in typical uniform form 


fi Quah [fe - geos 20f] (== Aw sin Q1, (9) 


where ow. - fof HQ 4 Lye 2] hs Yo 77 (1, —- 1.) w*/21, if N? ¢€ #°, and 

Uy do ty le | IQ, -1)w/2y, if N= @. 

The dynamic stability of the gyrotachometer, the motion of which is 
described by linear differential equation (9) with periodic coefficient, 


is provided by selection of the damping factor f (or of the relative 
damping factor ( = f/2Iy#o) according to the inequality [4] 


C2 L8LL ++ ging — (1 + 2qun*)'”). (10) 


Calculation of ( by formula (10) shows that the dynamic stability of the 
considered device is easily provided by the parameters of serial 
gyrotachometers over a wide range of angular velocities w, since quv°-» ¢ 
<i. 


Assuming that inequality (10) is fulfilled, let us represent the partial 
T-periodic solution of equation (9) in the form of a segment of the 
Fourier series (T = 7/f) [4]: 


[} — a, +- a, sin Qt +- b, cos QL -| ay sin 20¢ +. be cos 2024, oe 


Having substituted (11) into equation (9) and having set the 

coefficients at identical harmonics equal, we find a, -=a,-=b, = 0; a, = AAW; 
by = AAT A, = A(om—Q%;  Ag= — 4% %WQ—g/2); A=(H— 

+)? — 9/2) (mp — Q2) -}- 2h, Q (2b, 82 — gi 
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Thus, the output signal of the device can be described at fl? € #? and 
f. = w by the following mathematical expression 


f = Dw sin (82 +- e), (12) 


where |) = A(a; +0j)'"; €= arclg If the resonant operating mode of the 


gyrotachometer is used (#9 = Mies) formula (12) yields formula (4). 


Thus, a two-degree gyrotachometer with rotating gimbal suspension has 
properties of a one-cage rotary vibratory gyroscope (RVG) [2, 3]. The 
advantage of the considered device compared to a rotary vibratory 
gyroscope, besides those noted earlier, is the possibility of realizing 
a higher natural vibration frequency ¢), dependent on additional dynamic 


components of total stiffness cy=c 4 //©-\- (/,—/z) (Q?—w? sin*Qt), 
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[Text] Numerous studies of a gyroscope with flexible suspension | 1-5] 
have not devoted attention to the acceleration of its rotor, which is 

important in solving the problem of reducing the readiness time of the 
devices. 


A method is proposed in this paper and the characteristic features of a 
gyroscope upon acceleration of its rotor are studied. 


Analysis of the nonlinear equations of motion, which correspond to 
different kinematic layouts of gyroscopes with flexible suspension 
(layouts in the form of a rotor with one- and two-race and inertialless 
flexible suspensions), permits one [5, 6] to write the generalized 
mathematical model of a gyroscope with flexible suspension in the form 


VA OW + (bay + by) P= eM, ™ 


O — boapip + (Dasy? -I- Og4) 0 = eM: 


Y - wi Ay + byo'? = Elll gy, 
(2) 


” + 7 + by _ eMs, 


es) 
G 








where coefficients bj; are dependent on the inertial and elastic 
properties of the system. 


The differences of systems of differential equations (1) and (2) for 
different kinematic layouts of gyroscopes with flexible suspension are 
reflected in the specific type of coefficients bjj. 


As follows from system of equations (1) and (2), the motion of the rotor 
with respect to coordinates 7 and y is independent in linear 
approximation of motion with respect to the remaining coordinates. 

Thus, one can study the problem of acceleration of a gyroscope with 
flexible suspension by equations (2), while one can study the motion of 
the rotor with respect to the drive motor shaft by equations (1). 


Let us assume that coefficient b30, which characterizes the stiffness of 
the flexible suspension, is large and the inertial properties of the 
rotor can be disregarded with respect to angle y, i.e., one can assume 


that y = 0 in equations of motion (2). 


The last equation of systems (2) is reduced by replacement of variable 
y= 2, 7 = 2 to an equation with separable variables 2 n2? = Mg 
Integrating this equation at initial conditions 7(0) = 0, 7(0) = 0, we 
find 


meg CEE 1’ (3) 


where L = MgMc 0» 


Having substituted (3) into the first equation of system (2), we find 


Bing L et lt 
bao (e2Lt -|- 1)? : (4) 


(p= — 





Let us estimate the effect of the inertia of the rotor during its motion 
along coordinate y on the law of variation of angular velocity (3) and 
angle (4). 


Since there are no methods of constructing the general solution of 
system of nonlinear equations (2), they can be integrated numerically on 
the digital computer. Let us solve this system by the fourth-order 
Runge-Kutta method, reducing to the form 


2h 











Y = (Ms (Dyg — 1) — mg ++ meop® — byobgeP)/(Yyg ~ 1), 


P = (ig — meoy* + bsoQ)(bqg — 1). 
After substitution of variables y»—=x,, p=%, Y= 4%; P= 4%, we write 
Y =X. P > Xgy Xp = Xs, x,-- x, Finally, the canonical form of system (2) 


assumes the form 


x, == Xz, Xe = NX, 
X = (€Ms (4, — 1)— Mg + MoX3 — b59ga%o)/(Oq, — l), 


X= (Mg — Meo x3 + b392)/(b42 — |). 


(6) 


Let us assume that the engine torque is constant. The value of mg = mg) 


on the first time interval (7 < 0.870), and mg = Mg2 on the second 


TCK 
time interval (0.870 < ar 


System of equations (6) is integrated numerically on the SM-4 digital 
computer with values of the parameters of 70 = 1,500 s~'!, bya = 1 Nem, 
bso = 10-8 N-m, Moo = 20s-*, mg = 450 N-m and at initial conditions 
7(0) = 0, 7(0) = 0. 

The results, found analytically without regard to elastic deformations 
of the rotor with respect to the shaft axis by angle » and by the 


numerical method with regard to the flexible properties of the rotor, 
coincide. 


It is obvious from analysis of these results that the consideration of 
the inertia of the rotor results in additional vibrations along 
coordinate y. This must be taken into account when designing the 
suspension and in strength calculation of it. 

It frequently becomes necessary in practice to force acceleration of the 
rotor within a limited time ty: It is desirable that 7 approach 40 
smoothly (without jumps and vibrations), i.e., that it correspond to law 
(3). 


Function (3) can be approximated by a function of type 


y= rol 6) (7) 


2) 











so that at 


(8) 


We find with respect to the law of variation of the angular rotational 
velocity of the engine shaft (7) with regard to (8) from system of 


equations (2) 





0 aif ot” (9) 
m., — "2 —at Yo%"b42 — 9 42 ‘9 at 
gq = Meofy — @ ator Vv -|- 2MenVy — MeoyVee . (10) 


The graphs of function m,(t) at different values of ty are shown in the 


figure, where curve (1) corresponds to ty = 1 s, curve 2 corresponds to 


= 0.5 s, and curve 3 corresponds to t. = 0.1 s. 
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It is also of interest to study the motion of the rotor with respect to 
angular coordinates @ and #, which characterize its dynamic properties 
as a gage of navigation instruments. Let us first study the motion of 
the gyroscope in linear approximation. 


Let us rewrite system of equations (1) and (2) with regard to (7) in the 
form 


26 








P + Oisyo0 -+ [D,s7? + by)p = braver “0 — bisy? (—2e~™! + oy — 


— bisYo, {cos (yf) cos (te e™') — sin (Yo) sin ( Js *)| n 


\ 


+ bates [osu cos (Be) — sin sin 28 4] 


+- 2h3 | — ip a 00 ” voto |; 
a _ ; . . A (11) 
0 -~ Desyn'P {-[Posy? -- by, |9 = — Dose p —_ Desy? (—2e—7 -|- e204) 4 


-|- Dare, sin (yf) cos (te *) -|-cos (yl) sin (te o)| —bs,ye-™o, x 


a 
X sin (yo!) cos (1. e) +- cos (yo!) sin (4 “| ae 


— 2h, [0 a yo"? — Yor 2p]. 


The considered gyroscope based on a flexible suspension is a vibratory 
system. As studies show, vibrations of the rotor with respect to the 
drive motor shaft occur with period that is considerably less than the 
acceleration time by coordinate 7. Therefore, complex motion of the 
rotor, consisting of rotation together with the shaft and deviation from 
it, can be divided into "fast" (vibrations with respect to the shaft) 
and "slow" (rotation with respect to the shaft, including acceleration), 
and one can use N. N. Bogolyubov’s method of perturbations [7] to study 
the motion of the rotor upon acceleration, performing an averaging 
operation with respect to the "fast" time. The terms in system of 


equations (11), which contain the multiplier aaa (k = 1, 2) that 


characterizes the "slow" motion (acceleration), are assumed constant 
upon averaging, i.e., they are "frozen." Thus, we find a system of 
equations with coefficients, quasi-constant with respect to "fast" 
motions, in the acceleration phase for study of the motion of the rotor 
with respect to @ and ¥. 


Having used the method of averaging [8], let us find the solutions of 
equations (11) in the form 


2 
p -_ y Cer" 4- Dei", 
I=! (12) 


® 
2 


0 «xs v xy (Ce! — Dye~**1"'), 
j= 1 
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Having made the necessary transformations and having completed the 
averaging procedure, we find a system of equations of first 
approximation 


c 


Gyo Qa, (Cy {2 ype (Mtg 4- Mtg) - - i Ly Oye | Mga) -| 


+ y3 (— Qe A e 2b, -|- No.)}} -|- + yoo, (1 — m) “ 





_(% — 1. -_ “ 
x Sin (18 ¢ “\b1s -|- bg,) — > yy, (l —e ‘) COS (Teena! (Py5+-0o5)). (13) 
, e o. ar 
Dy = — gy Dy {2poe—* (Ng thts) + i Lyje X (bye + ban) +t 


-|- y; (— Qe—=! -|- ge") (M15 + by.)J} + 5-10, (i—-e “) x 


(Ye | | } 
x sin (ae e at (O45 -1- Og;) -f- “5 Yo, (1 — e*) cos (2 Tame (UE } M55)). 


The equations for variables Cy» and D» describe the attenuating 
vibrations of the system with second natural frequency and are of no 
practical interest in the given postulation of the problem [5]. 


As follows from analysis of system (13), the variable coefficients are 
determined by variation of the angular rotational velocity of the motor 
shaft during acceleration. Therefore, the solution of equations (13) 
can be divided into two steps: I--corresponding to acceleration of the 
motor shaft when system of equations (13) has variable coefficients and 
II--corresponding to motion of a system with angular natural rotational 


velocity close to constant (7>¢ = const). Equations (13) will have 
constant coefficients. 


The general solution of system (13) must be found to study the motion of 
the gyroscope in the transient mode. Since the equations of the system 
are unrelated and since C; and D; are complex conjugate values, it is 
sufficient to construct the solution of one of the equations of the 
system, for example, of the first equation. 


Let us rewrite it in the form 
Cre Cyl) -b ie) 1) — ig (0, (14) 


where 
a 
WU) Volta} tisde~ 


’ 
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& : - .? { . — at ) ‘ 
ho=—- 7 yoo, (1 -—— e~*) sin { . (P15 “be Pas) 


, 
j 


g (i) = — a Yo, (1 — e~*') cos G oe (D.5 -f Og5). 


Equation (14) is a first-order inhomogeneous linear equation with 


variable coefficients. Solving it by the method of variation of an 


arbitrary constant, we find in the first step 


C,=L+ iM, 


where 


& 


(I — Kyl + K) (Ky (cosa +- cos) + K, (sina — sin B) + 
+ K; (cos y + cos ) + Kg (cos y — cos p) -+ K (sin y— sing) — 
— K, (sin yp 4- sin @)] +- > VA,; 

M = + (1 — K,t— K,) [K; (sin +- sinB) + K, (cosB — cosa) + 
+ Ks (sin y + sin @) + Kg (sin y — sin @) ++ K; (cos =p -|- cos y) + 

+ Kg (cosy -+ cos q)] >A dy, 


and #, #, 7, and y are time functions, while K, (i = 1, 7) are 


(15) 


(16) 


coefficients that are dependent on the parameters of system (14) (they 
are not presented due to cumbersomeness), and A, and A» are sums that 


contain teims of higher orders of smallness with respect to t. 


Let us find from (15) the conditions on the second step (t ? t.) at t 


p 
= ty: C; = L* + iM*, 


The solution of equation (14) on the second step has the form 


& ; 
= Sta FES aM 





_ , is 
Ce |- (pw -{- 4) 4+- ——"*—_, 
Si9 -}- tSi4 
where 
$418 SyoS 

p= LF —- ee; oy = Mt —**. 
| 9 9 w+ 

Si, + Sj2 it Sho 


(17) 








Let us derive from (14) and (17) the law of motion of the rotor over the 
entire time interval from the moment of acceleration until it reaches 


the steady mode 


€ . ; 
c a (t < t,) 4 iM (f < t.) 4-¢ ? (Syeté ut—ty) 


+ i(M*— frase )]+ iSi6Si2 


2 2 - 
aT v Si9 Si + Si2 


Taking relations (12) into account, we find from equation 


; —_ — Sis(/—fy) & ] 
YP = 2C0S yYol x <t)+e {cos| —-> Su ll— ) x 


Lt Susie | +-sin] ~ su (¢—1,)|[M*— rusts, |} + 


2. 2 9 9 
$i) F Si2 Si; + Si9 





. “= Six(!—fp) : 
+- 2sin pot M(t<t) +e {cos |— + sult—#,)] x 





2 

x [Mr $1951 |—sin|— su (t— )|[4+— -|} 4. 
Si + Si0 . . sit + Si2 
Si6Si2 
sit + Sto 
e 

o. “> Sis(f—-fy) £ 

0 = —2sin yet |L(t<t,) +¢ {cos|— 5 su t—t,)] x 
*  _ $11516 | gin |_E _ qx _ __ 512516 __ || 

x{L s\, i | = | -™ ( A) [ar si ss Sto \\ T 


e 
— > Sil!—ty) 


+ 2.008 Yof [Mi ((<t)+e ? cos |— + su(t— ,)| x 


x [ate — tte) — sin ff 54) fr si |} + tee 


2 1g 
Si +5t9 


2 2 2 2 
Si, + Sio Si) + Si9 


(18) 


| 


Thus, the law of variation of angular velocity (7) can be taken as the 
initial law in determination of the required torque of the motor in the 


forced acceleration problem (10). 


The motion of the rotor with respect to coordinates # and @ is a 
transient vibration process that results in increased dynamic effects on 


the flexible elements of the suspension. This nature of motion 


should 


be taken into account when guaranteeing the strength of the suspension. 
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[Text] The dynamics and errors of vibratory gyroscopes (VG) with 
flexible suspension as a sensor of small angular deviations of the base 
or as an angular-rate sensor in the steady mode were studied in 
sufficient detail [1-3] in the case of a free rotor of vibratory 
gyroscopes. At the same time, problems of design of the angular-rate 
sensor (IUS) on a tuned rotor gyroscope (DNG) with the required speed 
[4] have hardly been studied. A model of the errors of the angular-rate 
sensors on tuned rotor gyroscopes with symmetrical flexible suspension 
is developed in this article. 


As shown in [4], a two-race DNG can serve as a sensor of small angular 
deviations of the base in the transient mode of the gyroscope, and as an 
angular-rate sensor in the steady mode. Errors caused by friction 
losses, asymmetry of the dynamic parameters, angular and translational 
accelerations and vibrations of the base, and also by cross-coupling due 
to inaccurate dynamic tuning, which as an effect only during an initial 
error distinct from zero, are inherent to DNG in both cases. The speed 
of the DNG reaches several tens of seconds [4], which does not satisfy 
current requirements on the angular-rate sensor. Compensation feedback 
(OS) with respect to moment, by selection of the transfer function of 
which one can increase both the speed and can reduce the cross-coupling 
of the IUS, can be introduced. 


The equations of motion of the tuned rotor gyroscope can be written in 
complex form 


ep -+ (¢, -I- iN) + (C+ ie,y) = —e,Q— §@682 — Ina Kun + My, 


re p=a+ ip, Q=2,+ iQ); e, =B,+C,; ee=K3+Ky& = om 
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= (Ag+ A, +C,—Byys MH =(As+2C)y; C= 255 + (A —B, 


—C))y Mo =M,+iM ®=—1; a, B 


are angles of deflection of the rotor of the tuned rotor gyroscope in a 
coordinate system unbound to the drive shaft, fy and 2, are projections 
of the vector of the angular rotational velocity of the base onto the 
axes of sensitivity of the tuned rotor gyroscope, Aj, Bj, and Cj are 
axial and equatorial moments of inertia of the center rings (j = 1, 2) 
of the flexible suspension (UP) and of the rotor (j = 3), kj are 
coefficients of the viscous frictional moments of the center rings 


(j = 1, 2) and of the rotor (j = 3), y is the angular rotational 
velocity of the drive motor shaft, Css is the angular stiffness of the 
flexible axis of the tuned rotor gyroscope, M; and Mo are the perturbing 
moments of the angular-rate sensor with respect to the axes of 
sensitivity, caused by parameters of both the tuned rotor gyroscope and 


of the elements of the feedback circuit, Tan is the control current in 


the torque motor (MD) circuit, the functional dependence of which on 
time t and on the parameters of the tuned rotor gyroscope is determined 
by the selected feedback circuit, and Kun is the characteristic slope of 


the torque motor. 


Taking into account that the output signal in the compensation operating 
mode of the angular-rate sensor is the current in the torque motor 
circuit, we find from expression (1) 


Tun = [— iegQ + €,2 + My + ep — (¢y + iH p —(C + ies}! (2) 


Hence, it follows that the angular drift speed {1} is measured by the 
angular-rate sensor with some error, determined both by the dynamics of 
the tuned rotor gyroscope and by moments Mj, determined by the 
characteristics of the angular-rate sensor. Let us compile a 
mathematical model of the errors of the angular-rate sensor in the 
steady operating mode, since the dynamic errors of the sensor are 
considerably dependent on the selected feedback circuit. To do this, 
let us introduce the following notations: mj is the mass of the center 
rings (j = 1, 2), of the rotor (j - 3) of the tuned rotor gyroscope and 


of the drive shaft (jj = 0), rej= (Teivst Teini rej} is the displacement 


vector of the centers of mass of the center rings (j = 1, 2) and of the 
rotor (j = 3), rjox is the nonintersection of the flexible suspension 


axes, and W={W,; Wy; W) is the vector of translational acceleration 
of the base. 
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The useful signal of the angular-rate sensor that characterizes the 
angular drift speed of the base is selected in the form lo== —(ieg/Ky 0 


The remaining elements of the right side of expression (2) are then an 
error of the sensor, which in the steady mode has the form 


I, = [M,—(C + iesy) Pi Kua - (3) 


Let us consider the components of moment Mo, which determine the errors 
of the angular-rate sensor. 


The moments that occur in the flexible suspension of the angular-rate 
sensor during translational accelerations of the base (Wy = const and 
W, = const) can be written in the form 


My = — (Msfeaxg 4- 0,50 re 44) W — 0,5 (tary g, - O,5rtgt ea aN, 
Ww - Ww, + a). (4) 


They characterize the error of the sensor, dependent on displacement of 
the centers of mass and nonintersection of the axes of the flexible 
suspension in the axial direction. 


The radial unbalance of the center rings and rotor of the tuned rotor 
gyroscope determine the errors of the angular-rate sensor in the case of 


action on the translational vibration sensor of the base We Wee cos mf: 
WV,= W 0 COS Nol; W,= Wa cos Nol at frequency relations =; No=N3=2y. 


The perturbing moments can be written as follows after averaging on a 
time interval that considerably exceeds the period of natural rotation 
of the drive shaft: 


(5) 


My Gy = 9,5 (MgFeaxa “+ 1 Ferys)W x0; 


Mey) = 9,25 (Mar yoe Morea 1 Mafcaxe) (Wyq + EW 0). (6) 


As follows from expressions (4)-(6), the errors of the angular-rate 
sensor are dependent on the modulus and mutual orientation of vectors 
rej (Jj = 1, 2, 3), which characterize the errors in manufacture and 
assembly of the sensitive element. The center rings, the mass of which 
is comparable to that of the rotor of the tuned rotor gyroscope, have a 
significant influence on the errors of the small angular-rate sensor. 
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It becomes necessary in this regard not only to adjust the races 
dynamically to ensure the resonant operating mode of the tuned rotor 
gyroscope, but also to balance them statically. 


A quadrature moment Myn? the cause of which is correlation between 


displacements of the rotor upon elastic defcrmations, determined by 
asymmetry of the spatial arrangement of the pliable elements of the 
suspension, by rotations of their principal stiffness axes, by different 
dimensions, by violation of spatial symmetry, and by instrument errors 
in manufacture and assembly of the suspension, also occurs under the 
action of translational accelerations and vibrations of the base in the 
angular-rate sensor. Studies show that tne effect of instrument errors 
in the quadrature moment of the angular-rate sensor increases 
considerably (tens of times) in the presence of pliable center rings. 


Having denoted the stiffnes matrix of the flexible suspension by lIci ji 
where c,j (j = 1, 2, 3) are coefficients of the translational] 
stiffnesses of the flexible suspension and Cj; (j = 4, 5, 6) are the 
coefficients of angular stiffnesses of the flexible suspension, let us 
write the components of the quadrature m ment. 


The moment caused by cross-couplings between the translationa] and 
angular displacements of the rotor are 


Myoy = 0.513 (C9 + Cog) (reayaW yCa2! + ires,3W C33). (7) 


It is characterized by stiffness coefficients C5» and Cg3 and by radial 
displacements of the center of mass of the rotor of the tuned rotor 


gyroscope. 


The moment caused by correlation between the translational displacements 
of the rotor in the axial and radial directions is: 


1 | Cys Cie Co3 Cos 


91 oe leans — 9G % e323 — ES — Ti” — 
ie a ic er (8) 








. 42 , 
-- 2.58 maY"resss | m,(Wy + iW) 
22 Coo — Ay? my 


Expression (8) determines the total effect of the axial and radial 
unbalance of the rotor and of linear accelerations in these directions 
on the error of the angular-rate sensor. 


The moment caused by cross-couplings between the angular displacements 
of the rotor is: 
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Ay C5, — Cos . 7 ‘gi 1. 
Mya3 = — wih (Fg (iMtg -f- 3) -1- M3 e32g) Wy +E (ryplttty (9) 
-}-M4) -|- Mafegna)W,), 


where ryo and rzo is the displacement of the point of suspension of the 
rotor of the axis of symmetry of the shaft with respect to the 
rotational axis. 


Moment (9), characterized by stiffness coefficients C54 and Css, occurs 
upon displacement of the center of mass of the system with respect to 
the rotational axis of the drive shaft. 


The moment caused by dynamic unbalance of the rotor by and by is: 


Myxss = 0,5A3 (6, -4- 5.) [(r.9 (tty A Mg) bE Leaggtta) Fy Ae trys (My + ty) - 110) 


-{- Fegygttts)W | (Ag --Ag)7!. 


It acts along the suspension axes of the tuned rotor gyroscope at 
nonperpendicularity of the axes of the flexible suspension to the 
rotational axis of the shaft. The noncoincidence of the flexible axes 
with the rotational plane of the drive shaft is caused by dynamic 
unbalance of the rotor or by instrument errors. 


Analysis of quadradure moment j,.,— 37 Muu: shows that the stiffness of 
pol 
the pliable elements of the center rings must be exceeded by an order or 
more with respect to the corresponding stiffnesses of the flexible 
elements of the suspension, when the center rings can be considered as 
absolutely rigid, to eliminate it. The arrangement of the flexible 
elements of each center ring should satisfy the symmetry with respect to 
the rotational axis of the drive shaft, while the entire flexible 
suspension, after rotation by 90°, should be identical to some 
suspension, symmetrical with respect to the plane perpendicular to the 
suspension axis. 


Besides the above moments, a flexible element proportional to the 
product of rotor displacements with respect to two mutually 
perpendicular axes acts on the angular-rate sensor due to inertial 
loads. Taking into account that displacement of the center of mass of 
the flexible system is directly proportional to the effective 
acceleration, the perturbing moment in the body axes system can be 
written as follows: 
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te Cy / (11) 


where m; is the reduced mass of the rotor of the tuned rotor gyroscope. 


As follows from expression (11), the translational stiffnesses of the 
flexible suspension in the axial and radial directions must coincide to 
eliminate moment My;. This can be provided by rotating the principal 
stiffness axes of the pliable elements of the flexible suspension. 


Besides translational accelerations and vibrations of the base, angular 
vibrations with double rotational frequency of the drive ,—Q cos Qyt; 
Q, --osin2\/, also affect the errors of the angular-rate sensor. The 


perturbing moment caused by them characterizes the nonidentity of the 
center rings of the flexible suspension and its systematic component has 
the form 


Moor = 0,5 (A, — B, + 3C, — Ag + By — 3C,)yQp. (12) 


According to expression (12), there is no drift upon angular vibration 


of the base at frequency 27 in an angular-rate sensor on a two-race 
tuned rotor gyroscope with parameters A; = Ao, B; = Bo, and C; = Cy». 
Thus, the symmetry of the inertial characteristics of the suspension of 
the sensitive element with respect to its axes makes the sensor 
insensitive to angular vibrations, the vector of which lies in a plane 
perpendicular to the rotational axis of the shaft. 


The residual stiffness (C # 0) and the viscous friction of the gaseous 
medium of the sensor make a considerable contribution to the errors of 
the angular-rate sensor. They cause a perturbing moment (see (3)) only 
upon deflection of the rotor from the zero position by angle ¢: 


M. = [2Cyp + (A, — By —C,) y? + iegyl. (13) 


Angle y is a consequence of the action of the torque motor on the rotor 
of the tuned rotor gyroscope in a compensating angular-rate sensor in 
the absence of translational rotation of the base. The current 
occurring in the circuit of the torque motor is determined by the error 
yo between the angular deflection sensor (IU) of the sensitive element 
and the rotor of the tuned rotor gyroscope in its initial undeflected 
position, the noise Vo of the angular deflection sensor, and Uo of the 
electric circuit of the feedback circuit. Substituting these values 
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into (13) with regard to application of signals Vo and Up to the angular 
deflections of the rotor .=V5 K~'\y; Up=L~'|Uo/(Woe(p) Kny)], where Kuy 


is the steepness of the angular deflection sensor, L™! is an operator of 
the inverse Laplace transform, and Woe (P) is the transfer function of 


the feedback circuit, we write 
M, = (2C yg + (A, — B, — Cy) y" + iesy) (Py + Ug + Vo): (14) 


This expression characterizes the errors of the angular-rate sensor, 
occurring in the presence of a feedback circuit. They can be reduced by 
selecting the elements of the feedback circuit and of the transfer 
function of the feedback. 


Thus, expressions (4)-(14) permit one to estimate the error of the 
angular-rate sensor if there are different perturbing factors. when 
they are substituted into expression (3), we find the mathematical 
models of the errors of the angular-rate sensor: 


Ly == 0,5 [2 (MgFesng 4- 9,50 Fey x) (WW, + Wy) + (Marien + 0,51tel cone) X 
x (W, + iWy) +- (Myesys +- My Fery1)W xo + 0,5 (sl 12% ¢ Feyxi + 
+ MoFere2) (Wyo +- Wo) + my ess Ess WyFeays + (eazs!ts cH Efe, + 


Ao. a Cra — Coa (( 


I “Ayd-As (Ca P49 (1g -F 1g) Ma eaga) Wey EP yg (My Le ") -1- 


5, 
+ myfeay:) W,) + var 2 (Fao (Ma 15) Fess MQW y i (Po io + 





+ M3) +- Teay3!tt,)W,) 4- (32 


. 
Cos 





: 2y7 
X Fie — 2-68 _Mat*resa3_ + ) my(Wy + iW,) +2 mvs (‘as ou Wr, +. 


. C33 — . 
+ joe Sn w,) — (A, — B, + 3C, — A, + B, — 3C,)yQq 4- 2 (2Cg 4- 


+ (A, ~~ B, —C,)y tT ieyy) (To +- V Mer I. Ug)) Kua + 8/, 


where 6I is the component of the output signal fo the sensor, dependent 
on the random variations of the parameters of the sensitive element and 


on the feedback circuit. 


This mathematical model permits one to estimate the error of the 
angular-rate sensor upon exposure to different factors, and to formulate 
according to given conditions the requirements on the static and dynamic 
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unbalance of the rotor and center rings of the sensitive element, to the 
elastic-mass characteristics of the suspension of a tuned rotor 
gyroscope, and to the parameters of the feedback circuit. 
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[Text] This article is devoted to description and study of the constant 
error component of a compensation pendulous accelometer with flexible 
suspension (MKA s UP), operating under spatial vibration conditions. 

The problems similar to that postulated have until now been solved 
without refining the structure of the control circuit of the 
compensation penduluous accelerometer and the final result of the study 
was an expression that describes the constant error component of the MKA 
under two-component vibration conditions [1]. 


The purpose of this paper is to find expressions that permit one to 
estimate numerically the constant error component of the MKA and also to 
determine the requirements on the parameters of the accelerometer, 
fulfillment of which permits one to prevent the appearance of this 
constant component, i.e., to make the MKA invariant to spatial 
vibrations. Let us turn to the matrix form of writing the equations of 
motion of a MKA with flexible suspension: 


(Mp? -| C)q = —mQ(t)—u—R(O, (1) 


where M and C are matrices that characterized the inertial and elastic 
properties of the MKA, respectively, q is the column vector of the 
generalized coordinates, m is the mass of the movable part of the MKA, p 
is a differentiation operator, Q = HW, (t) is a matrix that characterizes 


the influence of the horizontal component of acceleration, acting along 
the axis of sensitivity, on the MKA, u = U(p)q = DW(t)q is a matrix 
whose elements are determined by the action of the vertical and lateral 


components of acceleration (W(1)?== (0, W(t), W(t). 
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Assuming that all the elements of R(t) are small, let us find the 
solution of (1) by the sequential approximation method. The equation of 


the first approximation is (Mp?-+4-C-+U(p)) qos --m1Wy (4). Its solution 


is 
qe) ss —- mil (Mp? 4-C 4-U (py W, (/). (2) 


The equation of the second approximation is (Mp? |-C -|-U(p))g@—. DW (tq, 
The solution of this equation is 


g@) = -- DW (t)g™ (Mp? 4-C -|-U (py. (3) 
The desired constant error component of the MKA is 


(Atay) = (Au()) = (Ug), (4) 
where ())—=('(); ¢- =F fTo( .)d{ is a time averaging operator and T is 
the time during which averaging occurs. Substituting (3) into (4), we 
find 

(AU ytax) = (— UUW (t)q" (C + Uy)~'). (5) 
With regard to (2) 


(Atinax) == (W's (1) Sy (7) Wie (OD) + (We 8) S2 (7) Wi OM) (6) 


where 5 —1),.)),mil (Mp? |-C-|-U(p))-!(C-|- Up)“; Sg =U Dall (Mp? 4- 

iC -| U(pyy HC 4 Uy); DW(O = DW) + DWe (0. 
Expression (5) permits numerical estimation of the constant error 
component at known parameters of the MKA. It is obvious that the 


conditions of the invariance of the MKA to spatial vibrations have the 
form 


S, = 0, S, = U. (7) 
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To make the arguments more specific, let us introduce the generalized 
calculated model of a MKA with flexible suspension. Analysis of known 
layouts of compensation pendulous accelerometers shows that they can all 
be represented by a single block diagram (Figure 1). 
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KEY: 
1. Sensitive element 3. Amplifier-converting module 
2. Displacement sensor 4. Actuating member 
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Figure 2 


However, the control circuit may differ considerably in different 
accelerometers. Thus, a linear displacement converter (PLP) or angular 
displacement converter (PUP) can be used as a displacement sensor (IP), 
while a force converter or torque converter can be used as the actuating 
member (IE). The layouts of the flexible suspension (UP) may also be 
different. Let us consider two layouts when solving the problem: a 


ho 








flexible suspension (Figure 2, a and b) whose flexible elements 1 and 3 
connect the base 4 containing the sensitive element (ChE) of the 
accelerometer 2 by different methods. 


Limiting ourselves to consideration of a MKA with two degrees of freedom 
(X is linear displacement along axis Ox, @ is the angular displacement 
with respect to axis Oy), let us write the equations of motion 
corresponding to (1) in a coordinate system whose center coincides with 
the center of stiffness of the flexible suspension: 


(mp? -|- Cn) x -| mL p20 mW (key — Pi(d, 
mL p*x +4- (1, -}- ml.2) p? | Cy) 0 == — mLW (1) — Rucu (82y, -- 
— 8zsy2) —- M(t). 


(8) 


Here 1! — Rint? (p) (fh, (x 4-929) 4- 9229), where C, and Cy are the linear and 


angular stiffness of the flexible suspension, fyjny, Raniozs, Rao fy 472 


are the amplification factors of the linear displacement converter, 
angular displacement converter, force converter (PS) and torque 
converter (PM), respectively, Io + mL? is the moment of inertia of the 
movable part of the MKA with respect to axis Oy, 1 and b, are the 


distances to the axis of sensitivity of the linear displacement 
converter and to the effective axis of the force inverter, respectively, 
$(p) is the transfer function of the amplification inverting module, 
F(t) and M(t) are the force and moment caused by the influence of the 


spatial vibration of the MKA base on the stiffness of the flexible 
suspension, ;, Bo, yi, Yo, fi, f2 are constant coefficients that take into 
account the values >of 0 or 1 as a function of the composition of the 
control circuit (fi, fo, y1, yo, fi, f2 are equal to 1 if the control circuit 


contains the PLP, PUP, PS, PM, and UP, respectively, shown in Figure 2, 
a, and the flexible suspension shown in Figure 2, b. Otherwise the 
values of the coefficients are equal to 0). 


Let us rewrite equations (8) in the form 


A(p)x 4- B(p)0 = — mW, (1) — Fu, 


(9) 
C (p) x 4- D(p) 0 = — mLW, (t) — Mu), 


where 
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A(p) = mp? -}- Cy 4. RismikncD (p)pyy,; 
B(p) = mLp? -- Riki) (p) V1 (Byoz 4- B,oz5); 
C (p) = mLp® + Rnanknc (p) B, (v15z + 2525); an 


)(p) = (Ip 4- mL’) p? +. Cy + knanknc (p) (Boz + B,023) (y,5z + 


+- Y25zy). 


According to notations (10), the expressions for S; and S» in formula 
(6) assume the form 


= fombskunn®, [A,A (p)]-! x [a,a,, (ByozsC, — Bi%s) + OgQog X 
x (P44, —By9zsAp)], 
Sp == MP Rng [AA (p)]! (% (B1412% — Bf 11% 4 P0232 (f1A1C5 — (11) 
Avra) 1 9 (By%q (L A F122) —By%3@19— By5z2 (Ap (f1429-+ L) — Co%2))I> 


where % = LB(p)—D(p);_ % = C(p)—LA(p); a3 = Dyo— Coz; = 
= 1, —07A,; by = 2mLd-!. ' . _™ 


Expressions (6) with regard to (11) link the constant error component of 
the output signal of the generalized model of the MKA, caused by the 
action of spatial vibravion on the MKA, to the values of the vibration 
amplitude and parameters of the MKA. 


Let us consider as an example two design layouts of a MKA with flexible 
suspension. 


First layout. An accelerometer contains a sensitive element on a 
flexible suspension (Figure 2, a), linear displacement converter, 
amplification-converter module (UPB), and force inverter. Taking into 


account that /;=f:=yi=1, fo=f.—y.—0, for this accelerometer, we 


write 


S,=0, 


Se = mPbnca My [Agd (P91 [(Iqp? -+ Cy) (Ca29 + Cyc.) + al 


-+ Cyl (C yo, (L -}. 5.) -}- C\ Qo) -~-( (p) Runnki (8; a 1) Crcz (L -{ 


“+ A. — A497) 4- Cy (242 ~— 4 ,02)]. 


Second layout. The flexible suspension corresponds to Figure 2, b, 
while the control circuit consists of an angular displacement converter, 


en 











amplification-converter module, and torque inverter. For this 


accelerometer, h=—i=y,—0, fo=fPo=yo=1. Then 


S, = mbskiinnOzsC iL Wao, [A,A (p)]|7"', (13) 
S, = Mskinnz2P,C n [A,A (p)|-! [(/,p" -}- Cy ) ie -{- Cul -+. 
+ Rnanktic® (p) oz36220,5]. 


Expressions (13) show that a MKA, designed by the second layout, can not 
be made invariant to spatial vibrations due to variation of the 
parameters. This possibility is essentially present for an 
accelerometer designed by the first layout. Indeed, expression S; = 0 
means that the accelerometer will be invariant to lateral vibration at 
any values of its parameters. The condition of invariance of the MK4 to 
vertical vibrations (So = 0) reduces to imposition of additional 
requirements on the amplification factor of the amplification-convert ing 
module: 


WD (p) = ((/pp? + Cy) (Cnozay, | Cy-a,;) + Cul (Ca oz(L + Ayo) -+- 
--Cyaj,)] [RiicRiain (62— L) (Cy (L -| (ly9— Aya't7) -—— Cy (2449 7--044))]. 


Calculations made for a constant flexible element with rectangular 
profile permit one to write with sufficient accuracy 


D (p) & Cal [Ricknan (62 — L))—! + Ig(Cnozay, 4- Cyay,) p? X 
X [Ricknan (5z — L) CyLoz}. 


The physical impracticality of the transfer function, described by 
expression (14), limits the practical application of the result. 
However, if the system is subjected to low-frequency perturbations, a 


static regulator with transfer factor ()=(CyL|hnokinn (52—L) r can be 


used in the control circuit to achieve invariance to spatial vibrations. 
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[Text] One of the main difficulties in study of the dynamics of 
thermomechanical systems of liquid bodies and solids is related to the 
infinity of the number of degrees of freedom of the liquid and to the 
nonlinearity of the equations of hydromechanics. It can be overcome 
with approximation of an infinite thermomechanical system by an 
idealized model with finite number of degrees of freedom. The most 
preferable way of designing a finite-dimensional model includes 
replacement of the differential equations in partial derivatives by a 
system of ordinary differential equations of the Galerkin type. This 
direction is based on expansion of the solution of the initial equations 
to a finite series with respect to some, rather representative system of 
coordinate functions, which corresponds to design of an 
infinite-dimensional phase space of solutions onto a finite-dimensional 


space. 


For the temperature field 


Tir.j= Virol), rEg (1) 


em 


j=l 


n 


where fj; are functions of radius vector r, determined by their own 
projections in a suspension-bound coordinate system. 


The temperature field T(r, t) should satisfy the requirements of 
continuity at the interface of a solid and liquid body and the principle 
of minimum energy dispersion. Therefore, one can select the forms of fj 
as solution of a linearized equation of thermal conductivity 
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Vii(r)=90, rEQ (2) 


at inhomogeneous boundary conditions /;(r) = /;(r), i= i,m; r€aQ and of 


the linearized equation of thermal conductivity with respect to the 
natural forms of temperature distribution 


Vit) = df(t), rEeQ (3) 





at homogeneous boundary conditions /i(f)--0. i -m 4. j. rego Jil, om, 


where V? is a Laplace operator, 4j is the eigen-value, and fj are the 
forms of temperature distribution in region f}. 


As is known [1], the Laplace operator is negatively determined, i.e., 
its eigen-values are ).;<(). j=1!,~. This corresponds physically to 


averaging of the temperature through the region over time. We note that 
the longest are forms fj}, corresponding to the least eigen-values 4; in 
absolute value. This is explained by the fact that the eigen-values 4; 
determine the attenuation decrement of the corresponding form fj due to 
the diffuse nature of the process of thermal conductivity. 


The first group of expansion coefficients (1) wi, j = 1, m, 
corresponding to the forms of temperature distribution with 
inhomogeneous boundary conditions and given in the form of a time 
function, is the input effects of the thermomechanical model. The 





second group of coefficients #), j = m + 1, n, corresponding to the 
natural forms of temperature distribution with homogeneous boundary 
conditions, determines the variable states of the model. 


The region of liquid flow 2 is generally complex in shape, and, 
therefore, the linear boundary-value problems (2) and (3) can be solved 
only by using numerical methods. The finite element method, the main 
idea of which includes approximation of the continuous value by a 
discrete model, is the most convenient. This model is constructed on a 
set of piecewise continuous functions, determined on a finite number of 


E 
elements [2] o= LJa. where E is the total number of finite elements 
¢=21 


The region 11 occupied by the liquid can be divided into eight-nodal 
finite elements for a broad range of devices with floated spherical 
suspension of the sensitive element. 
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Solution of problems (2) and (3) on each element is represented in the 
form of the expansion f-—|N\}{j}y, where [VJ=[\i\>....\Vsf is a matrix 
of the linear functions of the shape of the element and {/}°= {/,/o/2... fe} 


is the column vector of the values of {f}* in the approximation nodes. 


Functions of form Ny can be written as follows: Ng~=0,125(-+ Een) (14-nnp) + 


--tts), 1,8 in the natural coordinate system ( -1<:, 1, t<1) 


Using the Galerkin method, let us find the approximate solution of 
differential equations (2) and (3). To do this, let us satisfy the 
condition: the difference between the approximate and accurate 
solutions should be orthogonal to the functions used in approximation. 
Minimizing the discrepancy of the approximate solution by the basic 
function in region {}., let us write for expressions (2) and (3), 
respectively 


' eT e ° 
\ IV) Wine dQ, {f}i= 9, i= 1,m; 
Q 


e 


(4) 


FLV y yt Ly dQ, 18 = Ay J INIT LY dQ, (8 (5) 
2, i» 





i=m-+1,n; j=Ii,n—m. 


Using integration by parts, let us reduce the order of the derivative in 
(4) and (5) and let us find a series of systems of linear algebraic 
equations with respect to the unknown values of the forms of temperature 


{f};, i = 1, m in the approximation nodes 
[A] {}; =0, i=Lm, (6) 


and also the generalized problem for eigen-values 


[A} (f}i =A (BI, i=m+.a, j= la—m. (7) 
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E E 
Here [Aj = [J [A]‘, [8] = LJ [B}’ are quadratic matrices of problems of 
e-=:] 


e--1 


linear algebra (6) and (7) whose dimension is equal to the number of 
approximation nodes N, where 





pAyt = ( (QUIT INIE LONI UNI 
~ dl dp op p> dv 00 
Re 
|_dINJ"" any) J wi” 2 ao 
+ p?sin? 0 = dy Oy } dQ, a, on ” 
e 


[By = \INj* [N]* dQ, are the quadratic matrices of the finite element of 
oF 

dimension 8 and p, §, and y are the spherical coordinates of the 

approximation nodes. 


The resulting matrices [A] and [B] are positive-determinate, sparsely 
filled and symmetrical with respect to the main diagonal. The first 
property permits one to use iterative methods of solution. 


The upper relaxation method, which permits one to operate only with 
nonzero coefficients of a densely packed [3] matrix [A], can be used to 
solve the system of linear algebraic equations (6). 


The value of the j-th variable on the k + 1-th iteration is determined 
by the formula 


j-! N 
_ k-4- oft 
he = fj -|- qajj » ajafs ‘_ p Qjq/s )s 
‘s=1 s==] 
where q is a relaxation parameter. 
The criterion of the end of the iteration process is max | /; a fil< 
} | 


<?, k<l,, j=1,N, where € is the accuracy of calculating the form and 


I; is the maximum number of iterations in solution of the system of 
linear algebraic equations. 
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Let us solve problem (7) by the power method [4, 5]. Let us assume that 
pj = 4-'j, and then 


LB) {7}5 = py (Al {f;. (8) 


Thus, determination of the minimal eigen-value of problem (7) is 
identical to determination of the maximum eigen-value for problem (8). 


Let us construct two sequences of vectors {X}k and {Y}k such that 


(Vn = IBN}, [AN{A}, =}, k= Too. (9) 


The sequence {X}, reduces [4] to the eigen-vector {f};, corresponding to 
the maximum value #; of problem (8). Let us find the approximations to 
the eigen-value from the Rayleigh relation 


= (Vege (AWN. (No k= Tree. 


Let us formulate the criterion of the end of the iteration process in 
(ait! — nt) (uty! <6, k-1,/7,, where 6 is a previously given small number 


and I» is a number that restricts the number of iterations. 


The eigen-vectors, corresponding to the higher eigen-values of jj, i = 





= 2, n - m, are determined by the above algorithm with regard to their 





orthogonality to each of the previous eigen-vectors {f})j, j = 1, i - 1. 
Let us assume for this at each k-th iteration of process (9) for the 
i-th eigen-vector 


i -| 
(X}n = (XJn-- May (Di 
j=! 
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where the unknown constants aj are calculated by the formula 


ee 


0 (XS, BUA, BUA) | j= THT. 


The first three eigen-forms in the cylindrical region (Figure a, b, and 
c, respectively) were calculated to check the constructed algorithm for 
calculating the forms of temperature distribution with homogeneous 
boundary conditions. The difference of the results of calculations from 
the analytical estimates from [1] does not exceed 5 percent. 


The developed algorithm for calculating the forms of temperature 
distribution permits one to construct effective mathematical models of 
the dynamics of semi-aggregate thermomechanical systems. 


BIBLIOGRAPHY 


1. Marchuk, G. I., "Metody vychislitelnoy matematiki"” [Methods of 
Computer Mathematics], Moscow, 1980. 


2. Segerlund, L., "Primeneniye metoda konechnykh elementov" 
[Application of the Finite Element Method], Moscow, 1979. 


3. Kureyev, A. S., "Solution of the Stationary Stokes Problem ina 
Quasi-Spherical Thin Layer Using the Finite Element Method," Kiev, 
1984., deposited at UkrNIINTI on 26 Nov 84, No 1940 Uk-84 Dep. 


4. Fadeyev, D. K. and V. N. Fadeyeva, "Vychislitelnyye metody lineynoy 
algebry" [Computer Methods of Linear Algebra], Moscow, 1963. 


5. Radysh, Yu. V., A. S. Kireyev, and O. P. Marosin, "Determination of 
the Forms of Liquid Flow in a Floated Suspension," MEKHANIKA 
GIROSKOPICHESKIKH SISTEM, No 5, 1986. 


a1 








UDC 531.313:532.58 
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candidate of physicomathematical sciences, Kiev Polytechnical Institute] 


[Text] A floated suspension (Figure 1) is a supported solid (float), 
immersed in a liquid and included together with it in the spherical 
cavity of a bearing solid (body). The region occupied by the liquid 
between the surface of the float and the surface of the body cavity Is 
is multiply connected and branched. Let us link to the float a closed 
surface [';, which separates the considered region into two subsystems: 
the float and liquid in region fl; with boundary [; and the liquid in 
single-connected region {ly with boundary [; UT... Surface [; consists 
of connectors in the liquid and in part of the surface of the solid. 


The purpose of the paper is to construct a mathematical model of the 
rotational motion of the float about a fixed point upon given rotation 
of the body and with a partially given field of the relative flow rate 
of the liquid on the connections, belonging to surface [';. 


Let us link the frame of reference and the coordinate system to the 
float. In this case, the rotation of the float will be transport 
[translator’s note: several words illegible] of body I» and the flow 
region of the liquid will not vary over time. 


The complexity of studying the dynamics of mechanical systems that 
contain solids, interacting with a viscous liquid, is related to the 
infinity of the number of degrees of freedom of the liquid in 
combination with the nonlinearity of the equations of hydromechanics. 
The flow nature assumes the motion of the liquid with considerable 
Reynolds numbers, which eliminates the possibility of using linearized 
equations. One of the methods of solving the problem is to approximate 
the real infinite-dimensional system by a model with restricted number 
of degrees of freedom, which reflects the significant aspects of the 
dynamics of this system [1, 2]. 
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Figure 1 


Let us assume the velocities of the mechanical system as linear 
combinations of the limited number of forms {,(r), q(t) (i= 1,n) of the 


possible velocity distributions of the float and liquid in region ly as 
the basis of kinematic concepts. These forms can be established from 
solut on of auxiliary boundary-value problems [3]. The speed of the 


liquid in subsystem 1 is y (f.4)-- fins. ef) (FEQ, UU) and that in 
subsystem 2 is V,(r,/)—¢; (rym (4) (rE, 1 PY UNS). mere 
e'.o' (i= 1,7) are the quasi-speeds of subsystems 1 and 2, respectively. 
Let e',a' (i — 1,3) be the components of the angular rotational velocity 
vector of the float, wi (i = 4, 6) be the components of the vector of 
the angular velocity of the body, let ¢'.m' (j—7,9) be the components of 


the vector of the flow rate of the liquid through the connectors of 
surface [;, which determine the transport of the moment of momentum from 


region fl; to region flo, and let ¢i, m' (i=—10,") be the quasi-speeds 


that correspond to forms of motion of the liquid, generated by different 
boundary conditions in a fixed coordinate system, bound to the float. 


The quasi-speeds yi (j—1!.3, m-+1.7). which are variables of state of the 
model, must be determined, while quasi-speeds vj (i = 4, m), given in 


the form of time functions, are input actions. 


Systems of ordinary differential equations 


l 
by (1) 


1 .., | 
J ’ foR o 
Lijé -|- | ine! 
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go! -| a, po! -|- [,j,enleat = P, (i, I k = T,n) (2) 


describe the motion of the float and liquid, respectively [4, 5]. The 
components of the positive definite numerical tensors of second rank of 
kinetic energy in regions ; and {2 are determined with regard to the 
rotating frame of reference by the expressions 


l —— a 
By = (ffdpQ (i= 1,3, j= 4,0), 
0, 


(3) 


gu = (ondeQ (i,j = Tn), 


Qs 


Here p is density. 


The components of the positive definite numerical tensor of second rank 
of the dissipative function of the liquid in region Sp with regard to a 
rotating frame of reference have the form 


sf] 


a = 2 A“ ° ont F e eee 
Vu ) Wii V'OdpvQ (i,j =Ta), (4) 


Here v is the kinematic viscosity of the liquid. 
The components of the tensor of third rank of dynamic connectedness 


between different forms of motion of the float and liquid are determined 
with regard to a rotating frame of reference by the expressions 


I 
Psy, = 0,5 { (fVF, + fivh) fidpQ (i= 1, 3; jik=T1, 3, 4, n); (5) 
0, 
| 1 aa nes pie 
Ving = Vin = (fyytfidpQ (i= 1,3; j= 40; k == 7,3). (6) 
0, 


Coefficients 1%, (i, j,k = j,n) are determined by expressions (5) and (6) 
after substitution of indices 1 by 2 and of forms f by ¢o. 
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| 2 
In equations (1) anc (2), Pj and Pj; are generalized forces corresponding 
to €! and w!, respectively, and applied to the liquid. 


Let us exclude surface (contact) hydrodynamic forces, adding equations 
(1) and (2) with indices i = 1, 3 and subtracting equations (2) with 
indices i = 4, 6. Disregarding the variation of the moment of momentum 
of the liquid in relative motion through the float ge! = 0 (i= 1,3, j= 4)n) 
and the principal moment of Coriolis forces of inertia acting on the 
rotating float from the direction of the liquid OI, elect = Oi = i, 3), 


passing through it, we find the equations of rotational motion of the 
float only with respect to quasi-speeds w 


") 


0 


, 3 . mes 0 
» Kea 1 Bij) w! + p2 (Vijn -F ijn) o'r | + » | gis! + ayo! + 
= Fa j=4 (7) 
eo 0 n 3 9 
’ P ’ a | . ; . J a 
-+- > I’; ;,@/@"| + 2 > > I';;,0/0* == 0 (i = 1, 3), 
k= juth=l 
0 2 0 2 2 
where ° : a “se an = 
Bij = Bij — Biya Ui Gis Vin = Vin — Pisin. 


The dynamic model (2) and (7) is closed with respect to variable states 
w! (i — 1,3, m+ 1,n) if the input actions w' (i = 4, m) are known given time 
functions. 

The vector of the relative angular velocity of the body in a rotating 


coordinate system, bound to the float, can be expressed by the vector of 
the absolute angular rotational velocity of the body in a fixed 


coordinate system «/, (j=4.6): 


(4-3 74 (8) 


° P i . 
@ <= ci) — @ (i = l, 3), 


where [c';] is a matrix of orthonomal transformation of a fixed 
coordinate system to a movable coordinate system. Differentiating 
expression (6), we find the components of the vector of the relative 
angular acceleration in a rotating coordinate system 


+3 — — of miptalt8 + el@it3—al (jf i bel. 
@) Cn; 7) CO, CM, 0) (1, }; k l, 3), (9) 


where e'«) are Levy-Chivit symbols. 








The dynamic equations of state in quasi-speeds (2) and (7) are 
supplemented by three kinematic equations of rotation of the body about 
a fixed point with respect to Euler angles when expressions (8) and (9) 
are used [6]. Matrix [c!j] is fully determined by the values of the 
Euler angles, which are independent parameters of orthogonal 
transformation. 


The structure of system of equations (2) and (7) was studied by method 
[5] for n = 22 forms of liquid flow with regard to the structural 
symmetry of region {l2. Analysis showed that the tensor of inertia (3) 
is restored from 18 independent coefficients and dissipation tensor (4) 
is restored from 96 independent coefficients, calculated by the above 
expressions. The time of calculating functionals (3)-(5) on a computer 
is reduced by more than an order of magnitude. 


For numerical integration of system (2) and (7), it must be solved with 


respect to the arbitrary unknown quasi-speeds. Matrix [gij] must be 
inverted for this: 


( | 2 2 2 2 — __ 
Bij + Bij — 8143, ~~ Bi j+3 + Bia sj43 (i,j = I, 3), 


Rij — Bi 4, (6 = 1,3, j=m+ ln), 
Bi = | . 


——— -_——-- 


Biy-— Bi jaa (1 -m4-1,n, j = 1.3), 





By (j= m+ Vn). 


a 


Theoretical analysis showed that [gjj;] is a nondegenerate positive 
definite symmetrical matrix. It turned out in study of the real values 


of coefficients [gj;] that the largest coefficients are located on the 
main diagonal of the matrix. The difference of the diagonal 
coefficients by five-six orders upon direct inversion of this poorly 
determined matrix [g;;] results in violation of the properties of 
symmetry and of positive definiteness to such an extent that numerical 
integration becomes unstable, while its results are uncertain. To avoid 
this, the inversion is made in three steps. 


n~ 


1. Transtion from matrix [gi;] to matrix [g*jj;] is: 
gi, = Titular (10) 


the order of the diagonal elements of which is identical. In expression 


(10) 7); -6,/Vg: (i,7-- 1,3; m j-i,n), Where §;; is a Kronecker symbol. 
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KEY: 
l. s*! =. 3 


This transformation guarantees conservation of symmetry, positive 
definiteness, and also equality of the elements of the main diagonal to 
unity. 


2. Inversion of matrix [g*jj] by the numerical method of 
orthogonalization [7]. 


3. Determination of matrix [gij]~! by the formulas 


“A 


gy! = Tye; y.- (11) 


One can ascertain the validity of relation (11) by a direct check. One 
can derive from (11) the equality gjj =T7;)'gij7i;, and having multiplied 
both the left and right sides of which by T, we find expression (10). 


The results of numerical simulation of the angular rotational speed of 
the float, which had cubic symmetry, at constant angular velocity of the 
body 


o(f)= 1 s7l (§ =4) (12) 


and zero initial conditions, are presented in Figure 2. Upon rotation 
of the body, the liquid in region {ly is gradually set into rotation, 
drawing the float after it, the angular velocity of which emerges to 
steady value 1s"! (curve 1). A similar pattern is observed (curve 2) 
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upon assignment by law (12) with i = 7 of a constant flow rate of the 
liquid through the connectors of the float, which determines the 
transport of the moment of momentum through surface I. 


6),¢°' (1) 








Figure 3 


1. gs! 2. s 


If the vector of transport of the moment of momentum of the liquid 
through surface [; is not colinear to the vector of the initial angular 
velocity of the float, the results of calculations indicate the presence 
of a gyroscopic effect during the transition period (Figure 3). 


Calculations made upon approximation of region {lx by 1,140 eight-node 
spatial end elements (1,800 nodes), occupy 700 kbytes of the main memory 
of the computer and require 2.5 hr on the YeS-1060 computer. 


The model permits numerical study of the dynamics of a floated 
suspension on the computer and optimization of it by the given figures 
of merit. 
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{Article by A. A. Leonets, candidate of technical sciences, Kiev 
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[Text] One of the promising directions of improving the accuracy of 
information measuring systems is compensation of the errors of measuring 
transducers (IP) by special processing of their signals. If the 
parameters of measuring transducers are unstable, a necessary condition 
of the effectiveness of this method is estimation of the errors of the 
measuring transducers during functioning of them. The principal 
capability of using bending oscillations of the body of an object to 
estimate the drift of measuring transducers is shown in [1]. 


Expressions that permit one to estimate the multiplicative and additive 
errors of measuring angular velocity transducers during functioning of a 
redundant measuring system, mounted on a rigid base, are found in this 
article. 


Postulation of problem. Let us assume that two measuring systems, 
consisting of three one-degree accelerometers and three one-degree 
measuring angular velocity transducers, are mounted on a rigid movable 
base. This type of system is described, for example in [2]. Let us 
denote the parameters of the systems by subscripts 1 and 2. The axes of 
sensitivity of measuring angular velocity transducers and accelerometers 
of the first and second systems form orthogonal coordinate systems 0;xyz 
and Osoxy, the unit vectors of which are mutually colinear. The 
position of point 0» in coordinate system 0;xyz is determined by radius 


vector R = (Rx, 0, 0)", Expressions must be found that permit one to 
estimate the additive and multiplicative errors of a measuring angular 
velocity transducer by the output signals of the measuring transducer. 
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To do this, let us present the output signals of the accelerometers 
a*¥,(k) and ato(k) and of the measuring angular velocity transducers 
#*,(k) and #*o(k) at k-th moment of time in the form 


ai (k) = (i: 4- AKg,)a(k) -4- Aa, + R. (k); (1) 


a(R) = (E + AKo,) (a(k)+ o (2) x (w(k) x R) + © x R)4- Aa (k)+ (2) 
-}- Eq, (R); 


oO; (R) = (E -|- AKo,)o (k) -} Ao, -f- Eu, (k); 


07 (k) = (E +- MKo,) (8) + Soy + Eo, (A), sl 


where E is a unit matrix; 


(AK;;), 0 0 
A jy => ) (AK 5 ;)y 0 » i=aea, j{j=1,2 
0 V (AK jj), 








are matrices of the errors of scale coefficients of the measuring 
transducer, = (dx, dy, Az)"; m= (Mx, @,, @2)™ are vectors of linear 


acceleration and of the angular speed of the base, ¢ is the angular 
acceleration of the speed of the base, and Aa; = const and Aw; = const, 
i = 1, 2 are vectors of displacements of the zeros of accelerometers and 
drift of the measuring angular velocity transducer. 

Disregarding second-order values, one can find the angular speed of the 
base from expression (3) in the form w(k)=(E—AKa,) (m*:(2) —Awi—Eqr: (A) ). 


Subtracting equation (1) from (2) and substituting the value e«(k), after 
transformations, we find 


Y (k) = H (k) X (k) + B(R)E(R), (4) 


where 
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Y (k) = a, (k) — a) (k) — @; X (@; X R)—e; X R; 
i= |H,: Hak, +t Hi Hax,,' Hax,,' E\; B= \ Hy H,i E|; 


























0 2wi Rx 20° R, 
Hy=|—Ro, —R,o; 0 |; 
— Ro; 0 — Rw, 
0 207° Rx 20°?R, 
Hoax, =| —00,R, —owiR, 0 ; 
— woo R, 0 —o7o'R, 
0 0 0 a) U 
H, =| 0 0 — e, Rs ° Haka, =| 0 a, 0 |; 
O eR, 0 0 0 ay 
a — (oF + of) R, 0 0 
Hax = 0 a, + @,0) Ry 0 
0 0 a+ 0, ©) Re 








X= Aors, ‘AK! AKa, | AKa, Aa, -- Aa, ("; 


E=|Eo, Ew! Ne, — ‘a, I’. 


Let us supplement equation (4) with a mathematical model of the vector 
of the desired parameters X, ordinarily determined during experimental 
studies of measuring transducers: 


X (k) =: D(k| kR— 1) X (k— 1) + BA), (5) 


where $(k|k - 1) is the transient matrix of dimension 15 « 15 and #(k) 
is white noise of intensity Q. 

Expressions (4) and (5) permit one to estimate the values of vector X, 
for example, using a Kalman filter: X(k|k)=X(k|k—1)+ K(A)LY (A)— 

- HX(k\|k—!)). Here X(k]k) is an estimate of vector X(k) at k-th 


moment of time, X(k|k - 1) is the predicted value of X(k) from the 
results of measurements on time interval [1, k - 1], and K(k) is the 
amplification factor of the filter. 
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When realizing the Kalman filter to guarantee the convergence of 
estimates to the value of vector X(k), it is necessary that the rank of 
the matrix of observability M[0, k] be equal to the dimensionality of 
the estimated vector (m = 15 in the given case): 


AM [a UJ == ['(k| 0) H' (0), "(Rk \ 1) H(A), ..., (kl k— Ll) A" (k — 1), 
H' (ky (0) D(R| 0), HC) @(k] 1), ..., D(R| k— 1) (k— 1), H (A). 


Analysis of expressions (5) and (6) shows that the components of vector 
X are observable through at least five measurement cycles at «)(i|i—-.1)= 


= const, w(i)#o(j) 40; a) #a(j) #0; 7, j—1,5; 14]. 


Taking into account that the additive and low-frequency multiplicative 
error components of the measuring transducer during a measurement cycle 
are considerably less than the high-frequency components, observation 
data Y(k) can first be accumulated to increase the ratio of useful 
signal/noise and to reduce the load of the computer upon realization of 
the Kalman filter. The elements of H(k), contained in the expression 
for determination of vector Y(k), are ordinarily periodic in nature and 
the positive values are compensated by negative values upon summation of 
them at different moments of time. Therefore, the signs of the 
measurement results should be taken into account to increase the 
signal/noise ratio upon summation. 


The process of accumulation must be divided into four cyclically 
repeatable steps to guarantee uniform convergence of the estimates of 
all the components of vector X to the true values. The group of 
elements of vector X is set into agreement to each step and is added 
with regard to the signs of the elements of matrix H(k) so that there is 
rapid convergence of estimates of the corresponding group of components 
of vector X. Assuming that X(k) = X(k - 1), the vector of observations 
in this case is determined from the relation 


(j4-4)n (i+t)n 
Vijy= VWWOYM=( Vv W , (i) Mi) X + 
i= fn | (—jntt (7) 
(i -iyn 
4+ NV W,DBWMEM, 


i-= ju -l 


where n is the number of observations accumulated on one step, and 
W, (i) = diag (signa; , signe” , signe, ); W, (i) = diag (1, Sigh( 1 sign(w? @ )); W, (i) = 


= diag (sign a, Signa’ ,signa’); W,(i) = E. 
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Expressions (7) and (6) permit one to estimate and compensate for the 
errors of the scale coefficients and drift of the measuring angular 
velocity transducers and accelerometers and thus to increase the 
accuracy of the measuring system. 
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[Text] The most effective methods of evaluating the state of complex 
engineering systems, for example gyromotors, are engineering diagnostic 
methods, based on decision-making theory. The purpose of the article is 
to solve the problem of classification of rotary systems (RS) with 
symmetrical ball-bearing (ShP) suspension. 


The state of a rotary system is determined mainly by the quality of the 
ball-bearing seat, i.e., by defects in manufacture and assembly of the 


~ 


ball bearing, which are described by random vector x = (xi, «.+, Xn). 


The ball bearing of the seat can be diagnosed by parameters xj, i=1, n 
only upon disassembly of the rotary system, which is not always possible 
in practice. The basic direction of modern diagnostics of rotary 

systems is nondestructive checking, which uses the determined functional 


relationships between the parameters of state xj, i = 1, n and the 
diagnostic observable signal z. The amplitude of the harmonic in the 
stiffness and vibration spectra of the rotary system will be this 
diagnostic signal z. 


For the considered rotary systems, the dependence of z on defects of the 
ball bearing has the form [1] 


2 = [(X p49)? -f (vyty)? -f Qryrgngr, cos x,]!°, (1) 


where x; and x» are the amplitudes of the harmonics of defects on the 
races of the first ball bearing, x3 and x4 are the amplitudes of the 
harmonics of defects on the races of the second ball bearing, and xs is 
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the phase shift between the harmonics of defects on the races of the 
first and second ball bearings. 


Variables xj, i = 1, 4 are random values, distributed by Rayleigh law 
with probability " dannity distribution [2] 





fa (u) = = ex (—Ser) u> 0: (2) 


and the random value of x5 has uniform density 


fu(u) =~, ue (0,2), (3) 


where u assigns the phase space of the random values of xj, i = 1, 5. 
All the random values of xj, i = 1, 5 are independent and accordingly 
the joint density f+ = (Uy, «+4 Us) for random vector x = (x;, ..+, x§) is 
determined by the formula [3] 


| Uy, 2+) Us) = [fa (4) 


Since z is a function of xj, i= 1, 5, z is also a random value. The 
density f(z) of random value z can be determined on the basis of formula 
(1) by integration of joint density (4). Since the values of density 


(2) contain unknown parameters Oy. , iz, 4, the joint density (4) will 


also contain them and this means that all ¢ i = 1, 4 will go into 


Xj ’ 
density f(z) in combination. The en ee ee become very cumbersome. 
Therefore, it is desirable to reduce the number of unknown parameters, 
which can be done on the basis of the engineering requirements on the 
ball-bearing seats. Using these requirements, one can show how the 
dispersions of the defects on the inner and outer races are related to 
each other for each ball bearing, i.e., one can assign the number 7;, 


i = 1, 4, which determine the ratios between 0,» i=-1, 4. Asa 
ms 


result, all ¢, are easily reduced to one unknown parameter @: 
si 
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—— 


Ox, =Vi9, i=1,4. The density f(z, 6) will then contain only one 


parameter @. 
To find expression f(z, @) in parametric form, let us find the 
distribution function z, using relations (1)-(4). By definition 


F, (8) = P {2 <s} = P£l(x1Xq)* -+ (x54)? + 2x1x5%5x, COS x,]'/2 <5}, 


where P is a symbol of probability and s is the phase space of random 
value z. Let us introduce in formula (1) the substitution of variables 


Bp = XyXq, Eg = XgXqp Es = COS Xg. al 


Since the random values xi, i = 1, 5 are independent, then §j, i= 1, 3 
are also independent. After replacement of (5), formula (1) assumes the 
form: 


2 = (t2 +82 + ikea)” (6) 


The densities of random values §,;, i = 1, 3 are as follows: 








tl ro l — tL we -}- ye dt ’ 
re) = ae fap en ao (ar te) 7 
0 
v, 0) = —2- | Lew|— ar +. 2 y }| au. 
I) = yz" \ 4 20° err y3/ J (8) 
. 0 
hy, (w) =[tVI—w yy", we(-- 1,1). (9) 


The distribution functions F,(s) is found by integration of densities 
(7)-(9): 


I’,(s) = | | | fe, (00% (0. 0) g, (ce) dudodie:, 


Je 
1). 
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where D; is the range of integration found on the basis of formula (6): 


D,:(u? 4- 0? 4 2uve)' 7=<s, 


Turning to cylindrical coordinates and having differentiated F:(s), we 
find 








\ l 
23 ‘ du ’ jv] 
0) = en laa * 
(2 ) ' Jo | a— ud ) (0? -- Qu -- IP | 
it | | yu" 0 —| 
- (10) 
- Z . u 
xX Ky manor EP Ts i dv, 
Vive (V° j- Quy |- Wy Vay gl? (0? j- Que -f- It 





where Kj is a modified Bessel function of third kind [7]. Having used 
the known approximations for Kp, we find 


j (2, 0) = pO ‘2 exp | - os 2] , (11) 


' 


where the approximation error is |¢,— [0 ', ¢ = 1,03806, p= (4% °)/ 3yr). 


A series of independent measurements of signal z is performed to 


estimate parameter @ and sample {7,!;_, is found. Having sample f2,} 


and using the method of moments [3], due to which the computation is 
simple, we determine 


oo 


, 


mn, IV: nny | 2/ (z, O) dz, Mm, My. (12) 


‘ 
i | () 


7? 


where m-, is the first sampling moment and my is the first theoretical 


moment. 


The theoretical moment m, can be calculated in analytical form: = m,—5¢- 192. 


Having set m, equal to the sampling moment mz, it is easy to find the 


I 


estimate @ of parameter @: @ = 1.0075. 








aA 


One can turn after calculation of estimate @ to solution of the 
classification problem. To do this, let us use the Bayes criterion [4], 
since it is the most effective. Let ts compile the probability ratio 


te) ee, Ne me ST. PP ed, 


(Cig — Coe) Po (13) 


where f;(z) and fo(z) are arbitrary densities z of classes $); (suitable 
products) and {ly (unsuitable products), respectively, P; and P» are the 
given probabilities that classes {); and {l, will appear, and {(Cij}, 


i, j = 1, 2 are elements of the loss matrix when making the decision. 
In our case the probability that an unsuitable product will appear is 
very small: P» € P;. However, the cost of passing a defective product 
is considerably greater than the cost of a false alarm: C,2 } Co;. 
According to [5], one can select dn = 1 in this situation. 


The main difficulty of using the Baves criterion is in the laboriousness 
of finding complex densities f;(z) and fo(z). The procedure of 
numerically finding these densities analytically on the basis of 
functional relationship (1) and of a priori description of classes {| 
and {ly is proposed in this paper. The classes are described in 
practice, using the technical documentation for the ball bearing. Based 
on the requirements on the ball-bearing seats, one can indicate the 


tolerance a,, i = 1, 4 for the i-th defect of ball bearing x, 
(permissible misalignments of the races, dimensions of the balls and so 
on). It is thus easy to assign classes §); and {lb in the following 
manner: 1) the product belongs to class {); at 0 < xj < @, and 2) the 


product belongs to class {ly at xj > a; (at least for one i, i= 1, 4). 
We find the conditional density f;(z) with this assignment of classes, 
using formula (1) and taking into account the constraints on variables 


x,: O<¢ x; £ @, 1 = 1, 4 accordingly, from the condition of 
normalization for densities (2) 


a 


Bi ( f.,(u)du = | 


one must find the normalizing coefficients 


2 ~~] 
, —- -- (_ ies BK . 
(i, | Il — exp | 5 l 
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After approximation, the conditional density f;(z) assumes the form 








A 
, i=| . 
h@=G \ Vana V (z, w) dw, 
1 T] ¥ 0° I 
i=! 
where 
bd (z, w) = Py (w) zHe(w) — exp |-—— ly (u’) Ziel) [}ty (u')— jt, (wm) jt,(u’) zits) |; 


2) A) —_ 
}t, (w) = ve" (1 - 1)"* Ive —(w+ I)! +g, R= 1, 


Having calculated the densities f(z) from (11) and f;(z) from (14), the 
conditional density f(z) can be determined by the formula 


a nr 


i, (2) = pe (f (2) — Pf, (2)). 


After substitution into (13), densities f;(z) and fo(z) were 


approximated at @ = 1, aj = 2.4426, i=1, 3, aj = 4.885, and i = 2, 4 
by simpler distributions while retaining high accuracy (|«| < 8.1-10°5) 
P (2) = 9,2 exp (— 7,2), fo (2) == P92 EXP (— T22). (15) 


where p, = 0,44564, 1, = 2/3, z€(0,0,), b, = 12,3; for f (z) and p,=2,119.107?. 
1, = 1/8, 2€(0, 0°) for Fy (z). 


Having substituted (15) into equation (13) and having solved it with 
respect to z, we find the number zo which separates classes Q; and Qo. 
If there are probabilities that classes P; = 0.9, Py = 0.1 and dy = 1 
will appear, Zo = 5.604999 will be found. Thus, the product is assumed 
suitable at z < zo and is considered unsuitable at z > Zo. 


The quality of recognition is evaluated by the probabilities of errors 
of first and second kind [6]. If the object belongs to class 1; and is 
assumed to be an object of class $l), this error is called an error of 

first kind. On the contrary, if the object belongs to class {5 and is 
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erroneously related to class {};, an error of second kind is made. The 
most complete characteristic of the quality of recognition is the 
average risk R, which takes both types of errors into account. The 


probabilities of errors of first kind a’. of second kind q and 
average risk R are determined in the following manner: 


tJ 


Ql = (Alerde, QF = \igerde, 


26 0 


R= P,CyQt +- PiC4,Q:. 


Since we used the value Jp = 1, then Co; = Py, Cys=P,, R=r(Q? + Q), 
where r= P,P,, r= 0,09, Qy = 1,05738%, Qh = 4,83152 %, R = 0,5889 %. 


Errors of the first and second kind and the average risk were small, 
which indicates high quality of recognition and permits one to recommend 
the suggested approach for solving problems of evaluating the state of 
symmetrical rotary systems. 
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[Text] Close attention has traditionally been devoted to study of the 
dynamics of precision gyroscopes. An object of detailed study is the 
influence of the parameters of ball bearing (ShP) imperfections on the 
position of the resonant zones of the devices [1]. The parametric 
effect on the vibrations of the devices as a perturbing factor have not 
been considered previously. Study of the dynamics and calculation of 
the natural frequencies of a surveyor’s gyrocompass upon the parametric 
effects of an imperfect ball bearing are the purpose of this paper. 


Let us take as the perturbing factor the deviations of the coordinates 
of the centers of curvature of the outer races in the axial 
cross-sections of the ball-bearing seat along the rotor axis from their 
true values when deriving the equations of motion of a surveyor’s 
gyrocompass, thich is the sensitive element (ChE) on a torsion bar 
suspension. Disregarding the interaction of the sensitive element with 
a rotating earth, let us assume that the rotor in the device makes only 
axial displacements and that its position in the device is ideal, i.e., 
there is no unbalance. 


Let us introduce into consideration a fixed coordinate system 0; X,Y; Z; 
(Figure 1), bound to the center of mass of the sensitive element without 
a rotor in the initial unperturbed position, and movable coordinate 
system OnXoYnZo, rigidly bound to the center of mass of the sensitive 
element without a rotor. The position of point Oo in coordinate system 
~ ~ 
O;X:¥1:Z; will be given by vector r; = ri(x1, Y1, 21) at an arbitrary 
moment of time. Let us rigidly link coordinate system Cf£7¢ to the 
center of mass of the rotor. The position of the rotor at arbitrary 
moment of time will be given in coordinate system OnXnYoZo. Let us 
determine its axial displacements by coordinate z. Let us give the 
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angular displacements of the sensitive element without a rotor by a 
system of angles of finite rotation a, 8, 7. The sequence of 
revolutions is shown in Figure 2. Let us determine the rotation of the 
rotor along the axis of its own moment of momentum by angle 'p’ 



































Figure 2 


Let 1 be the length of the torsion bar, letl; be the distance from the 
point of attachment of the torsion bar to the sensitive element to the 
center of mass of the sensitive element without a rotor, let lo be the 
distance from the center of mass of the sensitive element without a 


rotor to the center of mass of the rotor, let By HH (tox. yy, O20} and 
ime (Fro tig Cre be vectors of the angular velocities and accelerations of 
the sensitive element without a rotor in projections onto the axes of 


coordinate system OoXoYoZ, and let = {en, wy. wr} be the vector of the 
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angular velocities of the rotor in projections onto the axes of 
coordinate system Cf£9(: 


Ke, = 0 ly, U , 














is the torque of the sensitive element without a rotor with respect to 


pole Oo, 


Ke--]} 0 Ty Olle 








is the torque of the rotor with respect to pole C, To Tyo? and I, are 


the moments of inertia of the sensitive element without a rotor with 
respect to axes O9X0, OoYo, and OoZo, respectively, Tes Ty and ‘ are 


the moments of inertia of the rotor with respect to axes Ces Cy 


and 0, Cy, Ces ty (%mVO}s 7 (0,0, 2}5 x9 = —- Lp. 


and Co 


Let us use the Euler-Ishlinskiy method to compile the equations of 
motion of the sensitive element. Considering the interaction between 
the bodies--a sensitive element with a rotor and one without a rotor--we 


find the following equations of motion in vector form: 


dake aad — — > iad > 4 
a + W, X Ko, = Lo, — Mo, — (Fo sa r) x Q. 


Here M is the mass of the rotor, m is the mass of the sensitive 


= ~ 


(1) 


element 


without a rotor, @ is a vector of forces acting on the rotor, P is the 


ae 


vector of forces acting on the sensitive element without a rotor, g is 


7h 








-_ aad ~ pond 


free-fali acceleration, and W. , W. and M. , M. are vectors of the 
Oo C Oo C 


absolute accelerations and moments of the centers of mass of the 
sensitive element without a rotor and with a rotor, respectively. 


— 


Let us write the components of vector w; (see Figure 2): 


w,, = Bsiny + acosB cos y, w,, =B cosy — asin ycosB, 
ys (2) 
@,,= y + asin. 


The total potential energy of the sensitive element includes that of the 
sensitive element in the earth’s gravity field fl; and the energy of 
elastic deformation of the imperfect ball-bearing seat [lp [2]. With the 
adopted assumptions we write 


I], = Mg (x* — xo) + mgx,, 

2 =. + 

5 y Knw/ ’ 

j=! (3) 
wy = {(p + Ps — Da + Pa)? + Lz + (— 1) pg — (— i)! (05 — 2)? — 


— (py + Pp, — pi}? . 


Il, = 


Here pj are the parameters of imperfection of the ball bearing after 

V. F. Zhuravlev, z° is the tension in the gyromotor, n is the number of 
balls, K is Hertz’s coefficient, x* is the coordinate of the center of 
mass of the rotor in coordinate system 0;X,;Y;Z,: 


x* = x, + x,cosB cosy + zsinB. (4) 


When writing the potential energy, let us take into account that the 
introduced coordinates are related in the following manner [1]: 


x, —1—1, = —IcosA cosy — 1, cosh cosy, 
y, = —lsinu — I, (cosasiny + sinasinB cosy), 
2, = lsinAcosp — I, (sina sin y — sinB cos a cos y), 


where A and p# assign the position of the torsion bar in 0,;X,;¥;Z;. 
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Hence, we find with regard to the smallness of the angular displacements 


¥ = [Yat L(y + OB) 20 + [ar +h (ay — B22 + 1872+ Ly72. 9) 


Thus, the expression for the total] potential energy with regard to 
(3)-(5) will assume the form 


T= (M + m)g {ly +A (y + oB)P(20) + fz, + 4 (ey — B20 + 


. a 
+ 18/2 + Lyy"/2} — Mgxo(B* +- y)/2 + Mgfz + = Kn {{p, + 
j=1 


; (6) 


+Po— Ps + Py)? + le + (— 1)", — (— 1)! (pg — 27}? — 


5 
—Pi1—Pr+pr}° ; 


Let us find the expressions for the moments and forces on the right 


sides of (1). Denoting Mo, = {M,,, M,,, Mz,}, Mc = {Mz, M,, Me}, Oe = {Q:, Qn Qr}, 
Qo, = {Cea Quer Qe)» Poy = {Pros Pur Pe} we write 


Le, = Q = Qn = Qe, = QW, =" 
Ly, =(M + m) g {lz + 1, (ay —B) 4 — BVI + Mg (xf — 2), 
Lz, = Mgxoy —(M + m)g {lyr + (vy + By A + lvl, 


M; = M,, cosy, + M,,siny,, My = My,, cosy, — M,,siny,, (7) 


M = M,,, 


Py, =—(M+mgig +h (y+ op), Pz, =—(M +m) x 
x g (2, + l, (ay — B)/1, Qe — Mgb— (Kz -+- K apPs) ~— K pPs- 


Here 


call 
“Ge = (Kz + Kapps) 2 + Kp, (8) 


where Ks K_, and K, are the coefficients of expansion of the elastic 


p’ p 


force acting in an imperfect bal] bearing into a Taylor series. 
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Based on equations (1), (2), (6), and (7), we find the following 
linearized system of differential equations: 


Ay {oByy,2,2}" = {f,f,00f>f,}.- (9) 


Here i; = Q; is = LK pPgl(1y, -}- !,);  P = K yp,/m; fi. oo (M-+-m) x K,p,/(mM) _ 
~ BK psy + 1s), H = Tey. ly = by = Ins 





ii 1p 0 0 0 0 
Hop p*-4- By 0 0 — B, B, 
() 0) p?-D, D, () 0 
\, 0) () Cl, pe+c, 0 0 
0) C, 0) () p-+-C, - bal] 
Intl yp N, Y 0 Ne pr-+Ns, 











where p is a differentiation operator; 


Ky = K,4- Naps 1, = Hig, 4-15); Hy = Hl, +13); By = 
= (M +- m) gl,/[LC,, -|- 7s); Bp = [Mgl, 4+ (M 4+ m) gl, (+ 4)// 
My. 1s); By = [Ag — 1K o/(1y, 4- 1s))3 Cy = (M 4- m) g/ml); 
Cy = (AL -}- m) giy(ml) 4- Mg/m; Ny = — Bal -- Cy 4- g; 
No=- —C,-} IB, Na = — Bale 4- (m+ M)K,/(mM). 


Setting the determinant of matrix A; equal to zero, we find the 
following characteristic equation: 


p? +- d,p* + d,p® + d,p* + d,p? 4- d, = 0, (10) 


where 


d, == 2C, -+- Dy -}- Ny +- By -t- Hyg; dy = 2C,D, + Ci -+ 2C,\Ny + 
|-fDoN, -+ (2C, + D, + Ng) (By + 4,4,) + K,N2/m — D,C,l,—B,C, + 
4- ByN, 4 1,11,H,By;  d, = CD, 4- 2C,D,N, + CiN, + (2C,D,4+C)-+ 
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|-2C,N, --D,C,l, + DN 5) (Bz + H,H2) — (Cy + Ny) D Cyl, + 
- (By -}- HH, 4- Cy 4+ Dy) KyN2/m + (Ny — 1g, H2) BK /m 4+ 
|- ByC,N, — B,C, (Cy +- Dy + Ns) + (2C, + Dy) (BgN, + 12,1485); 
d, — C,(N,C, — B,C, 4+- KyN2/m) (Dz — Dy) + (CD, 4+ 2C,D,N3 -+ 
| CIN, — CD, — NgD, Cyl, + (Cy + Dy) KoN,/m) (By + 14, 11,) +- 
| (Cy -4- Dy) By Kg (Ny — let H,)/m 4- (C, + D,) (By — By) CaN g + 
4}. Buy (2D_ + Cy — Dyly) (Ny + lal lls); ds = (N3CT + C,KoNo/m) X 
x (Dg —Dyl,) (By -1- TH.) 4 (CBN Ko/m — Cy, 1B Ko/m 4- 
| CiBaN , -|- Cilell 11_B3 — ByC\C,N'5) (Dz, — Dyl;) + BsC Ne X 
x (C.D, — Dilly). 


The roots of characteristic equation (10) give us values of the natural 
frequencies. At 3 = y; = 0 


, — {25/g,(1 - 2iks @)}' (11) 


2,3 = {(Bi ~~ 2g + Mei — 22)” — 4 (g3 — 2g,e,)]'7}2}"", 


where 


By = Ny -l By |-Cy-le MMe: fg = - gC B, -1- BLN 4+ (Om 4M 
x Ky(C, -| H,11,)(mAl) -|- C\By | CH, H, — BN, | Ky ,/m-— CAR; 
Ra == — K,(C,B,-— C\R, — C,H VAL -- (BK g/m — B.C,) g. 


Let us calculate the values of the natural frequencies of the sensitive 
element with the following input data: M = 0.6 kg, m = 0.9 kg, 1 = 

= 0.16 m, 1; = 0.15 m, lo - 0.16 m, I, = 12.74-10°4 kg-m?, I, : 

= 13.72-10°' kg-m?, I, = 7.22-10°4 kg-m°, I, = 58.8-10°4 kg-m°, Io = 
9.8-10°4 kg-m2, H = 0.4 kg-m2/s, K, = 0.43 kg/pm, ps = 43.5 ym, and 


[os 2 
Ky = 0.086 ke/am?. 


Formulas (11) yield wo! = 7.75 —; wl = 815g ~!; wl = 8295 ¢ | 


Comparing the value of us to the value of nutational frequency 
o* = /V(1,, + 15), -|- 1.) = 107 s-!, we see that vs is less than 20 


- 





percent less than wt, which is in agreement with known experimental 
results. 
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The roots of characteristic equation (10) yield refined values of the 
natural frequencies of the sensitive element: o!!=-6,5 g-!; wll =&8 g§ |: 


, 
wa! 24,4 s I wt! &3 s~- Ie wl! 8204 s° |. 


Thus, one can calculate the fundamental natural frequencies with 
sufficient degree of accuracy by formulas (11). Consideration of the 
cross stiffness of the ball bearing Kp in expansion (8) permitted us to 


. : ; I 
determine the accurate values of the maximum resonant frequencies w 3 


and ty. Numerical integration of system (9) by the fourth-order 


Runge-Kutta method at nutational frequency ws = 81.5 s~! with 
perturbation amplitude of parameter p5; equal to 3.6 ym showed the 
presence of a second parametric (harmonic) resonance of the studied 
system. 
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One of the characteristic features of the dynamics of pendulum devices 
is the presence of unstable parametric vibrations upon vertical 
vibration of the base, which can be eliminated by introducing the 
required damping into the system [1]. Let us consider these vibrations 
in a pendulum, vibration-insulated in the horizontal plane, studying the 
case when the vibration frequencies exceed the resonant frequencies of 
the system. 





A body of mass M (the body of the device), connected to the base by 
springs that perform the function of vibration insulator, moves in the 
horizontal plane (see figure). A pendulum of length 1 and mass m is 
mounted in the body. The base completes verticel vibrations at 
acceleration Wy = ©W cos {lt, where € is a small value. 


The equations of motion of the system are written in the following 


manner : 


(M + m)x ++ mip cos p + h,x + cx — mig sin @ = 0, 
(1) 


lp 4 hep 4 mi (Io + x cos — xpsinp) + ml(g + W,) sing = 0, 








where x is the displacement of the body of the device, ¢ is the angle of 


rotation of the pendulum, Ij is the moment of inertia of the pendulum, c 
is the stiffness of the elastic coupling, and h; and hy are damping 
coefficients. 


Limiting ourselves to second-order values inclusively, we transform 
system (1) to the form 


D+ Ep + hi(1 + ry) + hig 'x = 0, (2) 
x + Chex + kx + bp = 0, 


r 


where k, = | mgl . k, = \/- : are partial frequencies, I = Ip + 
m 








y+ 


+ ml*?; b= qn ~i ty=g'W,; and {, ( are relative damping coefficients. 





Let us use the averaging method to solve the system of equations. Let 
us find the solutions in the form 


—lwjt 


2 iw ;t 
g=VCe'+De ", 


om (3) 


r. tw jt —to jt 
x= YH (Ce +De °), 
j=l 
a bw . 
where C) and Dj are complex unknowns, j = V—1; 1 = 77 1 is a P 
2 — Vig , 


coefficient of forms, #)o is the frequency of natural vibrations of 
system (2) in the absence of damping, and «; is the vibration frequency. 


let us introduce frequency tuning cy), having written 


S a- 


w; = wi (1 — en). (4) 


Having substituted (3) into system (2), we find by the method of [2] 








C, =A (e"™ (Riv, — R,d,)), 


D, =—A "td (R,¥, — Ryd,)), 


c = Aq! (e~ ft (R.d, — R,v,)), 


Dy = —A~' (e'™ (Rad, — Ryv,)), 


wher‘t 


A= dy,— dv; Ry = —ajA, — p (0 4. e~'8) Ay — 8)B;, Ry = 
= — gyAy — f,B;; d; = a;2iw,;; U; = 21W;;; a; = Rien yay; 
bk? 


2 2 


qi = kenny, a= 14kig Wy; x)= fy = Chaptyior,; 


l | -_ 
p= > Rit yo; $s = ER, i@y; Tyo = go 'eWV,; A, ae Cel 4 Dye a 


B, = Ce’?! — Dei", 


Analyzing expressions (5), we find that non-trivial solutions are 
possible in the cases 


Rye !e) = —n,C,—pD,, (Rye) = — uD, — pC, (Rye!) = 


0 =— NC, (R ef) = U,D,; 
b ) 2W = ° 


(Remit) = —n,C,, (Rei) = —u,Dy, (Rye) = — 1g, — PD 
(Rye!) = — ugD, — pC, 


c ) Wa + @, = 2: 
(Rest) = —n,C,—pD,, (Rye!) = —uyD,— Cy (Re-"*) = 
= —nC,—pD,, (Rye!) = — uD, — pC,; 
d) Wg —@, = 2: 
(Rye!) = —nC,—pCy; (Riel!) = —u,D,—pDy, (Rye "#) = 


= — ngCy—pCy, (Rye!) = — uD, —pD,, | 


where fig = dg + Sq} Ua = dg — Sa; & = 1,2. 








The following relations occur for all four cases 


( Rye) = — (9, + fi)Cy (Ree) = —Qi — f:)Dy,  (Rge") = 
= — (9+ /2) Cy, (Rae!#) = — (qa — fs) Da. 


We assume that ws > «;. Cases b and c are of interest with respect to 
the considered postulation of the problem ({} > wo). 


If condition 2wv» = ft is fulfilled 


AC, = [— 0, (a, 4+- s,) +d, (9g, + f,)) Cy, 


AD, = [v, (a, — $,)—d,(q,—/,))D), (6) 


AC, = — dy (do + fy) Cy + 2; [(@, + Sg) Cy + pDgl, 


AD, = d, (42 — [,) Dz — 4 [(Ag — $2) Dy + Cy). 


We find the conditions of stability of the solutions of system of 
equations (6), which is divided into two subsystems. Analysis of the 
first subsystem, which includes the first and second equations, shows 
that its solutions are stable. The condition of stability has the form 


ke w?, 
on ER, -+ Oe © ~ Ck, > 0. 


7 5 ; 
Wo9 — kp 2 9 








Let us consider the second subsystem, having written it as follows: 


AC, -+- uC, +- dpe = 0, UC, +- AD, + Up Ds =), (7) 


where ic = ds (qo + fe) — U4 (42 + 82)5 Hyp = — UP) Ug = Ps Hyp = — ds (42 — /2) +- 


+ eq 


The following conditions should be fulfilled for stability of vibrations 


Wie t+ Upp > 9, Uy cllay — Ugg p > 9. (8) 
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The first inequality can be transformed to the form 


trie _Bhska_ 5 9 
eer yomemeet, I D) 2 (9) 
Wy — Rp ky — 10 


Since w»5 > Ky and v;o < kd, inequality (9) is always fulfilled. 


The second inequality can be written in the form 
\9 .2 ' : 9 2 
(a4, 0Ro}t2 — ER,%,)°w2 + fey (#42)? (44%, — Hy)? —- Hip? > O. (10) 


Let the frequency tuning «7» be equal to zero. Inequality (10) is then 
simplified: 


(hy — Oy) 30 lo. a 
Wap | %6 —>—-— 5 -— + ER > > Rit. 


(5, — k5) ke (11) 


Thus, if inequality (11) is fulfilled, the vibrations will be stable. 
If damping in the vibration insulator is very small, inequality (11) is 
simplified considerably and the condition of stability is written 
thusly: 


(12) 
Accordingly, stability can also be guaranteed with small damping in the 
Vibration insulator by appropriate selection of damping in the pendulum. 


Let us now consider the case #» + #; = ft. System of equations (5) 
decomposes into two subsystems: 


AC, _ dC; -+ VD, = V, AD, +- A pDa -f- AC, om 1: (13) 
AD, +h, Dy oyl2 = 0, AC, + Aye Cy + AyD, = 0, (14) 
where Mig = ValQy + $1) — 42 (M1 4 AY hy = — Ua (Gy — 81) +49 (4 — hs Mog = Uc 
Nop = Uys Ay = Vals he = —Ays Ag = — 0) Mg = hs. 











The characteristic equations of system (13) and (14) have the form 


A?p? + A(Ae + Aap) P + AicMon — Ardy = 9; (15) 
A?p? + AA.) + Ave) PA Ap hae — Maks = 9. (16) 
Let us assume that €47; = €#9 = 0. The conditions of stability for 


equations (15) and (16) are then written in the following manner: 





» / 2 5 
2 k’ 10 | >0 
Rk) —__*—. } + Ck,| a ae a a ee 
Ek —_————_——— t- 7 > r 32 22 42 it) - Ry 
oi | 9 2 2 Rs ki — Diy 20 2 
ks — Wig M9 D) « 


hk. ; Wig ) __ tp { 
——-——— § ak, | 7 
40) 1920 (= ~~ 42 ER, }- Habe > “ane Re 0% 

20 


"4 a, Ck, w, — ke i “(ke iy) (Wy k’) 


The first condition is almost always fulfilled. Nonfulfillment of the 
second condition means the presence of unstable solutions. If damping 
in the vibration insulator is very small (¢ = 0), the second inequality 
assumes a rather simple form 


ky - (17) 


Since w1o < #99, inequality (17) is not fulfilled at smaller values of 
ryo than inequality (12), and accordingly the unstable vibrations at 
frequency {} = #;0 + #20 occur at smaller values of the amplitude of 
vibration acceleration of the base than at frequency {l = 2wgo. 


Analysis shows that a pendulum, vibration-insulated in the horizontal 
plane, can make unstable parametric vibrations upon vertical vibration 
of the base. When vibration insulators with small damping are used, the 
stability of vibrations can be guaranteed by selecting the appropriate 
damping coefficient in the pendulum. 
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[Text] The effect of uniform rotation of the base on evolution of the 
surface acoustic wave (PAV) and the possibility of developing a 
miniature angular velocity meter are considered in this article. 


Substrates, which have anisotropic elastic properties, are ordinarily 
used in converters based on PAV. But since the piezoelectric effect is 
only a smal] perturbation of the considered wave, which is primarily 
mechanical in nature [1], let us study the propagation of the wave along 
an isotropic substrate to simplify further computations, without 
disturbing the qualitative pattern of dynamics of the PAV. 


According to Newton’s third law, the equations of motion of the PAV are 
presented in the form 


Orn. Orn, On, 9 
) ~s75 — Cinnt 4 2 innStn Zp 7 p (Q,2,,u, — Qnty) = 9 (1) 
nol 


(i,n,k, l= 1, 3), 


where p is the density of the substrate material, C;pk] is the tensor of 
the elastic constants, €ijnk is Levy-Chi: it density, §) = {f), Wo, Mg} is 
the angular velocity vector, u = {u;, us, ug} is the elastic 
displacement vector, and x;XoX3 is a Cartesian coordinate system. 


The circular frequency of the PAV we is ordinarily much greater than 
angular velocity {l. Therefore, the third term of equation (1), which 
characterizes the Coriolis force of inertia, is on the order of 
smallness 4 = {)/w € 1, while the first term that describes centrifugal 
forces is on the order of 62. Let us further disregard consideration of 











the effect of the centrifugal term of equation (1) on the dynamics of 
the PAV because of its smallness. 


Let us represent the generating solution of system (1) by the expression 


u, = a, exp [jk (b,x, — vf)], (2) 


where 4, is amplitude, . is the phase velocity of wave propagation, k = 
= w/v is the wave number, and b; are direction cosines, corresponding to 
Cartesian axes xX}. 


Let the surface wave be propagated along axis x; and be attenuated along 
axis x3. Then, according to [1], b; = 1, be = 0, and bg = b at the 
selected orientation of propagation of the PAV. The value kb 
characterizes the depth of penetration of the surface wave into the 
substrate. 


After substituting solution (2) into system (1), we reduce it to the 
following form: 


(Vi — §,,¢v" — 2if = Fn Qn 0 = Q, 
(3) 


where Ty = Un = Pr ninst ig a Kristoffel tensor and 4;| is a Kronecker 


symbol. 


When a PAV is propagating in an isotropic substrate c,nk|, there are 
only two independent components c;,; = 4 + 2# and cq4 = #, where A and 4 
are elastic constants (Lame parameters) of the medium. The constraints 
imposed by symmetry on the tensor of elastic constants results in the 
following expression: 


( Cry 4- Cg” 0 (C1) —Cy,) 6? | 
Pi == | Q) C4, (1 4- 6°) Q 
(Cy, = Cqy)D? V Cag tH Cy,0” 


Let us consider the characteristic equation of system (3) 


(v2 4 wtb? — v2) (09 + v7? — vt) — (e? — en)? — Atv! = 0 (4) 





as an algebraic equation with respect to b at given value of v. In 
equation (4), v, — (¢,,/p)'”, v, =(c,,/p)"2 are the phase velocities of the 
longitudinal and transverse waves, respectively. 

It follows from the necessary condition of approach of the amplitude of 
the PAV to zero at x3 - » that b is the root of the characteristic 


equation (4), lying in the lower half-plane, solving which by the 
perturbation method [2] with accuracy up to ¢(&%3), we find 





27,2 
where g = 


Uv; — Uy 





Thus, the PAV is a linear combination of three partial waves: 


3 
up == VC, af” exp (ik x — 0d). (5) 


n—| 


where C, are the amplitudes of the partial waves and @‘"); are the 
polarization coefficients of the components of the partial wave, 
determined from system (3). 


Let us use the boundary conditions on the free surface of the substrate 
to find the phase velocity of the surface wave: 


au 


‘ k ee 
I au = Camnhl Ox) at Xa : (), (6) 


where T is the stress tensor. 


Having substituted solution (5) into boundary condition (6), we find 


dinGn i 0, (7 


C) 
Y 


where d., = ComniBe by”. 











_— 
i 


The condition of nontriviality of solution of expression ( 


det {din | = 9 (8) 


is considered as an implicit equation with respect to v. 


At 1 = 0, equation (8) is transformed to the form 


where vy is the phase propagation velocity of the PAV with a fixed base. 
This wave is Rayleigh and its phase velocity can be approximated by the 
expression [3] 


srt 
} 


If 2 # 0, then, solving equation (8) with accuracy up to o( 4") by the 
perturbation method, we find vu --v,-|- Av (I -|- adjy,, Where 





—— 237)? = FI =) 
( 204 ’ ae | (| _— °=*) (1 — =*) , 
| ee. 2n**? . 1/2 t', 

= ene ee 
{(l WO) (i — 7=*)| ' 


Since the maximum value of & of a real isotropic solid is 0.5, then 
0.87 < 9 < 0.96 and accordingly a ®% 2.4. 


It was established as a result of studies that there are slight changes 
(on the order of ¢°’) of the direction cosines b'‘'®’ to the value 


| 9\—1/2 3 l cay!” ; 
yl) Ab - — (1—apy'’a, AB™ == — > (1-18") upon rotation of the 


substrate at angular velocity !) and there are changes of the phase 
propagation velocity of the PAV vo, proportional to 4, by the value 


(9) 








As follows from [4, 5], there is a correlation between variation of the 
phase propagation velocity of the PAV Av, on the one hand, and the 
frequency of the autogenerators and phase signal of the delay lines and 
filters on the PAV, on the other hand, which is described by the 
following relations: 


ahi aoe (10) 


where fo and $) are the autogenerator frequency and phase of the signal 
of the delay line or filter on the PAV at v = vo and Af and A¢ are 
variation of the frequency and phase, determined by Av. 


After substitution of expression (9) into relation (10), we find 


where L is the length of the operating section of propagation of the 
PAV. 


Accordingly, rotation of the base results in variation of the 
autogenerator frequency and phase of the delay lines and filters on the 
PAV. Thus, L = 2-10°2 m, vo = 5-103 m/s, Af = 2.9-10°9 Hz, and Ad = 

= 4.6-10°'3 rad at M = 19-2 “Shr. 


Thus, there are physical prerequisites for development of low-accuracy 
miniature angular velocity meters that use a PAV, with frequency or 
phase outcome. We note in conclusion that the sensitivity of these 
inertial data converters with the frequency layout is dependent only on 
the elastic constants of the substrate, while that with the phase layout 
is also dependent on its dimensions. 
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[Text] Known specialized applications program packages for analysis of 
the dynamics of mechanical devices [1, 2] may not be used to design 
gyroscopic devices with regard to a real vibration-insulation system and 
the complex control circuit, which is ordinarily nonmechanical. This is 
caused by the constraint on the relationships on bodies in the Lagrange 
method of second kind, selected as the basic method for formulation of 
the equations of motion. 


This disadvantage can be overcome in the MIT-1 applications program 
package, which constructs a mathematical model based on general theorems 
of dynamics with representation of the device in the form of an oriented 
graph of relationships [3]. The general diagram of the proposed package 
is shown in Figure 1. The MTT-1 package permits one to formulate 
symbolic nonlinear differential equations of motion of a gyroscope and 
to solve them by using numerical methods. 


The input data for design are data on the structure of the device, the 
layout of shock absorbers and their parameters, the inertia-mass 
characteristics of the bodies making up the device, and the moments and 
forces acting on the components of the device. This information is 
represented in the form of ,erators of a specialized description 
language. This language is «u.tlined in more detail in [4]. 
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=< A aad 
Figure 1 
KEY : 
1. Description of device in input language of package 
2. Formulation of equations of motion 
3. With digital coefficients 
4. With symbolic coefficients 
5. Write (read) model to archive 
6. Mathematical model of device 
7. Construction of canonical form of digital model 
8. Description of function of right side of symbolic model 
9. Calculation of eigen-values and vectors 


10. Correction of parameters of mathematical model 
11. Correction of parameters of numerical integration 
12. Substitution of numerical values of parameters 
13. Numerical integration 


The mathematical model of the device, which has the following form, is 
formulated from this description 


MX + NX + KX = Df() + QP (t, X, X) + ” 


-+ Ff, (i, X, x, x) + F(t, Xx, x, X), 


where \j[n, nj. Nin, nj, K[n, nj. Din, m], Q[n, mi] are matrices with symbolic or 


symbolic-digital coefficients, X{n] is the vector of variables, f(t) is 
the vector of external perturbation, acting on the body from the 


direction of the base, P(t, X, X) is the vector of control, external 
moments and forces, and F,(t, X, X, X), and Fo(t, X, X, X) are nonlinear 


9h 








functions that contain paired and triple products of variables, their 
derivatives and external perturbation with symbolic or symbolic-digital 
coefficients. 


This type of mathematical model is related to the fact that expansion of 
the matrices of direction cosines with grouping of terms by orders of 
smallness into a series is fulfilled automatically during formulation of 
the equations. 


Selection of linearized equations as the main equations for design is 
related to the need to study the dynamics of complex gyroscopic devices 
with small computer resources. 


Special procedures of analytical transformations above power series of 
the following type (the MIT-1 package) were worked out to construct an 
algorithm for formulation of the equations 


n 
L = NV ka'pry?, 
ii 


where kj is a digital coefficient, @), %k, and 7, are variables and the 
first and second derivatives, and p = 0, 1, 2 is the exponent. 


The procedures of analytical transformations were subsequently modified 
(the MIT-1V2 package) for the case when 


3 
k, = k T] QQ), 
QI 


where k’; is a digital coefficient, and Q, are symbolic parameters that 
describe the inertial-mass characteristics of bodies. The MIT-1V2 
package exceeds the first version of MTT-1 in its capabilities, but is 
considerably inferior in speed. 


Retention of two similar program systems for formulation of the 
equations is related to the need to solve two design problems: 
selection of the parameters of the device and design of an optimal 
vibration-insulation system for a given device. In the first case, 
symbolic equations of motion must be compiled for determining the 
dependence of the coefficients in (1) on the parameters of the device. 
The MIT-1V2 package fulfills this. In the second case, with known 
digital parameters of the device, one must analyze the different layouts 
of the shock absorbers, one must select their parameters, and one must 
change the law of external vibration. The MTT-1V2 package can also 
solve this problem, but as experience shows, it is preferrable to use 
MTT-1, which operates two-threefold faster. 
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Consideration of the interaction between bodies is the most origina] 
part of the package. 


There are no constraints on the type of links in the formulation 
algorithm, and they are taken into account by substitution of codes of 
the moments and forces into the corresponding pro,jections of the 
equations. A specific mathematical expression of the moments and forces 
is contained in special matrices or is formulated automatically by the 
package for description of the device. These expressions are 
substituted into the equations at the stage of constructing the 
canonical forms. This permitted one to work out special procedures of 
correcting the parameters of the mathematical model without changing the 
equations of motion. The parameters to be corrected include coordinates 
of the points of attachment of the shock absorbers, the coefficients of 
stiffness and damping of the shock absorbers, and external 
perturbations. 











I) (2) | oy 
SnenrentT npuGopa | fy. yo Tp ly la fpened 14.10 i 
Z)hopuve 6,5 6,5 7,0 
3h Lapa Knast! paMka 3,7 3,4 3,6 
5)Kooxwyx tupo- ; . 
O.10KA 1,9 2,5 2,4 
(6)Porop 1,4 1,4 1,0 
KEY: 
1. Component of device 4. Outer gimbal 
2. N-m? 5. Housing of gyrounit 
3. Body 6. Rotor 


A canonical form of the following type is constructed from the generate1 
equations 


Y - FY, g, (2) 


where Y is a vector of the variables of state of the system. An 
accurate expression of function F(Y, t) can be found for equations with 
digital coefficients. It is formulated in the MTT-1V2 package by a 
special algorithm at each step of numerical integration. 


System of equations (2) is subsequently solved by using numerical 
methods, which are selected automatically by control operators, which 
are a constituent part of the data preparation language. 


The MITT-1 package and its second version MIT-1V2 perform the entire 
cycle of study of the mathematical model from its formulation by 
description of the device to numerical solution on a computer. This 
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permits the designer, changing the study control operators, to conduct 
various numerical experiments with a mathematical model of the device. 


The dynamics of a cushioned gyroscope with angular vibration of the base 
with regard to the real vibration-insulation system is studied in this 
article by using the MTT-1 package. A gyroscope in a gimbal suspension, 
mounted on spring absorbers, was selected as the object of study. Let 
us solve the classical] problem of estimating the efficiency of the 
vibration-insulation system under resonance conditions. Analytical 
solution of a similar problem is given in [5, 6], which permits one to 
present comparative estimates. 


A layout of a gyroscope in a gimbal suspension and the orientation of 
the axes of the coordinate system correspond to that adopted in [6]. We 
assume that the axes of all the coordinate systems are the main axes of 
inertia of their own bodies, the origins of all the systems coincide in 
the ini*‘ial position, and the inner gimbal is rotated by 30° about its 
own axis. The moments of inertia with respect to their own axes are 
presented in the table. The natural moment of momentum of the rotor is 
equal to 5.4-10°! N-m-s. The overall] dimensions of the gyroscope are 
O.1 x 0.1 m and the axial stiffness of the springs is Cj = 1,617.7 N/m 


(i = 1, 4). Damping moments with coefficients 9.8-10°§ N-m-s act along 
the axes of the suspension. 


Let us consider the motion of the body only along axis z. The dampened 
coefficient along this axis is equa] to 1.1765-10°% N-m-s. 


Let us compile complete nonlinear equations of motion of a gyroscope 
with regard to the nonlinearities up to third order of smallness 
inclusively, which are recorded in the archive. Using the method of 
step analysis of the mathematical model, realized in the package, let us 
increase the complexity of the model for the entire cycle of studies. 


“reid M 


1h, ; 

anvinl 

0,24 | | / 
Vo 


0,12 
Wal0b} -. 
pee) 


Figure 2 
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We find the frequencies of natural vibrations of the system in the first 
step by the linear model (it is found from general model (1) by 
eliminating part of the terms and this exclusion is performed 
automatically. With the selected input data, v; = 208.12 s~! and py = 

= 1,033.67 s~!. 


Let us determine the amplitude of the angular vibration velocity of the 
outer gimbal of a cushioned gyroscope with main resonance at frequency 
vo, Let us first give the perturbation along axis z in the form @; = 

= 4 sin (1,033.67 t) using the correction section. Integrating the 
linear system of equations by the graphs which are constructed 
automatically, we find the desired value: A % 10 s~', which is in good 
agreement with the results of analytical calculation performed by the 
formulas of [6]. 


The last step in analysis of a cushioned gyroscope is solution of the 
complete system of equations with nonlinearities up to third order of 
smallness inclusively. The results of the calculation are presented in 
Figure 2 (curve 1). 


The equations of a rigidly mounted gyroscope with the same law of 
perturbation of the base must be formulated to analyze the efficiency of 
the cushioning system. We should compile a similar task for a design, 
having taken from it a description the flexible elements. By 
integrating the linear system of equations, we determine the amplitude 
of the angular vibration velocity of the outer gimbal: A» = 310 s°!. 
The results of integrating the nonlinear equations are presented by 
curve 2 in Figure 2. Calculations show that the drift of a cushioned 
gyroscope is much less than that of a rigidly attached gyroscope. It is 
essentially invisible at the selected scale. The ratio between the 
drift of gyroscopes mounted rigidly and on cushions coincides with that 
presented in [6]. 


The use of the MTT-1 package for study of the drift of a cushioned 
gyroscope considerably reduces the analysis time. 
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[Article by V. N. Fedorov, candidate of technical sciences, Kiev 
Polytechnical Institute] 


[Text] It was assumed during study of the azimuth motion of a 
ground-based pendulum gyrocompass (MGK) during acceleration of its rotor a 
[1-3] that the angular momentum of the sensitive element (ChE) increases 
by linear law. Linear law can be realized rather simply in the presence 
of a controlled drive of the gyromotor. The dependence of the angular 
momentum on time in the case of uncontrolled acceleration is linear only 
during the initial segment of the acceleration curve [4]. The azimuth 
motion of the sensitive element of a pendulum gyrocompass is considered 
in this article with exponent al law of variation of the moment of 
momentum, which describes the orocess on the entire section of the 
uncontrolled acceleration. The constant My and linear time-dependent 
M;t uncontrolled moments with respect to the vertical axis of the 
sensitive element are taken as perturbations. 


The system of precession equations that describe the motion of the 
principle axis of the MGK at small angles of deflection from the plane 
of the meridian has the following form upon variation of the angular 
momentum 


Her -|- HU, + (mgl + HU,) B = 0, 


a (1) 
Hf + HB — HU,a = M, + Mit, 


where H is the angular momentum of the gyroscope, mgl is the pendulosity 
of the instrument, U; and Us are the horizontal and vertical components 
of the earth’s angular rotational velocity, respectively, @ and j are 
the angles of rotation of the sensitive element of the MGK in azimuth 
and in the vertical plane, and t is current time. 
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It is easy to find from system of equations (1) the equation of motion 
of the sensitive element with respect to the azimuth coordinate. With 
regard to mgl » HU;, it has the form 


Ha + 2Ha. 4- mgiUa = — 2HU, — M, “84 — yt Met (2) 


Let the angular momeut vary by the law // =//,,—(H,,—/,)e-“, where Ho 
and H, are the angular momentum at the beginning and end of acceleration 
and 4 is the attenuation index of the exponential function. 


Let us introduce independent variable z= (1 —/,/17;')e™ (2(0) Zo = 


mt 


= 1--1,l//;,'). Denoting the derivatives of the desired function with 


4 respect to the new variable by a’ and a", let us rewrite equation (2) as 
follows: 
2?(z— 1) a" 4-2(3z — la’ —U,mgl7H a = 
== QU,A7™'2 4- Mymgik? Hi,” (1 — 2)! + (3) 
+ My mgih 7? (1 — 2)! In(z,27!). 
Assuming [5] 
Q= 2th ry) (2), (4) 
where p —U,mgih~*lin' and i is an imaginary unit, let us write equation 


(3) in the form 


2(z-- 1)U" + (21 Vp + 3)z—(21Vp + I)JU' + 
4-(2i Vp — p)U = 2U,h'2-!P 4 (5) 
+ mgth 711? (My + M, A In z,) 2-8-1 (= 2)! — 
— M,mginHp2-" P-1 (1 — 2)" In 2. 


A homogeneous equation, corresponding to (5), is a hypergeometric 
equation (Gauss equation) and has the solution 


(6) 
U (2) = CU, (2) -+ CU, (2), 


101 








where C; and Co are arbitrary constants and U,(z)=-F(a, >, c, 2) and 
Uo(z)==F(a—c+1, b—c-+-l, 2--c. z) are hypergeometric functions with 
parameters a=il/; b= i)/p-+2; c=2iyp +1. 


The partial solution, found by using the method of variation of the 
arbitrary constants, can be written in the form 


Usact (z) = 2U,47" [U; (2) Q; (2) — U, (z) Q, (z)] + 
+ mgik Hn My + M,A™" In 2) [(U, (2) Re (z) — U2 (2) R; (2)] — 
— M,mgihkH?, WU, (2) S,(2)- Ue (2)S,(2)), 


(7) 


where 








, U | 9 (2) dz U, » (2) dz 
Q . (2 _ 1.2 R = 1,2 . 
12 (2) \ 2' +4 (1 — 2) A (2) 1.2(2) \ aq “Sie (8) 


In 2U, 9 (2) dz 


S1.2(2) = | 2’ (1 —2z)?A (2) ’ A (2) = (1 —c)z (1 —2) er" 





The general solution of equation (5) consists of the sum of the solution 
of the homogeneous equation, shown in the form of (6), and of partial 
solution (7). Integration constants C,; and C» are determined with the 
initial conditions 


2 = 25: U (a) = Ug; U' (z) =U". (9) 


After simple transformations, the solution of equation (5) can be 
written as follows: 


U (2) = Uga~' (25) [U2 (29) Uy (2) — Uj (zo) Us (2)] — 
— Up" (29) [Ug (29) U (2) — U, (29) Ug (2)} + 
+ 2U gh" [Qs (2o» 2) Us (2) — Q, (25s 2) Up (2)] + 
+ mgik~ Hi" (My -- M,A~" In 2,) [Rs (29, 2) U, (2) — 
— Ri (2s 2)Us (2)] — MymginH?, [S, (zy, 2)U, (2) — S, (zp, 2) Us (2)], 


(10) 


Where Q0(z,,2), Rj,2 (22), Si,2(%, 2) are integrals of type (8) with 


integration limits Zo and z. 
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Having determined on the basis of expression (4) the analytical 
correlation between (9) and initial conditions (0) = ao and f(0) = Bp 
of system of equations (1) 





- —a , = melBo “a ome ~ (Iba) | 
U,=%2z5%, Us= la riz) + Akg | 25 


let us write the desired solution of equation of motion (3) of the 
sensitive element with respect to the azimuth coordinate in the form 


© (2) = a1g2"Z9“A~' (25) [U2 (2) U; (2) — U; (2,) Ug (2)] + 
+ 2U,A~'za [U; (2)Q, (2p, 2) — U;(z)Q, (2, 2)] — 
z\a | ml _ Us 
- (é) 7A ee) E i ti —%) h a| x 


m 





X [Us (Zo) U3 (2) — U; (29) Us (2)) + 


-f- men (/ (M 0 += mt —+In 2) [U, (z) R, (Zp, 2)— Us (2) R, (2 z)) a 


— M,nglz°n Hy,” (U, (2) Sq (2q, 2) — Us (2) S; (Zo) 2)). 


The expansions of the hypergeometric functions and the integrals from 
them, for example, with respect to Chebyshev polynomials [6] and also 
direct integration of the equations of motion on a computer can be used 
upon numerical study of the motion of the MGK in the mode of an 
exponential increase of the moment of momentum. 


H,10"? H-M-c (1) 
a, 1 10° pad (2) 
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The results of integration of system of equations (1) at mgl = 2-104 
g:cm-s, U; = Us = 5.161-10°5 s“!, §(0) = 2.58-10-5 rad, a(0) = 5.2-10-? 
rad, Mo = 1.5-10°4 g-cm, and M; = 2.5-10°' g-cm-s~! in the interval t = 


= 0-700 s are presented in the figure. Solid curve 1 shows the 
dependence of the angle of rotation of the sensitive element at azimuth 
on time in the case when the angular momentum varies by exponential law 


(curve 2). The dashed curves show similar functions (@ = a(t) and H = 

= H(t)) for the same initial conditions and perturbations when observing 
condition h = (Hp - Ho)A, where 4 is the characteristic slope of 
acceleration of the rotor in acceleration law H = Ho + ht. 


It follows from the figure that the functions found in [1-3] can be used 
on the initial leg of the acceleration curve where it is described 
rather well by a linear function, upon study of the azimuth motion of 
the sensitive element of the MGK. The entire leg of the exponential 
accleration of the rotor is determined analytically by expression (10) 
of this article. 
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[Article by A. S. Apostlyuk, candidate of technical sciences, and V. N. 
Sheludko, candidate of technical sciences, Kiev Polytechnical Institute] 


[Text] Those numerical methods which permit one to integrate the 
equations of motion at a high rate can be used to construct effective 
discrete models of motion of mechanical systems, and the requirement of 
high integration accuracy is not compulsory for these methods. 


Let us construct a discrete model of motion of a mechanical system, 
consisting of an absolute solid, connected to base N by shock absorbers. 
Each shock absorber is represented in the form of three mutually 
perpendicular elementary springs, articulately attached to the base (see 
figure). Following [1], let us assume that the springs are arranged 
along the principal directions at any moment of time, and that 
displacement of the body, perpendicular to the axis of this spring, does 
not alter its reactions. 
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Let us introduce the conditionally fized coordinate system Of9¢{ and 
Cxyz, invariably bound to the body, and the origin C of coordinate 
system Cxyz is compatible with the center of mass of the body, while 
axes Cx, Cy and Cz will be considered its principal central axes of 
inertia. Let both coordinate systems coincide at the initial moment of 
time. Let us denote by a;;, Jij, and 7;j the cosines of the angles 

rm formed by the axis of the j-th spring (j = 1, 2, 2) of the i-th shock 

| absorber according to axes Of, O7, and Of. Let us select coordinates é, 
7, ¢ that characterize the translational motion of the body with respect 
to the base and Euler angles #, @, y that describe the rotation of the 
body about the center of mass as the generalized coordinates. 


In this case, the law of motion of the solid in the flexible suspension 
can be determined by using Newton’s second law and Euler-Poisson 
equations: 


me = F mij = F,; mt = F,; 
1104 — (Ig -— 15) yo, = ANS Loy — (Uy — 10,0, = Mo: 


(1) 


1400, — (Iq ~ I) 00 = Mb 


. 


cos (p M . 
t My; 0 = cosqw, — sin (y 


= -—— @. -|- —— 0." 
sinQ * | sin U y’ 


© = — sin pclg 0w, — cos pclg Voy + w,, 


where m is the mass of the body, I;, Io. and Ig are the moments of 
inertia of the body with respect to axes Cx, Cy, and Cz, respectively, 


=> 


r= (EW 0) is the vector of the absolute accel=ration of the center of 


mass, F;, Fo, and F3 are the projections of the main vector of the 
external forces F onto axes Of, Oy”, and I, ard ¢x, wy, and w; and M’,, 
M’o, and M’3 are projections of the vector of the instantaneous angular 


_~ 


velocity w and of the principal moment of external forces M onto the 
axes of fixed coordinate system Cxyz. 


_ 


The principal vector F is the geometric sum of reactions of the shock 
absorbers and of the gravity of the body, while the principal moment M 
is the sum of moments of these forces with respect to the center of 

— 
mass. The reaction Rj of the i-th shock absorber is a function 
(generally nonlinear) of deformation of u, and of its first derivative 
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uj. Let us represent it in the form of [1] Rj(u;, ui) =F (ui) Us (03). 


where Fj is flexible and U; is dissipative forces, and let us consider 
the order of calculation of the principal vector and of the principle 


moment in equations (1). If, for example, F;;= —k,ju;;, U;j=0 (i—1,... \: 


j = 1, 2, 3), where kjj is the stiffness coefficient of the j-th spring 
of the i-th shock absorber, then the principal vector can be determined 
in the following manner: 


N 3 
| % 
Fe—\ > 

i=| j=l 


Re, teyj. 


Assuming that the stiffness coefficients kj; and coordinates xj, yj, and 
zi; of the points of attaching the shock absorbers to the body (the 
latter are invarable in the bodcy-axes system Cxyz) are known, 


displacements A= (On, 5i2, 5i3) of these points with respect to the base 


is determined by vector equalities (see figure) 


— > —> 


6=M+Ar,—Pip (i=1,...,N) (2) 


where 9.=(19 is displacement of the center of mass with respect to 
the base, /; = (X;,y;,2;) are the constant radius vectors of the points of 
attachment in coordinate system Cxyz, Pin = (Eigs Nin Gin) are the radius 


vectors of the points of attachment at initial moment of time in 
coordinate system Of7(, and A is the matrix of transition from 
coordinate system Of7( to Cxyz. 


Deformation uj; is equal to the scalar product of displacement 6; by 
unit vector Dy = (ai; Bis Vis) » Airected along the axis of the spring: 
Uj; = 8, pi. Assuming that the principal directions of the shock absorbers 
are parallel to axes Of, O7, and O(, we write the projections of the 


~ 


main vector F onto the axes of fixed coordinate system Of 7( 


N 
Fy=— Yikiby (jf = 1,2, 3). (3) 


i=} 


107 








The expression 


. (4) 
Fi; = A'Fy; 


determines the elastic force of the ij-th spring (F’;j;) in coordinate 
system Cxyz (the superscript "T" denotes the transposition operation, 
while the prime indicates the assignment of the vector marked by it in 
the fixed coordinate system Cxyz). Finding the moment of this force 
with respect to the center of mass C by the formula 


Mi, =X Fip (3) 


we find the principal moment from the expression 


we=v YP XF... (6) 


Having been given in one way or another the forces and moments of 
figuring in the right sides of equations (1), for example, by formulas 
(2)-(6), one can generally integrate equations (1), using any of the 
standard numerical methods of integration of the system of ordinary 
differential equations. However, this approach is hardly suitable in 
real-time simulation. It may be feasible to use difference 
approximation, based on the principles of discrete variational 
calculation [2, 3]. Let us assume that equations (1) can be written by 
Lagrange functions L = T - Il (T and Il are the kinetic and potential 
energy, respectively), which is dependent on the generalized coordinates 


(their vector column q has the form . = [&, 9, ¢, # 8, p)), on 


velocities q and time t: 


d /oL OL 
— (|) —-=— =0. 
dt (F) 0q (7) 


It is known that equation (7) is a condition of the stationary state (of 
t 
the approach of the first variation to zero) of the integral *, ldt at 


the selected constant quantification step h = (t» - t;)/v (v is the 
number of intervals), the integral can be approximated 
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\’ L(t dy 9,) 4 = > Lyf, (8) 


p=! p=! 


where ty is the moment of the end of the #-th interval and qy, is the 
value of the vector of generalized coordinates at this moment of time, 


while the vector of generalized velocities a, is determined as follows: 


I= (Gy— Gh (w= dy... ¥). (9) 


t 
Sum (8) will approach symbol | *, Lat upon natural assumptions with 


respect to function L. Accordingly, a discrete analogue of (7), i.e., a 
difference relation that guarantees approach of the first variation of 
sum (8) to zero, will be the relation [2] 


1; / OL, dl. dL 
(<3) —-==0, p=l,..,¥, 


Oyj, Oy an (10) 


which, together with formula (9), forms a discrete model of the dynamic 
system described by equation (7). Since the method of constructing the 
Lagrange functions for dissipative systems is sufficiently complicated 
(see, for example, [4]), the dispersion of energy in the system can be 
considered in the first approximation, having adopted the value 
exp(-nt)L as the Lagrange function, i.e., having approximated equations 
(1) by the following discrete expressions rather than by relations (10) 
[5]: 





Om, (11) 


where n is a parameter that characterizes the damping in the system. 


Having used expressions (9) and (11), let us construct the desired 
discrete model of motion of a solid on a flexible suspension. To do 
this, let us determine sequentially the potential and kinetic energy of 
this mechanical system. 
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The potential energy of the system is: I]=[I]p-+I]p, where fe = 


N 3 
] 
= >» »> Ry jut, is the potential energy of 3N deformed elementary springs 


and IL, = mq{ is the potential energy of the body in a gravitational 


Let us consider this arrangement of the shock absorbers when the axes of 
the springs are parallel to axes Of, O7, and Of. Then 


N N N 
I .| dll orl 
= = hab r= Y babe Ge =); Miadie + og. (12) 


i=] i=] i=] 


We note that finding the partial derivatives of the potential energy 
with respect to coordinates #, @, and g is difficult, since it requires 
that the same derivates of matrix A be found. Therefore, taking into 
account that the derivates of Il for these coordinates are moments of the 
elastic forces of the shock absorbers with respect to the axes of the 
corresponding rotations of the body, we write the following equalities: 


am a os 3 cos 0; 
5p = Misin psin® + Mzcos@ sin 8 + Mscos 0, (13) 


dll 


’ es. oll ° 
3 = Micosp—Mrsing; 35> = M;. 


The partial derivatives of potential energy with respect to generalized 
velocities are equal to zero. 


The kinetic energy of the system can be represented in the form 


an +e + n? +) 4 + (Iw? + Io! + 1,2), 
and with regard to the relations 
oO, = sin 0 sing + 8 cos @, oO, = psin @cos p — Osing, 


, _ ip cos 6 +9, 


one can write 
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T => @ +- n? + ¢?) + > [/, (p sin ® sin @ + § cos @)? + 
| (14) 
+ 1, (xp sin @cos p — Osin @)* + J, (th cos 6 + @)?}. 


Further postulation of expressions of the kinetic and potential energy 
into formula (11) results in implicit difference equations of numerical 
integration. This can be avoided in the case of equality of moments of 
inertia I; and I9, which is found rather frequently in real designs of 


vibration protection systems. Then, differentiating (14), we find the 
following expressions: 


OT _ OT  , OT _ oT oT ss , 
EO Fy = Oa =O a” mn — 
aT oT 2 ei +4 ' } ll 
e = 0, =, = /,p’sin6 cos 0 — /, (p + pcos 8) psin 8, a = 0; (15) 


Tw Isp sin? 0 + /, ( + wp cos 8) cos ®, <= 1,8, 
ay 0 


or = 1, (9 + tp cos 8). 
dp 


After substitution of relations (12), (13) and (15) into (11) and also 


after certain transformations, we find the desired discrete model of 
motion of a solid in a flexible suspension: 


_anig fh Om 
B41 =ennh (E ~~ Th %) ’ G41 = F + hE is); 
A oll be 
nm an, ) Na = m+ hy, 
an (ph OU 
C41 = eh (6 _ mn) » Sear = Cy + Ala; 


ée 


1, sin* 8, [/;1p, sin’ 8, + /, (9 +} i cos 8,) (cos 8,— cos 8,41) +- 
oll oll 
+h (ae, cos 8:4; — 5y,)] , 


Vis = tp; + hip; 


O41 = e-"" {8, + ; If — 1) 2 cos , sin 8, — 





V4 = 
(16) 


ee . oll 
— 1,0, sin ® — 3, | 
s HP, sin 0, 00, 


014.1 = 0, + 10.413 


lil 








: : ere pe : h ail 
Py, = — Pi cos B14 + e-™ (7 + tp, cos 8, — T39,) 


Pi4t =; + AQ, 4 


We note that the expression I;sin"6\,;, which brings greater errors to 
the calculations at small] angles @, while it makes them possible at 


6 = 0, is found in the denominator in the difference equation for Piel: 
This can be avoided, having assumed that the movable coordinate system 
Cxyz is rotated at initial moment of time with respect to the fixed 
coordinate system by angle @, = 90° [6]. 


As studies showed, the use of difference equations (16), based on the 
principles of discrete variational calculation, is more effective in 
speed for simulation of the motion of a solid on a flexible suspension, 
described by equations (1), than the use of numerical integration of 
equations (1) by the Runge-Kutta method for the same purpose. 


BIBLIOGRAPHY 


1. Kolovskiy, M. Z., "Nelineynaya teoriya vibrozashchitnykh sistem" 
[Nonlinear Theory of Vibration-Protection Systems], Moscow, 1966. 


2. Logan, J. D., "First Integrals in the Discrete Variational 
Calculus," AEQUAT. MATH., No 19, 1973. 


3. Apostolyuk, A. S. V. B. Larin, and V. N. Sheludko, "On the Problem 
of Constructing Discrete Models of Mechanical Systems," 
MATEMATICHESKAYA FIZIKA I NELINEYNAYA MEKHANIKA, No 2 (36), 1984. 


4. Zaplatnyy, V. I. and N. Ye. Tkachenko, "Conetruction of the 
Lagrange Function for Dissipative Systems Br sed on Solution of the 
Inverse Problem of Variational Calculus," DOKLADY AKADEMII NAUK 
UKRAINSKOY SSR. SERIYA A, No 4, 1984. 


or 


Appel, P., "“Teoreticheskaya mekhanika. T. 2. Dinamika sistemy. 
Analiticheskaya mekhanika"” [Theoretical Mechanics. Vol 2. 
Dynamics of the System. Analytical Mechanics], Moscow, 1960. 


6. Sheludko, V. N., “On Digital Simulation of the Motion of a Solid,” 
"Modelirovaniye i upravleniye dvizheniem sistemy tverdykh tel" 
[Simulation and Control of the Motion of a System of Solids], Kiev, 
1984 (Preprint 84.6, AKADEMII NAUK UKRAINSKOY SSR. INSTITUT 
MATEMATIKT) . 


112 








UDC 007.62 


ALGORITHM FOR REAL-TIME CALCULATION OF REDUCED MOMENTS OF INERTIA OF 
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and V. S. Yevgenyev, candidate of technical sciences, Kiev Polytechnical 
Institute ] 


[Text] The reduced moments of inertia of the links must be calculated 
on a computer to control the motion of complex manipulation 
anthropomorphic-type systems with large number of degrees of freedom 
{1}. In view of the fact that the memory capacity and speed of modern 
computers do not permit this operation in real time, the problem of such 
equivalent representation of the links of a manipulator as solids 
arises, which would simplify considerably the solution of the postulated 
problem. The purpose of the article is to outline a method of 
calculating the reduced moments of inertia of the links on mini- and 
microcomputers using construction of a system of material points, 
equimoment to arbitrary solids. 


An equal moment system of the fewest number of points, which is equal to 
four, can be constructed to solve the considered problem [2]. The equal 
moment nature is guaranteed by the coincidence of the centers of inertia 
and by the equality of mass and tensors (ellipsoids) of inertia. If the 
principal central axes of inertia x, y, z are selected as the reduction 
axes, it is sufficient to require coincidence of only the axial moments 
of inertia. These conditions are expressed analytically by nine 
equations 


1 4 . 
q mM. 
y myx? = A, y my? = B, y m2? = C; 
1 I 1 


‘ ‘ 7 
YS mex = 0. Y myer = Dp mater = 0 1) 
1 j l 


4 4 e 
»> mx, = 0, y my, = 9, y) myz, = 0, 
1 ! I 
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where A, B, and C are the moments of inertia of the link with respect to 
coordinate planes yz, xz and xy, respectively. 


The desired system of four material points should belong to a 
second-order surface, which is an ellipsoid of inertia. Let us use the 
equimoment ellipsoid of inertia of the Legandre type, described by the 
following canonical expression, as this ellipsoid [2] 


2? l 
Tw" 7 ’ 


tls 


x? 
> es 


where M is the mass of the manipulator link. 


Since it is difficult to determine directly the coordinates (x;, y), 
Zi)y e+ey (X4, V4, 24) Of the four desired points, let us transform the 
coordinates according to the equalities 


_A B C | ; 
“= 77 Mp Y= HF Pw 42> 4 (¢ = 1,2, 3, 4). (2) 


If one assumes that the masses of the four material points are identical 
and equal to M/4, one can write equations (1) with regard to (2) with 
respect to the new coordinates #;, pj, qi in the form 


. (3) 
y Hip; = 9, p rig = 9, »> Hig, = OV; 
1 


4 4 4 


YH = 9, Ym =9, 24 = 0. 


It follows from equations (3) that the ellipsoid of inertia of the four 
material points is a sphere of unit radius. To simplify the calculating 
procedure, let us select the desired points at the vertices of a right 
tetrahedron, lying on the sphere. 


It is insufficient to determine 12 coordinates sj, pi, qi (i = 1, 2, 

3, 4) of a system of nine equations (3). Taking into account that the 
axes of symmetry for a homogeneous solid are the principal central axes 
of inertia and by locating one of the ,oints on axis z, we find 
additional relations: 
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h#=—%=0; g,=V3. (4) 
We note that if the solid is homogeneous, conditions (4) vary according 


to the geometry of the mass. 


Solving »lgebraic system of equations (3) with regard to conditions (4), 
we find 











0 
ty 0 I's 7 
2)+(3)(3)-f77} 
V1 V3 qo _ ti 
V3 | 
V3) V3} 
Its at Hy {7 
(»)- —| + ’ (») — | 5 
qs _ 4 | 
V3 3} 











| SH) 2 foo / DB oC) 
5 Ties he U; mM J » Ig >> 9; 2 |/ 3M : | 5i97| 
fy 2A OB. 4/ Ce | (5) 
r= 3 —-Var—h BAT I 
Jp / 2a (WB. _ |/ © 
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Formulas (5) permit one to determine the coordinates of the desired 
system of four point masses, equimoment to the manipulator link in the 
principal central axes of inertia x, y, z, bound to it. One can turn to 
the physical axes of the suspension {, 7, ( by parallel transfer of 
coordinate system xyz to values £9, 90, (09 and by rotation by Euler 
angles a, 8 and 7. The reduced moment of inertia of the link with 
respect to the axes of its suspension can then be calculated by the 
formulas 


Ty = [yy + yy + 1405, + ANE: 
In == Uy + Tyg + Lattin + MM 


Ie, = 140, + Ty, + 1,03, + MG. 


(6) 
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Here ajj are the direction cosines that determine the angular 
orientation of axes x, y, z and that are dependent on the method of 
being given the Euler angles [3] and Ix, Iy, and Iz are the moments of 
inertia of the link with respect to the corresponding axes. 


The relationship between these moments of inertia and A, B, and C used 
earlier is determined by the identities 


A = 0,5(/, +-/,—1,); B= 0,5 (/, -+ I. — Ty); 
C= 0,5(/, +- ly i). 


The derived formulas permit one to calculate the reduced moments of 
inertia of multilink manipulation systems according to the following 
algorithm. 


1. Find the mass M of links, the position of their centers of inertia, 
and the distances fo, #0, and (9 from these centers to the suspension 
axes of the links. 


2. Assign the Euler angles a, #, and 7 and calculate the direction 
cosines 4;; for each link of the manipulator. 


3. Determine the axial moments of inertia of links Ix, Iy, and Iz and 
find A, B, C by formulas (7). 


4. Calculate the coordinates of four points by mass M/4 for each link 
by formulas (5). 


5. Find the moments of inertia of equimoment systems of material 
points, which replace the links of the manipulator, with respect to the 
principal central axes of inertia. 


6. Calculate the reduced moments of inertia of the links with respect 
to their suspension axes by formulas (6). 


Computer programs for calculating the reduced moments of inertia of the 
simplest type of links--parallelepiped, cylinder, truncated cone and 
triangular pyramid--were worked out according to the proposed 
algorithms. The computation time on the SM 1810.41 microcomputer is 
2:10-% s, which indicates the effectiveness of the algorithm. 
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{Article by B. I. Genkin, candidate of technical sciences, N. I. 
Nagulin, candidate of technical sciences, and A. G. Arkhipov, senior 
scientific associate, Severodonetsk Branch of Kharkov Institute of Radio 
Electronics] 


[Text] The timely problems of modern instrument building include the 
development of devices, capable of operating stably at low-frequency and 
infralow-frequency, for example, seismic effects. One of the promising 
ways of solving this problem is to use vibration insulators with 
stiffness correctors [1]. These correctors permit one to optimize the 
dynamic elastic characteristic of vibration insulators, without 
influencing their static capacity. Deviations of the parameters of the 
stiffness correctors from optimal values are inevitable in view of the 
different types of random factors. Therefore, real dynamic elastic 
characteristics of vibration insulators with stiffness correctors are 
always nonlinear. The given paper is devoted to analysis of some types 
of nonlinear vibration-protective systems with stiffness correctors. 


Let us consider a system with one degree of freedom, consisting of a 

vibration-protection object of mass m, attached to a base, vibrating by 
law h(t) = ho sin wt by parallel-connected main spring of stiffness c, a 
viscous damper with viscosity coefficient gs and stiffness corrector, the 


characteristics of which are: Ce is the stiffness of the corrector 


spring, d is the maximum deformation of the corrector spring, and L is a 
geometric parameter. Let us take the displacement x of the vibration- 
protection object with respect to the position of equilibrium as the 
generalized coordinate. 


Let us introduce the dimensionless coordinate £ = x/L and dimensionless 
parameters @ = s,/s and f = d/L. The dynemic flexible characteristic 


F(x, c, Cy d, L) of the vibration insulator, consisting of the main 
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spring and stiffness corrector, is represented in the form F(x, c, c,, 


kK 
d, L) = cLf(&, @, 8) where f(£, a, 8) is the dimensionless dynamic 
elastic characteristic of the vibration insulator. Let us compile a 
differential equation with respect to vibrations of the vibration- 
protection object and let us reduce it to the form 


E+ Qnk + Rf (E, a, B) = jo* sin of, (1) 
where p= p/2m; k= Vel; y =h,!L. 
To find the approximate steady state solution of equation (1), let us 


use one of the linearization methods [2, 3]. Let us find this solution 
in the form 


E = Asin(wt — 9). (2) 
Let us replace nonlinear function f(f, a, #) by the approximate function 


i (Ea, B) © 9°(A, a, BE, (3) 


where q?(A, a, 3) are linearization coefficients. Having substituted 
(2) and (3) into (1), after transformations, we find 





A = yp*/V(q?(A, a, B) — p*}? + 48%p*; (4) 


» = arctg [25p/{9%(A, a, B) — p* ll, (5) 


where p = #/k and 6 = n/k. 


Relation (4) is an equation with respect to amplitude A and describes 
the amplitude-frequency characteristic (AChKh) of the system. Relation 
(5) describes the phase-frequency characteristic (FChKh) of the system. 


Let us consider two types of stiffness correctors: crank (Figure 1, a) 


and cam (Figure 1, b) [1]. The dimensionless eiastic characteristic of 
the vibration insulator containing these correctors is: 


[(§,a,6) = {1—a+a(l—py/Vi — EAE. (6) 
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Having used Ya. Z. Tsypkin’s formula [2], we find the linearization 
coefficient: 


g*(A, o,f) =1—a+ Fa(l—py4— Ay +—ay"4}, (7) 
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A computer program that realizes with step by frequency parameter p the 
search for the roots of equations (4), where q?(A, a, #) has the form of 
(7), with subsequent calculation of the values of ¢ and also of the 
amplitudes and phases of the absolute vibrations of the vibration- 
protection object, was worked out to study the characteristic features 
of the considered nonlinear systems. Analysis of the results of 
numerical experiment, conducted by using this program, showed the 
following. Equation (4) can have only one real root at small values of 
parameter 7, which characterizes the intensity of the exciting effect, 
and the amplitude-frequency characteristic of the vibration-protection 
object hardly differs from that found when using linear theory. A range 
of values of p, at which equation (4) has up to three real roots, 
appears with an increase of 7, the amplitude-frequency characteristic 
assumes a form typical for nonlinear systems: the resonant peak (curve 
1, found at 6 = 0.1 in Figure 2) is deformed, and additional branches of 
the amplitude-frequency characteristic (curves 2 and 3 in Figure 2) 
appear, to which the vibrations with large amplitudes correspond. Study 
of the stability showed that curve 2 corresponds to stable modes, while 
curve 3 corresponds to unstable modes of induced vibrations. 
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Figure 2 


With an increase of parameter é which characterizes the friction in the 
system, additional branches of the amplitude-frequency characteristic 
have a tendency toward fusion and displacement to a higher frequency 
region. Figure 2 was constructed for cases when the amplitude of the 
induced effect is independent of frequency (7 = const). Under real 
conditions, the amplitudes of higher harmonics of the kinematic induced 
effect is ordinarily much less than the amplitudes of several of the 
first harmonics. Thus, displacement of additional branches of the 
amplitude-frequency characteristic to a higher frequency range means the 
elimination of the danger of the corresponding large-amp]litude 
vibrations occurring [3]. Curve 4 in Figure 2 was found at é = 0.2. 
There are no additional branches of the amplitude-frequency 
characteristic. At the same time, an increase of friction results in 
some deterioration of the quality of vibration protection in the working 
frequency band. Accordingly, there is the problem of optimization of 
parameter é. 


A system with crank stiffness corrector on a low-frequency vibration 
bench, which permits one to increase the vibration amplitude of the 
vibrator to 15 mm, was tested for an experimental check of the results. 
The system to be tested contained four air dampers with replaceable 
pistons. Energy dissipation in the system was varied by mounting the 
pistons with openings of different area. Stable vibrations with large 
amplitude, corresponding to the positive branches of the amplitude- 
frequency characteristics, were observed at low friction. There were no 
nonlinear effects at high friction and also at small vibration 
amplitudes of the vibrator of the vibration table. 


Studies showed that the possibilities of large-amplitude vibrations in 
the working frequency band, which occur due to nonlinearity of the 
elastic characteristics of these systems, must be taken into account in 
design of vibration-protection systems with stiffness correctors. To 
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eliminate the undesirable nonlinear effects, one should try to see that 
parameter 7 have a rather small value, while parameter ‘ have a rather 
large value. A decrease of 7 at given ho is possible only due to an 
increase of L, i.e., due to an increase of the overall dimensions of the 
stiffness corrector. An increase of é results in deterioration of the 
vibration protection quality in the working frequency band. Thus, 
selection of the parameters of vibration-protection systems with 
stiffness correctors is an optimization problem. 
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[Text] Calculation and design of modern floated gyroscopic devices are 
related to postulation and solution of convective heat transfer problems 
in a working liquid. Typical forms of the suspensions of the sensitive 
elements are a cylinder in a cylinder, and sphere in a sphere; 
therefore, the desire to maintain the greatest number of symmetries in 
the design results in study of axi- and central-symmetrical problems. 
The basis for using numerical-analytical methods in solving these 
problems is inforamtion about the existence of the solution in general 
and about its behavior at different moments of time. 


Theorems on single-valued solvability as a whole and on the asymptotic 
stability of the solution of a system of quasi-linear evolutionary 
equations that describe the thermoconvective dynamics of a viscous 
incompressible liquid, which fills an axisymmetrical cavity in a solid 
coaxial to it, are presented in this paper. The problem to be studied 
is formulated in terms of the current #-temperature 6 function. 

1. If the liquid fills a cavity, which is the figure of revolution of a 
rectangle 7-={(mz)|U<R<r<R, fyz<ito} about axis 0z in a solid 
having the shape of a figure of revolution about axis Oz of triangle 
== ((r, 2) )0<r<Re. Hy<e<Hy(Ry<Ro, My<y<hg<l2)}, then the convective 


heat and mass transfer in it is described, as is known [1-3], by a 
system of equations 


Dw, — rd (ap, r-2Dip)/d (r, 2) —- vD*p = — yrO,, (7, 2,0 07; (1) 
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9 [0, —r-'0(p, @V/A(r, 2] — 2A = 0, (rz, NEU (2) 


with initial conditions of type 


P= Wolr2) (r.2)€n; @=@,(r,2), (r.2)€: 10 (3) 
and with boundary conditions 
b=Dp=0, (r,2)€dn; @=0, (r,z2)€OTl; t>0; (4) 
[0] = [Ad@/dn] = 0, (r,2)€dn; *>0, (5) 
where V* = 1~'0/0r(rd/dr) + 0/02"; D == rd/Or (r—'0/Ar) + 0%/d2*, D* ex Dp: 8(9, X)/a(r, 2) = 


dg/dr AX/z — Opldz AX/dr; n7=nx [0,T]; 11" = [xU (In) x[0, 7]; and n is the 
external normal to the boundary dx = 4/x; [rT] is a jump on the @z 
function t(r, z) € ff, equal to (t"—+’)|,, (the constrictions 7 by 7 and 


ll, respectively, are denoted by 7’ and r"), v, 7, A’, A", o” and o’ are 
positive constants that characterize the physical properties of the 
liquid and solid (j’ # A", o’ #0¢"), and t is time. 


It should be noted that rectangles were selected as 1 and Il in problem 
(1)-(5) only for clarity. The theorems presented below are also valid 
when 7 and Il are arbitrary constrained areas with rather smooth 
boundaries and 7 lags behing axis 0z by RD O. 


Conditions y = D# = 0 are not equivalent to classical conditions of 
adhesion of the liquid to a solid wall dz (these conditions are 
discussed in detail in [3]). The conditions of adhesion will be 
fulfilled if problem (1)-(5) is altered somewhat, having replaced 
equation (1) and condition (4) by the equalities 


D (r-*Dip,) — D [r-'0 (1p, r-2Dp)/0 (r, 2)] — D(r-2D*p) = — D(r-'®,), 
Plan = Mp/On|,, = Dip |,, = 9. 


The single-valued solvabiiity as a whole (with respect to t) of this 
problem is proved by the same scheme as problems (1)-(5). 
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2. Let Lp(7) be a Banach space, consisting of all functions y(r, z), 
measured on 7, and added with respect to 7 with superscript p > 1; Lp,(7) 
is equipped with evaluation 


Pll = ({ |p|’ rdrdz)'” (Il Iles = Il “lla 


Following [2, 4], let us introduce Hilbert spaces L, (x), L3 (I), W8 (x) 
W733 (1), WR (x), WY(n) With scalar products 


(9, )), = { pprdrdz; (8,t), = { o@trdrdz; 
i 1 
(PP) .2= | (vey) rdrdz; (8,1) 4 = \ o*yOytrdrdz: 


(PV) .0 = \ r—*Depr—Diprdrdz; 


(P P)3).2 = \ v (r—2DQ) y (r—2D yp) rdrdz 


and which correspond to valuations |j+\\,, |I+Ilj) II. (é = 1,2, 3). 


The space W{'}(n) (i= 1,2) consists of all elements of the Sobolev space 


Wio(x), which have zero sign on gx, while W‘3)5.9(7) consists of all 
functions g of space W%>o(7), for which the sign Dy is equal to zero on 


dx together with sign » (since R > 0, the valuations [Illi (= 1, 2,3) 


are equivalent to those of the corresponding Sobolev spaces). Space 
w' ')9.9(") consists of all elements of space W!'»(fl) with zero sign on 


Let us denote by V2i(n’), Viio(n’), Vib(n7), Vi8(n") (VE2q17)) the Banach 
spaces of functions p(r,z,f) (0(r,2,4)) from C°(0, T; Ly(m))NL,(0, T; WS$'} (x); 
VEO (x7) W5(0, T; Le (x); C(O, T; WE (m))NW2(0, T; W29 (x); C°(0, T; W328 (nm) N 


N L,(0,7; WS (m)) (C°(0, T; LO(U)) NL, (0, T; w3) (1l))) respectively, 
supplied with valuations 


‘ = max Ts 
IP lyso¢ar) ont Pilla +llvvlll, 
= ) . 
| \—p Ivitinty | Ivt.2¢n7) I 1; lle? 


I lvzdeaty = |Plysgeaty + MAX | lhayn +I Ve Mar 
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T 
1/2 
leagenty = (uae IV lnn + {I ly 9d) 


(|e Iy1.0n17) = nex lls + IV gr) 


T 
Where ||p|ir=(p qr (9, War = | (@ )), 44; H,(0,T; Hy(x)) is an ordinary 
0 


space [4] of distributions from [0, T] in He(*) with properties of space 
H, (with respect to t). 


Finally, let us denote by V3.1725 9 (42) the space of functions g(r, z, 
t) from V3-0,.9(#-), for which 


T—h 


ho! J Ilr, 2,¢+ h)—p(r,z,0) Ifo) 2x Ut gaa > 0. 


The spaces V39(II"), Vi(1I7) and yii/(q17) < yi®qyt) are determined in 


similar fashion. 


3. Let us call the pair of functions {#, @}, belonging to space 
V2 (07) x Vo'0 (II) and which satisfies the following integral identities the 
generalized solution of problem (1)-(5) 


| rlpDydede| “* + \ |r '@,Dp— 99 (Wp, DWV, 2) + 


(= 
n yom 


(6) 
+ vr~' \yDwpy¢@)] drdzdt = ) Oq,drdzdt; 
al 


{ a0 rdrdz|—* 


II 


4+- | [—ro0t, —od (wp, 0)/0(r, 2) + 


t 
{=0 nt (7) 


+r (y8yt)] drdzdt = 0 





Vt, from [0, T] and arbitrary pair {y, 7} from Vo'io(n") x Vi (I"). 


It is easy to show that the classical solution of problem (1)-(5) 
satisfies relations (6) and (7) and that the determination is correct, 
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in other words, all the integrals in (6) and (7) have meaning (since 

R > 0, i.e., the axis of symmetry does not belong to a liquid-filled 
cavity). 

Theorem 1. If the initial data of problem (1)-(5) are such that 

tba € WE (20), (,€L,(I1), it has a unique generalized solution (wp, 0)¢€ V2? (17) y 
x yyl/2(") at arbitrary final solution T. Moreover, the following 


integral relations are valid for (#, @) 


Qo ( V*rdrdz _ ++ v | w*rdrdzdt = y { Ou*rdrdzdt, 
at 


nf nt 


~" \ | w/r |?rdrdz _ v \ |v (w/r) |? rdrdzdt = y ) Od/dr (w/r) drdzdt, 
t nt nt 
go! | o0*rdrdz |\=) + + i. (v0)? rdrdzdt = 0 
il 


‘nt 


and the valuations 
I W30nty S C, II li2y.0 + Cyt" . r || 9, (In; 11.0 ny <C, || 8, n> 


in which constants Cj (j = 1, 3) are positive and are dependent only on 


7 


the coefficients of equations (1) and (2), and V = (v', v2) is a 
velocity vector. 


The answer to the question on the behavior of the solution of system 
(1)-(5) at t - » yields the following confirmation. 


Theorem 2. The generalized solution {#, @} of the quasi-linear 
non-stationary problem (1)-(5) is asymptotically stable at W!2)9.9(m) x 
" Lo (fl), 1.e., 


Cf 
Olly <C, | 8, le 9 


C, (Cy — 2Cs) —2C,t Cal 
et 4 _Cr( 7 O12, (e- * —e *), 





I) wb (205 2¢ <S Coll Po Win) 2 


where the positive constants Cj (j = 4, 8) are dependent only on the 
coefficients of eauations (1) and (2). 
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4. The methods of [1-3] based on quantification with respect to t with 
subsequent linearization and splitting of equations (1) and (2) 
according to the following scheme are effective upon numerical solution 
of the problem of thermal convection of the liquid in gyroscopes: 


Dap, — Drpr—1 — t [7 (pr_1, 7 ?Dp,)/0 (r, 2) + vD*p,| = — tryd0,/dr; (8) 


o (8, — Op—1) — TIO (pe-1, ,)/9 (7, 2) + AY7O,] = 0, (9) 


where fr is a time step, k = 1, m is the number of the time layer, and 
tk =kt, 8, = O(t,), Py = P(ln)- 


Let us apply the Galerkin method with a finite element Hermitian basis 
on triangulation [', of domain fl to system (8) and (9). The approximate 


solution of #> and @h in domains x'=[Je and II'=|[Je will be found 
ecn ecll 

in the form tp"(r, z, t)= N(r, z) p(t), @"(r, z, t)=N(r,2)O(/), Where 

p(O=(p,, Po -- Wm)» Q(t) = {O,, Og): Om,}. are the desired values of the 

corresponding functions and of their derivatives at the nodal points of 

of domains 7h and lh, m,;=m,(C,) (i= 1.2), e is a triangular element of 

domain fl, and N(r,2)={N,...,Nm} is a vector of basic functions [2, 3]. 

Using the known procedure [1-3] of the finite element method in 

combination with the Galerkin method and taking into account conditions 


(3)-(5), we find a system of linear algebraic equations of type Ax = B, 
where A is a square matrix with band structure, and B is a vector, and 


the desired values of vectors #(t) or B(t) emerge as the unknown vector 
x. The given procedure can be generalized rather easily to solution of 
other initial boundary value problems for parabolic-type equations. 
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Tsisarzh, post-graduate student] 


[Text] An effective solution of the problem of improving the accuracy 
in design of follow-up systems is to use combination control together 
with methods of invariance theory. The problem of designing an optimal 
follow-up system with combination control in the presence of random 
noise, based on the design system developed in [1], is solved in this 
article. 


Let the block diagram of the follow-up system with combination control 
have the form shown in the figure. The system includes the controlled 
object with transfer function Wo(s), a regulator W,(s), master signal 
meter q(t) with transfer function W(s) and filter Wa (s). The problem of 


optimal design includes determination of the transfer functions of the 
regulator W,(s) and of the filter Wg (s), which supply a minimum of the 


following functional in the class of stable closed systems 
1 =r (e?(t)) + c(U? (d) (1) 


at given Wo(s) and W(s), at known statistical characteristics--noise 
R{t) and n(t) and with arbitrary variation of the master signal q(t). 








Let us write the equations for a Laplace error transform for the 
combination follow-up system and control signal according to the block 
diagram 


e(s) = [1 — F, (s) — F(s) Wp (s)] 9(s) + Wo(s) Fa(s)n(s) + 
+[1—F,(Q) RS), (2) 


U (s) = [Fi (s) Wo (s) + Fa(S)) 9 (8) + Fa(s)n(s) + F, (9) Wo! (8)R (5), 


where 


F, (s) = W, (s) Wo(s) [1 + Wy (8) W, (8); (3) 


F,(s) = W(s) Wo (s) [1 + Wi (s) Wo(9I™. 


As is known [2], the transfer function of the closed system should be 
equal to unity to guarantee the absolute invariance of the system with 
respect to the master effect q(t). However, it is impossible to realize 
an infinite bandwidth in real systems; therefore, let us require that 
the following equality is fulfilled in the system to be designed 


F, (8) + Fy(s) Wo(s) = E (s), (4) 


where E(s) is some transfer function, close to unity in the proposed 
range of variation frequencies of the master effect. 


If the following ratio is added to equality (4) according to [1] 


a (s) F, (s) + B(s) Fy (s) = ® (s), (5) 


where a(s) and §(s) are some polynomials of s, the system of equations 
(4) and (5) permits one to express transfer functions F;(s) and Fo(s) 
and functional (1) to be minimized by a single function 9(s): 


> pe) . LE (s) B(s) — Wo (s) D(s)} 
Fi() = iB) at) Wi ; 


_ 19 (5) —a(s) £(9) 
*s9)= “Beja Wel’ al 











+joo 2 
a= £. D (s) — a (5) E (5) 
J=/,+/,4+ j \ || Beyer We W, (s) S, (S) + 
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—6 Bi) — (5) ¥, (5) 
+)e0 2 
+ > { | E (s) B(s) — Wo (s) ® (s) Wo" (s)| Sr (s) ‘ 


+|1— E (s) B(s) — Wo (s) ® (s) 





Sr (s)| ds + 





B (s) — a (s) Wp (s) 


—jo 








O(s)—a(s) E(s) |? 
+| B (s) — a (s) We (s) Sn 9] ds, 
where 
r He +je 
—joo te 


are components that are dependent only on function E(s) and the spectral 
density of the master signal S,(s) and Sp(s) and Spj(s) are the spectral 


noise densities of the error and master signal] meters. 


Function $(s), which transforms the first variation of functional (6) to 
zero and which has bands only in the left half-plane s, is determined by 
the relation 


D(s) = [By(s) + By (s)] D' (s). (7) 
Here D(s) is a function determined by factorization of the expression 


[rW,(s) W,(—s) +c] [A(s) A(— 5)I™ [S, (8) + Sa(s)] = D(s)D(—3); (8) 


Bo(s) is an entire part (polynomial of s), B,(s) is a proper fraction, 
having bands only in the left half-plane s in the expansion 


T (s) D~' (s) = By (s) + By. (s) + B_(s); 


T (s) = E (8) [A(s)A(—s)]7" [Wo (5) Wo (— 8) +c] [ax (8) Sp (8) + (9) 
+ B(s)Wo' (s) Sr(s)] —rW,(—s) A (— 5) Sr (5); 
A (s) = B(s) —@(s) Wo (5): (10) 
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A minimum value of functional (6) is achieved 


+joo 
Imin = 1, +1 ++ ( {B_ (6) B_(—s) + 
—je0 
E (s) E (—s) [c 4- rWo (s) Wo (— 5)] 
* [Sp (s) +S, (S)} Wo (s) Wo(—5) Sr(s) S, (Ss) + 


Sp (8) S, (8) [E() + E(— 9) 
+ 1Sr(s) — Sp) +5, 








r*Wo(s) Wo (—s) Sie (5) 
~ (Spl) +S, (51 (ro As) Wo (— 8) 4-¢] 





ds, s = jw. 


Example. Let it be required to design the structure of a regulator 
W,(s) and of a filter Wz (s) of the control system of an inertial object 


with transfer function \,(s)-=(1-4-7Tos)-', combination control of which is 
achieved by signals from the error meter e(/) = (t)—ac(t)4-R(1) and free 


gyroscope q(t) + n(t). The system should have first-order astaticism 
with respect to the master effect, i.e., there should be no static error 
in the system at constant master signal. This requirement is satisfied 
most simply when the transfer function is selected in the form E(s) = 
= (1 + T,s)>! » Let us represent the noise of the error meter R(t) in 


the form of a random steady process of the white noise type with 
spectral density Sp = R?. The initial error no and its systematic drift 


n;t will be considered as the gyroscope error, i.e., 


n(t) = ny + nyt. ian? 


The spectral representation of the determinant process (11), according 
to [3], has the form 











Sa(0) = lin lstay + OH aR | [ss + the }- 


(nj — ngs?) 
(— s)? (s)? 





Since the solution of the design problem is independent of the type of 
polynomials a(s) and J(s) [1], it is convenient to assume that a@(s) = 0 
and §(s) = 1, and then A(s) = 1. Substituting the necessary data into 
(7)-(10), we find 
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D(s) = (6 + ys) (Rs* + ys + m4) (1 + Tys)7's~?, 


Here 


B=(C+r); y=oTy y=(2Rn, +m); T (s)=RE(C— 
—1rT 9s — y*s*) (1 — Ts) ' (1 + Tes); 
By (s) + By (s) + B_(s) = T (s)D(—s)' = 


R? (c — rT 9s — y*s*) (—s)8 
(1+ 7s) (8 — ys) (Rs* — ys +1) 


A A (Ass+Ay) 
= By + >; * + (Rs? —ys +m) ’ 


i — 
(+75) | W— ys) 
_ —1 (bp +- 6,5) (1 + Tos) s? 
0) = 1808) + Be NOY” = TE TNO + w) RE + FM) 











where 4,=Ry+A,Te; 6,=RyT,; Ay = R* OTe—y) Te! (TE + yTe + RY. 


The final purpose is to determine the transfer functions of the 
regulator W,(s) and filter Wa ( s), which are related to the introduced 


function $(s) by the relations 


LE (s) B (s) — Wo (s) ® (s)] | 
W1(S) = WHIP 1 — Ets + We (8) O(N) — 2) a 


[D (s) — a (s) E (s)] 
Wo (S) = We iW, (1D (8) — 2 9) + BUS) 11 — E(5))) 








Substituting the necessary data into (12), we write 


s (1 -++ To) (69 4- 515) . 
T?. (0 + ys) (Rs* + ys + ) + 5 (bo + 615) 


(1 -+ Tos) [fys* + fx8 + fol 
Ws (5) = SUF yey (RSF us my) +8 (bo FOS) 





Wo (s) = 





where f, = 7 (RO + yy-—A,)— Ry; . fy = Te (yO + ym); fo = Te OM. 


The minimum value of the variance of error of the designed follow-up 
system will occur in the absence of constraints on the mean power of the 
control signal, at c = 0, r = 1, and will comprise 
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+Joo 
(€*()min = Ly + == | {BX (8) B_(—s) + Sa(s) Sq (8) (Sais) + 


—joo 


+ Sy (97 (1 + E (8) E(— 8) — E (3) —E(—s)]} ds, 5 = jo. 


Substituting the necessary data into (13), after integration, we find 





(€?(f))min = I, + 





RT | (y +-,T,) no + n?RT, 





y (R + 97g + 1,T,)* 
4 Te ly (y + 247) — Rey}* + ny (y + 147 o)* 
(R +97, + T,)* 


The relations found in the article permit one to design physically 
realizable structures and parameters of the components of combination 
follow-up systems that guarantee the maximum achievable accuracy in the 
class of linear stationary systems with combination control. 
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Polytechnical Institute] 


[Text] Let us consider a controlled object with perturbations acting on 
it, described by a linear matrix differential equation of type 


x (1) = Ax() + Bu(t) + F, (1) 


x (t,) = Xo, x€ R", u€ R”, 


determined in the open domain U C R™, A and B in (1) are constant n « n 
and n x m matrices that characterize the dynamics of the object and the 
effectiveness of the control members, respectively, x(t) is a vector of 
the phase state with components x;(t), ..., Xn(t), u(t) is a control 
vector with components u;(t), ..., Um(t), F is a n-dimensional vector of 
constant perturbations, and xy is a n-dimensional vector of the initial 
state. 








A block diagram of the system, corresponding to equations (1), is 
presented in the figure. 
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Solution of the problem of designing a regulator consists in 
determination of the control law u(t), which establishes the values of 
the control actions as a function of the phase vector [1]. Let us use 
the canonical representation of initial system (1) to design this 
regulator, as was done in [2] for the special case. 


Design. Let us assume that the phase vector x(t) is measured and that 


system (1) is controllable. After substitution of the variables 
according to the formula 


z(t) = Px(t), (2) 


where P is a constant n « n matrix, system (1) can be written in the 
following canonical form: 


where Ao is am * m matrix, Aoo is am * n-m matrix, E; and Eo are 
unjit n - m * n - m and m « m matrices, respectively, and z;(t) and 
z9(t) are vector columns of dimension m, n - m. 


0 &£, |z, (¢) | 0 || IF, 
a) + | 


2, (t) 


2, (t) 


4 (3) 


’ 





















































As Are Zy (f) Ey F, 


We note that the linear nondegenerate transform of coordinates (2) in 


the space of states of the linear system corresponds to ordinary methods 


of transformation of the block diagrams in automatic control theory. 


find the law of control u(t), let us differentiate system (3) and let us 


substitute the variables y;(/)=z,(f), yo(1) =z: (t), y¥3(/)= zo(t), then 


















































, 0 
vy, () 0 E Ol fy, @ 
Vs (1) 0 Ax, Asell Ilys (4) E 


System (4) still does not yield a simple method of finding the law of 
control. To achieve this, let us introduce the new control vector 


C(y, u) by the formula 


&(t) = Any; (t) + Asay (t) +: r (f)- (5) 
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Let 


E(t) = — Kyy, () — Kaya (t) — Kays), (6) 


and Kj (i = 1, 2, 3) are m x m, mx mand m ~« n - m matrices, 
respectively. 


We find from expression (5) with regard to equality (6) 


u(t) = — Kyy, (Q) —(Ke + Aa) ya) — (Kg + Aca) ¥o (0. (7) 


Having substituted expression (7) into systems (4), we find a closed 
object + reguiator system of type 


0 E 0 y, (4) 






































y; (4) 8) 
Vo (t) || = 0 0 E Vo (f) ||» 
Ys (0) —K, _ Ky — K;|| Ys (f) 


where blocks K; and Ky have dimension m * m, while block K3 has 
dimension m x n-m. Integration of expression (7) in variables zj(t) 
(i = 1, 2) yields 


T 
u(i)=—K, \ z,dt — (K, + Ay) Z, (t) — (K; + Ae») Zz (t). (9) 
0 


Let us represent the matrix of transform P as a block matrix to derive 
the law of control (9) in the initial variables. Expression (2) then 


assumes the form 2; (t) =P1)x; (4) -+P10x2(t), z2(t) = Pox) (1) + PooX0(L). 


Let us substitute the values of variables into expression (9) and, 
reducing similar terms, we find 


; 
\ x,dé 
u(t) =—WKyPuiKsProll |r | -WKe+ Aad Pur + Kat Ad (10) 
| x,d¢ 


0 
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K Pylx,()— (K, + As;) Pas + (Ks + Az) P49] Xs (4), 


or in more compact form 


T 
u(t) = —K,x()— Kg | xdt, (11) 
0 
where 
K, = || Ky + Ag} Ks + Aaa ll P: (12) 
Ky =| KiPui K:Prsll (13) 


P;; is am * m matrix and P;> is am *« n - m matrix. 


The derived law of control of the regulator uses no information about 
the perturbation vector. 


The first term of expression (11) is the law of control of proportional 
regulation and can be found independently of the perturbations acting on 
the object. 


Let us again return to system (4). Our problem with respect to vector 
y(t)==(yi(¢), y(t), ys(t)) includes the fact of guaranteeing y(t) - 0 at 


t ~ » from any initial states y(0). These conditions can be fulfilled 
if pair {A, B} is controllable [1]. The derived law of control permits 
one to change the dynamics of closed system (8), selecting coefficients 


Kj (i = 1, 2, 3) of the law of control of type (6). 


Having substituted (11) into initial system (1), we find a closed 
perturbed system, consisting of an object and regulator: 


T 
x (t) = [A — BK] x(t) — BK, | xdi + F. (14) 


0 


The law of control (11), which guarantees given dynamics to closed 
system (14), can be designed by using canonical representation of the 
initial system in the form of (8). 
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[Text] The continual increase of technical requirements on the accuracy 
of spatial orientation of moving objects makes the problem of working 
out algorithms for formulation of the components of the phase 
orientation vector timely with respect to signals of the discrete 
spatial position sensor, mounted on a movable object, the body of which 
is deformed. This causes the appearance of signals of the sensor with 
amplitude comparable to the permissible orientation error. An algorithm 
is known that permits one to formulate current values of the phase 
orientation coordinates of the body of the object in case of small 
deviations from the required spatial position and in the case of small 
absolute values of angular rotational velocity [1] (the spatial rotation 
of the object about the center of mass is described with sufficient 
accuracy by three mutually independent plane revolutions with respect to 
the principal axes of inertia of the body of the object). This 
algorithm assumes that the acceleration of the object ¢€., €,, caused by 


y B 
the effect of the controlling and perturbing moments, is given. 


Therefore, it is suggested that the parameters fy and ts be identified 
periodically by an algorithm according to the data in [1]. It is 


obvious that interruption of the orientation process, caused by the need 
to identify parameters fy and eR? is insufficient. 


Let us study the effect of the constant error of a priori data on the 
accuracy of formulation of current values of phase coordinates by the 
mentioned algorithm. 


Lh1 








The equations of plane rotations of the object with respect to one of 
its principal axes and the corresponding equations of the filter have 
the form 


X = AX + Be; (1) 


X = AX + Be + K(Z—ExX); sas 


~f{9 IF). p_f%]. s_ 7 . 
A=| 9 e=l' |: e=(M,+M,)/~"; Z=EX +8; 


Xx" = (X,, X,); X = (X,, X;); E= (1, 9); kK’ =(Ai, K,). (3) 


Here My and M, are the controlling and perturbing moments, respectively, 
applied to the object along the control axis, I is the moment of inertia 


of the object with respect to the control axis, X and X are the phase 
vector and estimation of it, respectively, X; and X» are the angle and 
angular rotational velocity of the body of the object with respect to 
the control axis, K is the filter feedback coefficient, § is the noise 
component of the data to be measured, ¢€ is rotational acceleration of 


the object, € = € + S€ is information about the angular acceleration of 
the object to be used in the filter, and fe is the constant error. 


Since the influence of the error §¢€ on the accuracy of estimating the 
phase vector must be determined, let us subsequently assume that the 
noise component of the data to be measured is equal to zero. 
Eliminating Z from equation (2) and subtracting (2) from equation (1), 
we find 


X = (A— KE) X — Bee. (4) 


Since Bée = const, for the steady mode from equation (4) follows 


X = —(A—KE)"'Bbe. (5) 


Finally from (5) we find 


1he 








: ' 7 
Xiyer = K, 98 Xayer = Re 66. (6) 


Thus, the output of filter (2) will contain an estimation error caused 
by inaccuracy of the a priori-given data on the external moment applied 
to the object with respect to the control axis. It follows from 
relations (6) that the error in estimating the phase ccordinates is 
proportional to the data error é¢€ about the acceleration of the object. 
Estimation error (6) can be minimized using a refinement of the current 
value of € with respect to the orientation device [1]. On the other 
hand, due to the inevitable measurement errors, contained in the 
instrument readings, the constant component of estimation error (6) can 
not be completely eliminated within the considered algorithm. This 
deficiency can be corrected in the following manner. Let us write the 
equation of motion in the object in the form 


0 1 0 0 
K=AX4+Be A=f0 014s Ba]a fs X= Xe Xd 1) 
0 0 0. 0 
where X3 = é€. 
In scalar form X, = X,, X,= X42, X,=0. 
The corresponding equation of the filter has the form 
X = AX |- Be 4-K(Z--EX); Z=EX +8 a 


E = (1,0,0); X7 = (X,, Xa Xq)s K7 = (Ky, Kas Ka); B75 == (By 0. 0). 


Using the known method of [2], it is easy to ascertain the total 
observability of system (7) and (8). The equation for it with respect 
to the estimation error, which is found by subtracting (8) from equation 


(7) has the form ~—(4—Kf)X, It is obvious that X = 0 in the steady 
mode, i.e., X;=X), Vo=No, \3= ¢€. Thus, filter (8), unlike filter 
(2), has no constant components of the error in estimation of the phase 


orientation vector. 


It should be noted that actuating members, having static relay 
characteristic and operating in the pulsed mode, are ordinarily used in 
real control systems of movable objects. The perturbing moment is a 
slowly variable time function and can be assumed constant with 
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sufficient accuracy with regard to the length of the accurate 
orientation mode. In this case ¢ is a piecewise constant value. 


Accordingly, the real value of the error of a priori data will prevent 
discontinuities of second kind at moments of variation of state of the 
actuating member. It is obvious that this operating mode of the control 
system causes fluctuations of the observed values of the phase vector 
with respect to its real value. The amplitude of these fluctuations 
will be determined mainly by the inertial properties of the observer, 
while the frequency will be determined by the frequency of variation of 
state of the actuating member. The dynamic characteristics of these 
fluctuations are minimized during design of an observer [1]. The motion 
of the specific control system, into the circuit of which observer (8) 
is introduced, was studied by the mathematical modeling method. 


The model of the control process included models of the orientation 
device, filter, control signal shaping module, actuating member and 
object. The model of the orientation device was programmed in the form 


of operators US =KR*FI, UD=UD*EXR(—H/TP)+US*(1.—EXP(—H/TP)), 


where US and FI are the output and input values of the static 
characteristic of the device, respectively, KR is the characteristic 
slope of the characteristic of the device, UD is the output of the 
orientation device, and H and TR are the integration step and time 
constant of the device. 


The model of the filter was programmed in the form 


X1 = X(1), 
X (1) = X (1) + H« (X (2) + K (1) « (FIF— XI), 
X (2) = X (2) + H*(X (3) + E+ K (2) « (FIF—X])), 
X (3) = X (3) + HK (3) + (FIF — X1), 


where X is the file of the estimates of angle, angular velocity and ée, 
K is the file of the filter feedback coefficients, FIF is the output of 
the orientation device, reduced to the dimensionality of the angle, and 
E = € + é€ is the anticipated value of the angular acceleration of the 
object on a given integration step. 


The model of the control signal shaping module (BFUS) is a mathematical 
description of a three-position relay with hysteresis, the input of 
which is a linear combination of estimates of the angle and angular 
velocity of rotation of the object. The model of the BFUS was 
programmed in the following form: 


SIG = KI « X(1) 4+ K2 « X (2), 
IF (ABS (SIG)-LT-SO) F = 0, 
IF (ABS (SIG)-GT-SS) F = SIGN(1-, SIG), 


14h 








where SO and SS are the start and response thresholds, respectively, of 
the relay, and F is the output of the BFUS. 


A model of the actuating member (of the electric flywheel motor--EMD) 
was represented by operators M = -MY¥* *F+\IS,. MD=MD*EXP(—H/TD)+ 
+ M*(1.—EXP(—H/TD)), where M is the anticipated value of the control 


moment at the given integration step, MY and MS are values of the 
starting moment of the EMD and of the moment of resistance, 
respectively, and MD is the moment of the EMD at a given integration 
step. 


According to equation (1), the model of the object was represented in 
the form of operators 


Fl = Fl + FIT+H + (MD/I +H ++ 2/2., 
FIT = FIT + (MD/I) «H, 


where FIT and FI are the angular velocity and angular rotational 
coordinate of the object, respectively, and I is the moment of inertia 
of the object with respect to the control axis. 


The results of simulation showed that estimates of the components of the 


phase orientation vector X; and Xe fluctuate with respect to the 
estimated values with amplitude not exceeding 15 percent. There is no 
constant error component. The accuracy of orientation of the object, 
controlled by using estimates of the orientation vector, formulated by 
filter (8), and by using ideally accurate values of the phase 
orientation vector, essentially coincide. 
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[Text] Let us consider in this paper the generalization of the Lame 
theorem [1] for the problem of the theory of elasticity of a rotating 
medium. These problems were analyzed earlier only by the perturbation 
method [2, 3] or by other apporximation methods. 


The effect of rotation of an elastic medium as a whole on vibrations is 
refiected in the equations of motion, having connected after d’Alambert 
the inertial loads in a rotating coordinate system: 


(A + 2n) 9 (y-4)— py x(y xu) ~ p(F 4+ 0+ Rx U4 QXu) + 0(Q%, (1) 


where Q= {Q,,Q,,Q2,} is the angular rotational velocity of a body as an 


— 


absolute solid, u—{U,, Uy, u,} is the vector of elastic displacements, F 
is the load vector, and / and » are Lame constants. 


Let us assume that servocoupling in the form of a tor iich 
compensates the effect of elastic vibrations on rotation at angular 
velocity--a given time function--is applied to the body. 


According to Helmholtz’s theorem, representation in the form of the sum 


of the potential and eddy parts of the field u= uy 5 uy, v™ ti == Q, 


V-u, --() is valid, and this representation is unique. It is equivalent 


~ 


a= 


to % = VP, uy = vx, where yg and # are the scalar and vector 
potentials. In similar fashion, F = yl +yA. Substitution of the 
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“., 


expansions into (1) yields (A + 2p) <P (Y-m) — ny X(V xa,). for the left 
side of the equations of motion. We find A=y(y-)—y~x with regard 
to the identity (A + 2m) Au, — pdu, (Vx). 


The terms of the last equation are an expansion of the left side onto 
the potential and eddy components: (4/4 +21) y-(Aqg)— py» (Ayp). 


Let us write the right side of equations (1) in the form of a similar 
expansion. Let us made the transformations on the example of the 


~ ”~ 


component (H * u)x, where Hf — 229, -|- a, 


(H XU) _ Hyu, ~ Huy = (— Hy py — HAp,)sx T 
+ (— H,9 + Hy Px)sy + (Hy T H Wx): 


— 


In view of V-¥ = 0--the constraints which can be imposed on the vector 
potential (it will be determined with accuracy up to a time function), 


we write (ix u), =~ (H,»p,— Hy py— Hp,),-1- (Hy py— 16 1- HyPx)ry “b (Hp, -\- Lyp +- A py)s2- 


We represent the other components of H « u in similar fashion. 
Integrating and omitting the additive time function, we find 


Ag = w+ Hew- 
ceAg = p+ H-p+, (2) 


> — 


c?A) = p + Hxp— Ho +A, 





where c, = a 2H C, = |/ rj is the speed of sound. 





Let the boundary value problem be solved for a fixed medium and let the 
functions aF = Fe be found, which determine the solution = ,,(t)F,, 
p-- pF; at U=}y(t)Fy, A=}() Fy. Then upon rotation 


D, H-p=c + fo, 
Pot p IXPo 0 (3) 


p+Hxp—Hp, =cikp +f 
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or in the form of a quaterion equation 


q+ Hoqg=Aq+é, 


- (4) 
where g= (pop); H = (0H); g = (f/f). 


Let us consider the special case of rotation about a fixed axis. 
Without limiting generality, let us assume that , =2,Q,=,=—0, 


H; = H. Then (3) can be represented in the form of the direct sum of 
the subsystem 


Po — cin) — Hp, = fy, 


Pe— CHP + Hp, = i. 


and of a subsystem similar to a Foucoult pendulum 


Py — cixPy — Hp, = fy 
P,— ciyp, + Hp, =f,. 


In the case of the plane problem of elasticity theory [1] 














E _—— 
(here D =@—Q%; H, = 20,4: +0.0,; 2, =(q—yyp) 7% t= 7S 
and E and v are Young’s modulus of first kind and Poisson’s coefficient) 
similar results can be found with respect to functions = (Us + Uy), 


Y=-> (U,;—u,), if the operations of taking the divergence of the rotor 


as well are applied to (5). Then 
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Ae = De — Hey, 


(6 
(1 —v) Ay = Dy + He, 


and if one assumes that e=p,(i)F,, 7=49,(F,, then from (6) follows 


Dp, — fq, + x, = 9 ve 
Dq, + Hyp, + (1 — v) x9, = 0. 


It is easy to find the solution of (7) similar to Foucoult’s pendulum at 
v = 0 and at arbitrary function {(t). The characteristic equation of 


system (7) has the form (Xx —— Q — w) (1 — v) y—2?— @) — 40%? = 0, at { = 


= const, while the solution is 
p(t) = P eft" +k. c, 9, (t) = Qew" +K. ¢. 


This solution is specifically applicable to the problem of vibrations of 
a disk in its own plane [1]. 


In the case of fluctuations of a thin ductile ring in its own plane, the 
solution of the equations of motion [4] 


(00 +0)+ 1 0(00—1) = Do— He a 


— 0° (dw + v) + — (dv — w) = Dw + H,v 


(here —~ 8. _ pp pp FY . ._ 1 /h\? J and S are the moment of 
2, = ;? t= kt; k? = RS? & = 5 (x) 

inertia and cross-sectional area of the ring, R and h are the radius of 
the mid-line and thickness of the ring, v and w are components of the 
vector of displacements of the points of the mid-line in polar 


coordinates, and @ = -) can be found in the form 


uv = p, (t)cosng + q, (t) sinng, 
(9) 
w = q" (t)cosng + p, (4) sinng. 
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Having substituted (9) into (8) and having omitted subscript n, we find 


Dp' — Hq" + ap’ + ap* = 0, 
Dq' — Hp" + a,q' —ag* = 0, (10) 


Dg" + Hp’ + a,q"— aq’ = 0, 
Dp” + Hq’ 4- agp” + ap’ = 0, 


where a,=n(I+—-); gant+—; a=n(nt+—), 


It is easy to analyze the dependence of the solution of system (10) on 
the small parameter ¢€ by the perturbation method at {) = const [5]. 


Thus, if the classical problem of elasticity theory is solved, the 
problem reduces to analysis of an ordinary quaternion equation (4) upon 
rotation. Therefore, the results permit one to construct 
generalizations of known problems on surface and volumetric waves, and 
vibrations of disks, cylinders, and rings. 
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[Text] A non-autonomous system, described by differential equations of 
the following type, is given 


dx/dt = f (t,x), (1) 


where x=(x,,..,x,)’ is a n-dimensional phase vector, {(/,x)=(f, (t,x), ... 
.,fa(tx))’ is a vector function, and {,(t,x)EC'(RXR",R), i= (1,...,2), r>2:; 
A={|(:4:] is a finite time segment, and to < t; < +o. 


Let us consider along with them the following system of equations: 


h(x) =0, R=(1,...,L) baa; (2a) 


dx,,/dt = f,, (t,x), m=(L+1,...,n). (2b) 
Let x’(t) be a solution of system (2) for 


xin (ty) € Xo (Xo CR" EA), (3) 


while x(t) is the solution of systems (1) for 
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X(f)= x), (fo), m=(14+1,...,0), 


(4) 
Xi(4)€Xo, XoEll, i =(1,..., 2), 


where Il is some set and X.°, C R®° is a set of initial conditions. 


Let us find the conditions at which the following inequality is 
fulfilled for t € A. 


P(x (t), x(t) = Ix (Q—x (|| <e, (5) 


where A.—[/,; fi]; lo<fe<hi; e>0 is a sufficiently small number. 


Definition 1. The vector function x’(t) is called the manifold of 
system (1) at initial conditions (4) if condition (5) is fulfiiled. 


Definition 2. Points (t, x) of n + 1-dimensional space, which form some 
set E®*! and for which condition (5) is valid, is called the attracting 
€-neighborhood of manifold x’(t) of system (1). 

System (1) can be considered as a one-parameter group of transformations 
g':X"t'_. x"t! of phase space X9t!, generated by vector field f(t, x), 
which maps a [translator’s note: one word omitted] set X96 (x; (t,)€ Xo) 
into some finite X®, [translator’s note: part of item missing]. Let 
X®; be measurable according to Jordan [1]. 

Definition 3. Manifold system (1) has g-invariance of finite set X", 
with respect to some initial X%) if: 1) there exists such a moment of 
time t_ € A that all the integral curves emerging from Xjyf>t,, belong 
to Entl, i.e., Xie bt; 2) uXi<v, where # is a n-dimensional Jordan 
measure, v > O is a sufficiently small number, which is a Jordan measure 


of the cross-section of set E®*! by hyperplane T = ((t, x):t - t; = 0). 


Let us introduce the following geometric constructions into 
consideration. Let there be given hypersurface S, = (xq, t,x) Xn — Pr (t,x) = 
= 0), where xk is a derivative of function x, at point (t, x), be given 
in space R®**, Hypersurface Vy, of space R®°*! will be the intersection 
Sk and hyperplane y, = ((x,,{,x):%,=0) + Let us be given the set of 
dimensionality y+ 2Gn*? (x,t, x) Rt! such that git’?]S,# mS. Vr€ ‘we 
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Let Git? (U, =U" (t-%) be some set of dimensionality n + 1. Hyperplane 
Uk divided G"**, into two sets G*k and G~, such that if S, Gj Gt, 
then fx(t,x)>O y(t, x)eMy; if S,natea,, then f, (1, x1.<Oyi(t, x) € 


€ Il-,~, where It, and Il-, are sets by which [®*!; is divided by 
intersection Vx. Set Il*, is determined as an orthogonal projection of 


each point (xq, 4x) E(SaN Ge’) onto Ux, for which x, > 0; Iv, is 
determined in similar fashion, but with the condition x, < 0. 


Let us consider hypersurfaces V*, and Vx such that if x(t) € (V*k, Vk), 
x’(t) € VE, then 


| dx, (t)/dt| = |dx,()/dt| yte A, (6) 


We shall say that point (t, x) is higher than V_ with respect to Ux, if 
its coordinate x, is greater than the corresponding coordinate x’, of 
point (t, x’) € Vx, found by intersection of straight line P with V,, 


where (/,U,) =90, pEP, vu, €V,, (t,x)€P. In similar fashion, t(x) is found 


below V_ with respect to Ux, if x~ < x’~. If all the points of some 
set A lie higher (lower) than Vx, then A is higher (lower) than V, with 
respect to Ux. Let us also introduce the definition of motion of 
integral curve x(t) to some hypersurface B with respect to Ux, if all 
the points (t, x) € x(t) are located above (below) B and in this case 
Xk(t) is a decreasing (increasing) function of t. Let us require that 
V*x lie above Vx, and let V°, lie below Vx, with respect to U,. The 
following existence theorem is valid. 


n+ | 


i l 
Theorem. Let {]} lk" =II, where Il # @ and (¢€4, VV. # SVEN, k=(1,... 


for t € A. Then x(t)€B"t!:; upon fulfillment of the following 
requirements: 

a) if the arbitrary point (t, x) of curve x(t) belongs to IIx, it 
is higher than V*;,, and if (t, x) € H*,, it is lower than V°,: 

b) the inequalities |/x(t,x)|>|dxf/dt|,k=(1,....0) where xk(t) © V; 
are valid for point (t, x) € x(t), located above (below) V*,(V",). 
Proof. Let us consider the condition when (t, x) € x(t) belongs to Il-,. 
According to a), point (t, x) lies above V*,. Let us project curved 


x(t) orthogonally onto U,. We find part of the cylindrical hypersurface 
Ck in R®*!, the generatrices of which are orthogonal to Ux, while the 
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directrix is a projection of x(f) onto Ux. Let C,NVi =x" (0, C.NV, = 
= xk(t). The points belonging to x(f),x+(), x*() differ only by 
coordinate xx, and x,(/)>x}(t)>x*(f). However, dx,/dt =f, (t,x)<0 and, 
accordingly, x(t) moves toward V_, and V*, with respect to U,. Let 


x} (t) — xt (0< 6, yle A, where 6, is some positive number. Let us consider 


the difference 


t 
Xp — xt = x, (t') — xk (') + ( (fy (t,x) — dxt/dt) dt, ty € ("5 tL. 


ad 


According to a) and b): x, (t')—xt(')>0, f, (7, x) — dxt/dt <0, while 
X, (t)—x}+()>0. Therefore, either x(t’,) € Vt, at some t’, or 


Him (tn ()— 8 ())=0. This means that moment of time tk, is found such 


that ya<t<+ 0 |x ()—xt()|<6,, k=(1,..,0. the path of the 

arguments is the same for (t, x)Ex(MENS . Let whe Ae (le = max (4, ... , 
15) | xp (t) — x4 (fC) | <8, R=(I,....0. Let us consider the limits 6,+>0, te >t. 
We find x()€V,, k=(1....,) within the limit. Accordingly, the 


equations of subsystem (2a) are satisfied. Equalities (2b) are 
fulfilled automatically due to the fact that x(t) is a solution of 


system (1), while equations (2b) coincide with n - 1 last equations (1). 


Thus , t. = to, x(t) = x°(t) at 6, = 0 and one can write lim x (f) = x’ (1) 


or limw(/,6)=0, where a(t, 6) =x(Q—x’' (A, 6 = (6,,..., 6), fe —t,)'. 


$-+0 
Having selected vector 4 as sufficiently small, we find # with small 
coordinates #\(f) fort, < t ¢ t;. Let us assume that e = || w(/, 5) Ilo. 


where || @(¢, 5) Ilpb = a — )- The following conditions are then 


| Pe n 
fulfilled: |xj(t)—*,(|<e ylede i=(1,....") i.e., x(QER"™' The 


theorem is proved.. 


Let us transform condition b) to a more convenient form. Let the 





Jacobian ts #0 y(t,x)€ll, then according to [2], we find from 
Ir eee X) 
subsystem (2a) Xp=Dn (ly Kyo ee Xh—dy Xhpy 2 Ando k=(I,..., l) 
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The complete differentials are 


ap a 
dx, =—* at + Si ay, (7) 
m 


m=] 


Writing subsystem (2b) by differentials and substituting them into (7), 
we find 


ag _ ¢ 
dx, = (te +} oy (t, x)) dt. 


m=| 


Now condition b) can be written as follows: 


LR GwOI< lila] yi xell, 





n (8) 
. ay, | On, 
where I’, (, x) = =~ +- y Ie | (¢, x). 
ms=|,mgxtk 
The equations of hypersurfaces V*, and V~, have the form 
(9) 


fn (t,x) | =|1,(4x)|. 


The following important corollary follows from the theorem. In order 
for non-autonomous system (1) to have manifold properties with respect 
to x’(t) according to k equations, it is necessary that 


of} (t, x) | 


Ox! 


| <0, j=2r+1,r=(0,1,2,...), R=(1,...,0. (10) 
k x=x’(l) 





For proof, let us expand function f,(t, x) in the neighborhood x’ (t) 
into a series. Let us gather together the terms by odd powers x, - 
- x’,(t) and taking into account that fx(t, x’(t)) = 0, we find 


y af 
i, (t,x) = y — 


j=-2r+1 


(x, — x4 (t))’ + R(x —x’ (O), 





x=x’ (ft) 


r = (0, 1, 2,...), 
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where R(x - x’(t)) are higher-order terms. Hence, condition (10) also 
follows. 


It is obvious that contraction x(t) toward the manifold is reduced, 
since only in condition (8) the sign of the inequality is replaced by an 
equality. Since f,(t, x) approaches zero along the manifold, then the 
greater f,(t, x) is, the narrower the attracting ¢-neighborhood is. 


Intersection of hyperplane T = ((¢,x):4—i/,=0) and of set E®*! forms a 


closed sphere B"9(x’(t), €) [3], measurable according to Jordan. Let 

p#B2 = § and then one can assume that X9, = B® and, accordingly, px, = 

= v, Where v = #. The property of p-invariance for systems having 
manifolds can be used in design of controlled dynamic systems, invariant 
to initial perturbations. 





Example. 
dd R -+- CY -0,5pv?S (g 1 . 
“dh mv? sin 8 ~ VV2 Rs +h cte 8; (11) 
dm _ B 
dh vsin®’ (12) 
where 





a R —Cxo-0,5pv°S — (v 4- g/v) mv sin 0 \+ 
0,5R 4- A-0,5p07S | 


Table 1 














0 100,0000 0,3190 
12 99,9226 0,2315 
4 99,8357 0,1898 
16 99,7483 0,1853 
18 99,6678 0,1991 
10 99,5914 0,2088 
»1 99,5510 0,1769 
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Table 2 





h | 2f,/00)p=2 fh) 





2 —49,17 
d —48,49 
6 —45,10 
8 —41,35 
0 —39 , 83 
I —46,78 


Table 3 














0 100,0000 0,3252 
2 99,9290 0.2301 
4 99,8348 0.1894 
6 99,7474 0.1856 
8 99,6670 0,1994 
: 99.5900 0,2085 


1 
l 
l 
l 
l 
2 
2 99,5500 0,175) 


Functions (R, p, v, Cxo, C%y, A, g,...) are complex functions of h; 
therefore, spline interpolation was carried out. The initial system was 
first reduced to dimensionless form and "0100, 606[0,03; 0,32] was 
integrated at h; = 1, ho = 2.1. The results of integration with initial 
condition with respect to #9, equal to 0.319, are given in Table 1. 


To determine the manifold, let us set f;(h, m, @) equal to O and let us 
find 


~ a, Ng’ — (N3m'? — 4Ngm’? (Ny — Ngm’2))'/? ; 
=—* Nam’? (13) 





where 


(R + C%.0,5pv2S)? 
N,= “4 = (R — Cxo-0,5pu?S); 








Ny 


v3 v |’ 


R + CZ.0,5pv2S)2 , . 
_ Sf + Cy me (+ £): 





Ns = (fr — Regn) (0 SApv*S + 0,52). 
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Having substituted (13) into (12), we find the equation for the manifold 
(@’(h), m’(h)). Using condition (10), we write 





0 
ai} 9=6'(hy) <0. 


m=m’ (h) 


The results of estimating the order of this value are presented in 
Table 2. 


Let us compile Table 3 for the manifold. As one can see from Tables 1 
and 3, the difference in absolute value between variables m and m’, @ 
and @’ does not exceed 0.001. 


The absolute difference @, m at point h = 2.1 did not exceed 0.0075 upon 
variation of @> in therange of [0.03; 0.32]. Thus, the convergence to 
the manifold is also good. 
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Tsisarzh, post-graduate student, and V. M. Rechkin, post-graduate 
student, Kiev Polytechnical Institute] 


[Text] The purpose of this article is to study the effect of the delay 
T on the maximum bandwidth of a digital follow-up system (TsSS). 


To solve the postulated problem, let us introduce into consideration the 
discrete transfer function of the error of the digital follow-up system, 
which, according to the expressions found in [1], has the form 


Wo (2, 0) = go (2,0) g7* (2, 0) = [1 + W5(2, 0)"; (1) 


W*(z,0) =(1—27')D(2)Z[W,, (p)e I. (2) 


Here qy(z, 0) and qo(z, 0) are the corresponding z-transformations of 
the control signal; and errors of the digital follow-up system, Wo(z, 0) 
are the transfer function (digital follow-up system) according to the 


error signal at moments of closing the pulsed element, W plz, 0) is the 
transfer function of an open control circuit, D(z) is the discrete 
transfer function of the control program, Z is a symbol of 
z-transformation of the reduced continuous part W(P) = Wo(p)p°!, and 


Wo(p) is the transfer function of the continuous control object of the 
digital follow-up system. 


Let us use the method of the modified z-transformation when accounting 
for the delay 7 in the digital follow-up system [2]. This method, with 
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quantification time Tc > (rt + 6), yields the following expression for 
the transfer function of an open control circuit: 


W3 (z, 0) = (z— 1)27°W* (2, e), (3) 


where 


W*(z,e)=Zel[Wulp)ls &= 870: (4) 


and Zz, is a symbol of modified z-transformation. Let us estimate the 


bandwidth of the digital follow-up system at known delay 17, having 
considered the case of a digital follow-up system, of sufficient 


interest for practice, with continuous control object Wo(p)=(1+7:p)-! 


and with physically realizable control algorithm D(z), the transfer 
function of which, when using a bilinear operator z = (1 + w)(1 - w), is 
represented in the form 


b (1+ Tew) (1+ Tow) 
Dw) = ——e+T) (5) 





where T3, T,, and Ts are the time constants of the control algorithm. 
The transfer function of open control circuit (3) of the digital 
follow-up system with regard to (4) and (5) will then be described by 
the following expression: 


(1 —w) (1 + Tyo) (1+ Tew) (1+ Tye) , 
w (1+ ©) (1+ Ty) (T+ Tw) 





WV, (w) = (6) 


where 7; =(1+d—2d)(l—d)"';  T,=(144)(1—d)"'; d = exn(—a); a=T,77' 
€ = 6T"!59. 


One can show that the transfer function of type 


k(1 —w)-(1 + 7,@) 


VW; (w) = oUt Tow) (7) 





which can be found in a digital follow-up system by selecting the time 
constants of the control algorithm (T, = 1; To = Ts), will be optimal. 
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It follows from analysis of a stable digital follow-up system with 
transfer function (7) that the maximum permissible cutoff frequency will 
be determined by the expression 


2 1 (e%” —1) 


Where 0= 6t7'; y= q77"', 
The maximum permissible cutoff frequency in the digital follow-up system 
is wep (0) =2r-! arctge > at & = 0 (To = rT). For convenience of further 


analysis , let us introduce into consideration the proportionality 
constant m ~@cp(6)@p(0)-1, which, with regard to (8), will have the 


form 





l 
i 1 ce” — 1) |arct ly]. 


(1—e~”) 


Formula (9) permits one to reach conclusions about the considerable 
influence of the time constant T; of the control object Wo(p) on the 
proportionality constant with delay r in the system. The dependence of 
(9) on parameters ¢ and v is shown in the figure. 











p25 
m 
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The maximum permissible cutoff frequency of the digital follow-up system 
(at given delay rT) can be increased with an increase of quantification 
time To by a small value é in a digital follow-up system in some 
practical situations. 


1. If the dynmamics of the control object Wo(p) is characterized by a 
"small" time constant (v } 1) (considerably less than the delay in the 
system), the bandwidth of the digital follow-up system can be increased 
due to the fact that the proportionality constant (9) will have an 
extreme value with respect to parameter g. The lower the inertia of the 
control object (T; € 7), the less the value é should be. 


2. If the inertia of the control object is large (T; } 7), an increase 
of the quantification time To by value é (To = 7 + 4) does not result in 
an increase of the cutoff frequency of the digital follow-up system. 


The relations found in the article permit one to design digital control 
algorithms for follow-up systems (with known delay in them) and to 
determine the optimal quantification time To, at which the maximum 
transmission frequency of the digital follow-up systems is provided. 
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{sb} [Abstracts of articles appearing in “Mechanics of Gyroscopic 
Systems,’ Issue 8, 1989] 


{txt} [Text ] 
{ud}UDC 629.12.053.11 


On Dynamics of One-Rotor Gyrocompass With Tuned Rotor Gyroscope. V. V. 
Avrutov, M. A. Pavlovskiy, and L. M. Ryzhkov, p 3 


The dynamics is studied and the accuracy of a one-rotor correctable 
gyrocompass, based on a two-axis indicating stabilized platform and 
two-race tuned rotor gyroscope (DNG), is analyzed. Separation of the 
complete model of the device into compass and vibration permitted us to 
determine the characteristic features of the dynamics of a gyrocompass 
with DNG. Analysis of the compass motion showed the dependence of the 
instrument error on the drift, time constant and residual stiffness of 
the DNG, and also on the static errors of the follow-up systems. The 
two-axis suspension of the platform led to differences of the 
expressions of intercardinal deviations of the studied device and of a 
one-rotor gyrocompass with rigid torsion bar suspension. 1 Figure, 4 
references. 


{ud}UDC 531.383 


Selection of Geometric Dimensions of Flexible Elements of Gimbal 
Suspension of Tuned Rotor Gyroscopes. I. V. Balabanov and A. N. 
Motornyy, p 10 


A method of selecting the geometric dimensions of the flexible elements 
of a tuned rotor gyroscope upon fulfillment of impact and cyclic 

strength conditions and stiffness requirements on the flexible elements 
is presented. The comparative characteristic of the DNG with flexible 
elements of variable and constant cross-section is given. 3 References. 
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{ud}UDC 531.383 


Study of Dynamics of Two-Degree Gyrotachometer With Rotating Gimbal 
Suspension. V. S. Yevgenyev, A. V. Yeroshenko, and S. G. Bublik, p 14 


The dynamics of a two-degree gyrotachometer with induced rotation, 
mounted on a uniformly rotating base, with different relations of 
angular velocities, is studied. It is shown that this device has the 
properties of a one-race rotary vibratory gyroscope (RVG). The 
conditions of dynamic tuning of the gyrotachometer are found and its 
advantage is compared to the RVG are indicated. 2 Figures, 4 
references. 


{ud}UDC 531.383 


Dynamic Effects in Gyroscope Mounted on Flexible Suspension Upon 
Acceleration of Its Rotor. A. V. Zbrutskiy, L. Yu. Akinfiyeva, S. Ya. 
Svistunov, and I. B. Kushnirenko, p 17 


The motion of a gyroscope with flexible suspension on a rotating base 
during acceleration of its rotor is considered. Separation of complex 
motion into "fast" and ‘slow’ with subsequent use of the averaging 
operation is suggested to study the dynamics of the gyroscope. A 
solution of a system of differential equations of motion of a rotor on a 
flexible suspension, which describes the motion of the rotor over the 
entire range from acceleration to reaching the established mode, is 
constructed. 1 Figure, 8 references. 


{ud}UDC 531.383 


Errors of Angular-Rate Sensor on Tuned Rotor Gyroscope. A. V. Zbrutskiy 
and S. A. Shakhov, p 24 


A mathematical model of the errors of an angular-rate sensor on a tuned 
rotor gyroscope with two-race symmetrical flexible suspension is found. 
The effect of errors in manufacture and errors of the feedback circuit 

on the accuracy of the readings of the instrument. 4 References. 


{ud}UDC 531.768 


Error of Pendulous Compensation Accelerometer With Flexible Suspension 
During Operation Under nee Vibration Conditions. A. M. Ionin and \V. 
M. Slyusar, p 29 


The mechanics of the occurrence of the constant component of the error a 
pendulous compensation accelerometer (MKA) with flexible suspension, 
exposed to spatial vibration conditions, is studied. Expressions are 
found that permit numerical estimation of the constant error component 
of the MKA. The requirements on the parameters of one of the 
accelerometer circuits, fulfillment of which permits one to make the MKA 
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invariant to spatial vibrations, are formulated. 2 Figures, 1 
reference. 


{ud}UDC 532.516 


On Calculation of Temperature Field in Floated Suspension. A. S. 
Kireyev and A. N. Mikhaylovskaya, p 34 


Finite element postulation of the problem of determining a 
three-dimensional temperature field is presented. A calculating 


algorithm is described and the results of calculation of several eigen 
forms of temperature distribution in the cylindrical region are 
presented. 1 Figure, 5 references. 


{ud}UDC 531.313:532.58 


Study of Rotation of Solid in Floated Suspension. A. S. hireyev, Yu. V. 
Radysh, and A. I. Yurokin, p 38 


The kinematic model of rotation of a solid in non-canonical geometry in 
a floated suspension, which takes into account the nonlinear nature of 
the equations of motion of a viscous liquid, is worked out. The 
characteristic features of integration of the derived equations are 
considered. The results of numerical calculations are presented. 3 
Figures, 7 references. 


{ud}UDC 621.391 


Using Redundant Information to Estimate Errors of Measuring Angular-Rate 
Converters. A. A. Leonets, p 44 


Expressions are found that permit one to estimate the multiplicative and 
additive components of the errors of measuring angular-rate converters 
during functioning of the measuring system, consisting of two measuring 
converter modules of the same type, mounted on a rigid object. 2 
References. 


{ud}UDC 621.391:531.383 


Bayes Approach to Problems of Estimating Condition of Rotor Systems. V\V. 
Ye. Petrenko and A. N. Belyakov, p 47 


The distribution density of the signal probabilities was found in 
parametric form on the basis of nonlinear dependence of the diagnostic 
signal on defects of the ball-bearing seat of a rotor system. The 
problem of approximating this density is solved. Estimation of the 
unknown parameter of signal density is found by the given sample. 
Approximations for conditional densities of class of suitable and 
unsuitable products are constructed. The threshold of separation of 
classes is found as a result of using the Bayes criterion and the 
probabilities of errors of first kind, of second kind and the average 
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risk in making the decision are calculated. The proposed approach can 
be used to solve problems of nondestructive diagnostics of symmetrical 
rotor systems. 7 References. 


{ud}UDC 531.383 


Dynamics of Surveying Gyrocompass. V. Ye. Petrenko, S. A. Zakharenko, 
and A. Ye. Ponomarenko, p 52 


Approximate analytical functions that permit one to calculate the 
natural frequencies of a surveying gyrocompass during parametric 
perturbation due to imperfection of the ball bearings are found. 2 
Figures, 2 references. 


{ud}UDC 539.30/32 


On Parametric Vibrations of Vibration-Isolated Pendulum. L. M. Ryzhkov, 
p 57 


The vibrations of a pendulum, vibration-insulated in the horizontal 
plane, during vertical vibration of the base are considered. It is 
shown that unstable vibrations may occur upon excitation at frequencies 
which exceed the resonance frequencies of the system. Expressions are 
found that permit one to select the damping coefficients in the pendulum 
and in the vibration-protection system, which guarantee stability of 
fluctuations. 1 Figure, 2 references. 


{ud}UDC 531.383 


Gyroscopic Effect in Surface Acoustic Waves. S. A. Sarapulov and S. P. 
Kisilenko, p 62 


The effect of uniform rotation of a base on the evolution of surface 
acoustic waves was studied. The functions that link the variation of 
frequency and phase of the autogenerators and delay lines or filters on 
surface acoustic waves to angular velocity are established. 5 
References. 


{ud}UDC 513.36:681.3.06 


Computer Study of Drift of Cushioned Gyroscope During Angular Vibration 
of Base. S. Ya. Svistunov and Ye. V. Semikina, p 65 


The results of using the MTT-1 applications program package in study of 
the dynamics of gyroscopes with regard to a vibration-insulation system 
are presented. The behavior of a three-degree cushioned gyroscope 
during angular vibration of the base, the drift of the gyroscope in 
resonance modes, on shock absorbers and with rigid mounting on the 
object are analyzed. The results of analytical studies and the results 
of operation of the package are compared. 1 Table, 2 figures, 6 
references. 
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{ud}UDC 531.383 


On Azimuth Motion of Pendulous Gyrocompass in Exponential Acceleration 
of Its Rotor. VV. N. Fedorov, p 70 


The motion of the sensitive element of a ground-based pendulous 
gyrocompass in the horizontal plane during exponential acceleration of 
its rotor is considered. The law of azimuth motion of the sensitive 
element is found as a solution of the hypergeometric Gauss equation at 
non-zero initial conditions. 1 Figure, 6 references. 


{ud}UDC 517.94518.6 


Design of Discrete Model of Motion of Solid in Flexible Suspension. A. 
S. Apostolyuk and \V. N. Sheludko, p 74 


A procedure for designing a digital model of motion of a solid ona 
flexible suspension, based on the use of the discrete analogue of 
differential Lagrange equations of second kind, is proposed. The model 
permits one to take into account the nonlinear nature of vibrations and 
can be used in real-time simulation of vibration-protection systems. 1 
Figure, 6 references. 


fud}UDC 007.62 


Algorithm for Real-Time Calculation of Reduced Moments of Inertia of 
Multilink Manipulation Systems. S. G. Bublik and V. S. Yevgenyev, p 80 


A method of computer calculation of the reduced moments of inertia of 
the links of anthropomorphic-manipulators using the design of a system 
of four material points, equal moment to arbitrary solids, is outlined. 
Realization of the algorithm on a microcomputer showed its effectiveness 
and the possibility of use in control systems of complex multilink 
manipulation systems. 3 references. 


{ud}UDC 628.517 


Nonlinear Effects in Vibration-Protection Systems With Stiffness 
Correctors. B. I. Genkin, N. I. Nagulin, and A. G. Arkhipov, p 83 


Nonlinear vibrations of systems with stiffness correctors are studied. 
The possibility of vibrations with large amplitudes occurring in these 
systems in the working frequency band is shown. Ways of correcting 
these vibrations are discussed. 2 Figures, 3 references. 


{ud}UDC 517.946.9 


Mathematical Problems of Dynamics of Viscous Liquid in Gyroscopy. G. A. 
Legeyda, p 87 
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Problems of single-valued solvability as a whole and of the asymptotic 
stability of the three-dimensional problem of thermoconvection in the 
presence of axial symmetry are studied. Estimates for generalized 
solution of {¥#, 8} in the corresponding norms are found. The essence of 
the numerical analytical method of solution, based on finite difference 
quantification of the equations over time t and the use of the Galerkin 
method with finite element Hermitian basis with respect to spatial 
variables is outlined briefly. 4 References. 


{ud}UDC 62.502 


Design of Optimal Follow-Up System With Combination Control. V. V. 
Lyakin and \. V. Tsisarzh, p 92 


The problem of designing the transfer functions of elements of a 
follow-up system with combination control according to the master signal 
in the presence of random noise is solved. A formula is found for 
calculation of the variance of the error of the system to be designed. 
An example of designing a combination follow-up system is presented. 1 
Figure, 3 references. 


{ud}UDC 681.51 


Design of Regulator According to Given Characteristic Polynomial for 
Objects With Continuous Perturbations. Yu. V. Morozov, p 96 


Algebraic design of a regulator using the canonical form of representing 
the initial system of linear differential equations that describes the 
motion of the control object with regard to continuous perturbations is 
presented. The synthesized design has the property of isodromy. 1 
Figure, 2 references. 


{ud}UDC 629.7.05(075.8) 


Algorithm for Formulation of Orientation Vector of Moving Object, 
Invariant to Constant Error of A Priori Data. M. A. Pavlovskiy, G. Ye. 
Anupriyenko, and A. N. Klimenko, p 99 


It is shown that the known second-order estimator for the control system 
of orientation of a solid, making small revolutions of the object with 
respect to the principal! axes of inertia (the plane problem), forms 
current values of the phase vector with constant error at constant error 
of a priori data on the mass-size characteristics of the object and of 
the modulus of the moment applied to the body of the object. A method 
of synthesizing the algorithm for a third-order estimator, free of the 
indicated deficiency, is presented. 2 References. 


{ud}UDC 581.383 


General Solution of Problem of Theory of Elasticity of Rotating Medium. 
S. A. Sarapulov, p 103 
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The possibility of generalizing the Lame theorem for the proilem of the 
theory of the elasticity of a rotating medium is shown. Special cases 

of rotation about a fixed axis, the plane problem of elasticity theory, 
and vibrations of a thin ring are considered. 5 References. 


{ud}UDC 531.01:513.83 


Manifold Properties of One Class of Non-Autonomous Systems. N. S. Sivov 
and M. K. Sparavalo, p 107 


A new class of non-autonomous systems, which have special properties in 
a certain range, which are called manifold properties, is considered. 
The existence theorem of systems of this class is proved. The concept 
of the independence of a finite set of the system, characterized by a 
Jordan measure, with respect to some initial concept, is formulated. 
This concept is called p-invariance. 3 Tables, 3 references. 


{ud}UDC 531.768 


Optimization of Quantification Period in Digital Follow-Up System With 
Given Delay. V. M. Slyusar, V. V. Tsisarzh, and V. M. Rechkin, p 112 


The effect of delay in a digital follow-up system on the maximum 
bandwidth is considered. Relations are found that permit one to 
determine at given delay the optimal value of the quantification period 
in a system, at which the maximum cutoff frequency of the digital 
control system is provided. 1 Figure, 2 references. 
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